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^ i ABSTRACT 

Q_i! A novel approach to study electroweak physics at one-loop level in generic 

o ■ SU(2) L x U(1)y theories is introduced. It separates the 1-loop corrections 

into two pieces: process specific ones from vertex and box contributions, and 
universal ones from contributions to the gauge boson propagators. The latter 
are parametrized in terms of four effective form factors e 2 (q 2 ), s 2 (q 2 ), g 2 z {q 2 ) 
and gwiq 2 ) corresponding to the 77, 7Z, ZZ and WW propagators. Under the 
assumption that only the Standard Model contributes to the process specific 
corrections, the magnitudes of the four form factors are determined at q 2 = 
and at q 2 = m 2 z by fitting to all available precision experiments. These values 
are then compared systematically with predictions of SU(2) L x U(l)y theories. 
In all fits a s {m z ) and a(m z ) are treated as external parameters in order to keep 
the interpretation as flexible as possible. The treatment of the electroweak data 
is presented in detail together with the relevant theoretical formulae used to 
interpret the data. No deviation from the Standard Model has been identified. 
Ranges of the top quark and Higgs boson masses are derived as functions of 
a s (m z ) and a{m 2 z ). Also discussed are consequences of the recent precision 
measurement of the left-right asymmetry at SLC as well as the impact of a 
top quark mass and an improved W mass measurement. 
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1 Introduction 

The Standard Model (SM) of the electroweak interactions has been with us for nearly 
two decades. Despite the general belief that it should be an effective theory valid at 
energies below the Fermi scale, so far no unambiguous sign of physics beyond the SM 
has been found nor any clue to the origin of the underlying gauge symmetry breaking 
mechanism. On the other hand, the accuracy of the experiments testing the electroweak 
theory has improved significantly in the past decade both in low energy neutral current 
experiments and in high energy collider experiments on the W and Z boson properties. 
The precision of these experiments has reached the level, where meaningful searches for 
new physics through the investigation of quantum effects can be carried out. 

The effects may be significant, if there are new particles with masses as light as weak 
bosons, or if many new particles contribute constructively, or if there exist new strong 
interactions among them. Even in the absence of such a signal, constraints on certain 
new physics possibilities can be derived and tightened in future precision experiments. 
With this motivation to study electroweak radiative corrections several groups have made 
efforts towards comprehensive and systematic analyses HlHUf . 



In this report a novel approach to confront electroweak data and theory is presented 
with the aim of a systematic look for new physics effects. In the following, the conditions 
imposed on the electroweak analysis scheme are outlined. 

Since it is the aim to search for new physics effects in the electroweak precision data, a 
model-independent framework to analyse the data is required. As both the experimental 
accuracy and the new physics effects looked for are of similar size as the SM radiative 
effects, it is essential to take account of the SM radiative effects as accurately as possible. 



For testing grand unification of the three gauge couplings [pJ|-p3| the fits should be stud- 
ied quantitatively as a function of a s . Furthermore, the level of precision accessible in 
the near future is such that the present uncertainty in the hadronic vacuum polarization 
contribution to the running of the effective QED coupling constant a(q 2 ) severely limits 
the ability to study new physics through quantum effects. In order to assess the effects of 
possible future improvements in the e + e~ hadroproduction experiments at low and inter- 
mediate energies, the consequences of varying a(m 2 z ) should be examined quantitatively. 
During the course of this study, sometimes the published results of earlier theoretical 
analysis could not be reproduced easily. This happened in most cases because not all the 
details of the assumptions and approximations underlying the analysis have been clearly 
stated in the literature. The quantum effects studied are so sensitive to details of the exact 
treatment of higher order effects and to uncertainties in the analysis that equally sensible 
looking assumptions often lead to a significant numerical difference. We therefore make 
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every effort to render the report self-contained so that all our results can be reproduced 
unambiguously. 

In order to comply with all the above requirements, our comprehensive analysis of 
electroweak precision experiments is performed according to the following steps, by sys- 
tematically strengthening the underlying theoretical assumptions. 

1: All electroweak data are expressed in terms of model- independent param- 
eters. 

For the choice of model-independent parameters, we basically follow the strategy of 



ref. |24| for low energy neutral current experiments, and that of the LEP electroweak 
working group || for the Z parameters. In addition, the W boson mass, the fine struc- 
ture constant a and the Fermi coupling constant Gf are used as inputs of the analysis. 
Some of these parameters are directly related to experimental observables up to correc- 
tions due to known physics, such as the external QED bremsstrahlung effects and the 
quark-parton model, and uncertainties in these correction factors are included as part of 
the errors of the experimentally measured parameters. 

2: The mo del- independent parameters are then expressed in terms of the 
pole positions of the W and Z propagators, and the S-matrix elements of four 
external fermions, quarks or leptons, which are approximated as products of 
two external standard V ± A currents and the scalar transition form factors. 

All electroweak precision measurements that have been performed so far can be expressed 
in terms of the S'-matrix elements of quark and lepton processes whose masses are neg- 
ligible compared to the weak boson masses. To an excellent approximation, chirality-flip 
terms in the loop amplitudes can be neglected and the relevant S'-matrix elements can 
be expressed in terms of the scalar product of the standard V ± A currents multiplied 
by transition form factors depending on the flavors and chiralities of the currents as well 
as the momentum transfer of the process under consideration. External QCD and QED 
corrections can hence be applied exactly as in the SM, and electroweak models can be 
confronted with experiment, once the transition form factors are determined in a particu- 
lar model. The dependence of the fit on the QCD parameter a s and quark masses is taken 
into account by introducing appropriate external parameters. Up to this stage, our anal- 
ysis is quite general, as the formulae are valid for any electroweak model respecting the 
flavor and chirality conservation laws of the SM, that is, for all new physics contributions 
which can interfere with the leading SM amplitudes. 

Although one may attempt to constrain these model-independent transition form fac- 
tors directly by experiment, we find it impractical, since the number of independent 
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transition form factors exceeds by far the effective number of degrees of freedom provided 
by precision measurements. Hence, we perform the quantitative comparison of data with 
theory in a more restricted class of models which are minimal extensions from the SM, i.e. 
those models which respect the SU(2) L x U(l)y gauge symmetry broken spontaneously 
down to U(1)em- 

3: The transition form factors are expanded perturbatively in SU(2)l x U(1)y 
gauge couplings, and radiative effects are classified either as the universal 
gauge boson propagator corrections or as the process specific vertex and box 
corrections. The universal propagator correction factors are then parameter- 
ized by four charge form factors, e 2 (q 2 ), s 2 (q 2 ), g 2 z (q 2 ) and ffyy(q 2 ), correspond- 
ing to the 77, jZ, ZZ and WW propagator degrees of freedom. 

The restriction to the electroweak gauge group SU(2)l x U(1)y implies at the tree level 
that all fermions, quarks and leptons, couple to the electroweak gauge bosons universally 
with the same coupling constant as long as they have common SU(2)l x U(1)y quantum 
numbers. This universality of the gauge boson coupling to quarks and leptons can in 
general be violated at the quantum level, because the gauge symmetry breaks sponta- 
neously down to U(1)em- It has been widely recognized, however, that this universality 
of the couplings holds true even at one-loop level in a wider class of models where new 
particles affect the precision experiments only via their effects on the electroweak gauge 
boson propagators [fLHTOH. This class of new physics effects is often called oblique |],£| or 
propagator corrections, or those satisfying generalized universality [|l(|. This concept 
of universality can be generalized to certain vertex corrections with non-standard weak 



boson interactions [TT| • It is also often useful in theories with non-standard vertex and box 
corrections, such as the supersymmetric SM (SUSY-SM), since the propagator corrections 
can be larger than the vertex/box ones: propagator corrections can be significant either 
because of a large multiplicity of contributing particles or by the presence of a relatively 
light new particle. 

When confronting the electroweak theory with experiment, we adopt this distinction 
between new physics contributions to the gauge boson propagators and those to the rest, 
where we allow the most general contributions in the former, whereas we consider only 
the SM contributions to the latter (vertex and box corrections). 

4: By assuming that the well-known SM contributions dominate the process 
specific vertex and box corrections, apart from the Zhjjbj, vertex for which new 
physics contributions are allowed, we determine from precision experiments the 
four universal charge form factors at the typical momentum transfer scales, 
q 2 = and m|. 
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The new physics contributions may either prevent our ability to fit the experimental data 
within our approach, or lead to non-standard values of the fitted four charge form factors 
and the Zb L b L vertex form factor, Sb{q 2 )- At this stage, the whole body of electroweak 
precision data can be expressed in terms of the two weak boson masses m w and m z , and 
these five form factors, that is, the four universal charge form factors and 5b{q 2 )- Although 
the form factors could be determined at any point on the momentum scale q 2 , they are 
actually measured with adequate precision only at two specific q 2 ranges, namely all four 
charge form factors at q 2 = or q 2 <C m|, while s 2 (q 2 ), g 2 z {q 2 ) and Sb(q 2 ) at q 2 = m 2 z . 
Hence, there are just 9 electroweak parameters measured accurately enough to be used 
for testing theories: m w and m z , e 2 (0), s 2 (0), g z (0) and <7^(0), s 2 (m|), g z (m 2 z ) and 
8 b (m\). 

Apart from the vertex form factor 5b(m 2 z ) the remaining 8 parameters characterise 
the universal propagator corrections. On the experimental side, the three quantities a, 
Gp and m z are measured so accurately that it is justifiable to treat them as constants: 
a = 1/137.0359895 and G F = 1.16639 x l(T 5 GeV~ 2 from the PDG listing ||], and 
m z = 91.187 GeV from the LEP results [pq| . Among the 8 universal parameters above, 



e 2 (0) = 4ira and m z are fixed immediately, while Gp fixes the ratio g^ v -(0)/m^ v , once 
we assume the SM dominance of the vertex and box corrections (5q) to the muon decay 
lifetime. Since the gauge boson properties are fixed at tree level by only three parameters 
in general models with the SU(2) x U(l) symmetry broken by a vacuum expectation 
value, the remaining 5 universal parameters serve to test the theory at the quantum level 
(see Table 1). We therefore first determine from precision experiments the 5 parameters, 
s 2 (0), #l(0)> 9w(®)i ^ 2 ( m z) an d #f ( m l) ; together with 6b(m 2 z ), and then confront their 
values with various theoretical predictions. 

In the fit to the Z boson parameters the strong coupling constant a s (m z ) is treated 
as external parameter which can be varied within certain limits. In this way the analysis 
remains transparent and easy to update. The fitted electroweak parameters s 2 (m?7), 
(?|(m|) and <5 2 (m|) are thus presented as parametrizations in a s (see Table 2). 

When the new physics scale is significantly higher than the scale (< J m z ) of precision 
measurements, new physics contributions to the running of the charge form factors can 
be neglected. 

5: By assuming further that the running of the charge form factors between 
q 2 = and q 2 = m 2 z are governed only by SM physics, three universal param- 
eters sensitive to radiative effects can be determined. We adopt a modified 
version of the S, T, U parameters of ref. ]4] by including the SM radiative 
effects as well as new physics contributions. 
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Among the 5 universal parameters, the values of s 2 (0)/a and g z (0) can then be calculated 
from s 2 (m 2 z ) / a(m 2 z ) and gz( m %)> respectively, using SM physics only. There are then 3 
remaining universal free parameters which correspond to the parameters S, T and U of 
ref. [|J, ex, e 2 and e 3 of ref. 0, or the other related triplets of parameters in refs. 0. 
When the scale of new physics that couples to gauge boson propagators is near to the 
weak boson masses, its signal can be identified as an anomalous running of the charge 
form factors. This point has been stressed in refs. [12[ in connection with the possible 
existence of light SUSY particles. It has also been pointed out that when new physics 
effects to the electroweak gauge boson sector are parametrized by the four dimension-six 



operators of ref. [10|, there occurs anomalous running of the charge form factors [11]. The 
triplet parametrizations are then no longer sufficient to account for new physics degrees of 
freedom, and all 5 parameters in Table 1 should be regarded as free. Several alternative 
approaches to the same problem have been proposed in refs. |T2"|,[H], 



Note that in order to obtain the charge form factors from the three known parameters 
a, Gp-i w, z and the radiative parameters S, T and U that are calculable in a given model, 
the effective QED coupling at the Z mass scale, a(m 2 z ) is needed. Its value is calcula- 
ble from a in the SM but suffers from uncertainty in the hadronic vacuum polarization 



contribution [p7H29fl. The effect of this uncertainty on the final results turns out to be 
non- negligible. In order to gauge the effects due to this uncertainty quantitatively, we 
introduce 5 a = l/a(m z ) — 128.72 as external parameter and allow it to vary in the fit (see 
Table 2). It is then straightforward to examine the effects of shifts in the a(m 2 z ) value 
and the impact of future improvements in its measurement. 

In the minimal SM, the three universal parameters S, T, U and the Zb L b L vertex form 
factor Sb{m z ) depend on just two parameters: m t and m H . 

6: Finally, by assuming that no new physics contributes significantly to 
electroweak precision experiments, we can express all the radiative effects in 
terms of the two parameters of the minimal SM, m t and m H . The x 2 curves 
of the global fit are shown as a function of these two parameters, for several 
values of a s {m z ) and a{m 2 z ). 

The preferred range of m t is presented as a function of m H , a s , 5 a , that of m H as a 
function of m t , a s , 5 a , that of function of m t , m H} 5 a , and that of 5 a cLS cL function 

of m t , m H , a s . The chosen value for the parameter 5 a is essential, since it is not well 
constrained by the present precision measurements alone. 

A clear advantage of this approach is that we can test the electroweak theory at three 
qualitatively distinct levels. If we cannot fit all the data at a given q 2 with common form 
factor values, we should either look for new physics affecting the relevant vertex/box 



7 



corrections significantly or else we should introduce new tree level interactions such as 
those induced by an exchange of a new heavy boson, or from new strong interactions 
that bind common constituents of quarks and leptons. If the 'universality' in terms of 
the above four charge form factors holds, but their g 2 -dependence does not agree with 
the expectations of the Standard Model, we may anticipate a new physics scale very near 
to the present experimental limit ||12|| , or effective higher dimensional interactions among 
the gauge bosons []X0|, P^lfl - New physics contributions which decouple at low energies can 
thus be identified by their anomalous running of the charge form factors. If the running of 
the form factors is found to be consistent with the SM, then our approach reduces to the 
standard three parameter analyses or those with three plus one parameter [ 12 . 1_4 



when including the Zb L b L vertex parameter 5b(m 2 z ) as well. Deviation from the SM is still 
possible, since the SM has only two free parameters, rrit and m H . In this case sensitivity 
to those new physics contributions which do not decouple at low energies remains. 

As emphasized at the beginning of this section, we present at all stages of our quan- 
titative analysis the best-fit values of the model parameters, including a parametrization 
of the x 2 goodness of the fit around its minimum as functions of the external parameters 
a s = <y s (m z )y[s and 5 a = l/a(m|) — 128.72 ~ (l/a(m|) — l/«)hadrons + 3.88. One can 
examine consequences of possible future improvements in the measurement of a s |3(| and 
those of hadronic contribution to 5 a by adding to the quoted \ 2 function terms of the 
form [(a s - aj)/(Aa s )] 2 and [(S a - I^)/(A5 a )] 2 . 

The paper is organized as follows. In section 2, we present our formalism in detail. The 
helicity amplitudes are stated for general four-fermion processes in terms of the universal 
charge form factors and process-dependent vertex and box corrections. Definitions of 
the form factors and the S, T, U parameters are given and their SM values are shown. 
Section 3 contains theoretical formulae for the electroweak observables expressed in terms 
of the helicity amplitudes of section 2, with QCD/QED corrections. Numerical predictions 
are also given for wide ranges of the form factor values, and also in the minimal SM. In 
section 4, we present our model-independent parametrizations of all experimental data 
and, confront them with our theoretical predictions. The universal charge form factors 
and 5fe(m|) are determined by assuming SM dominance in the remaining vertex and box 
corrections. Section 5 presents a systematic analysis of the electroweak data by gradually 
tightening the theoretical assumptions. First the running of the charge form factors g z {q 2 ) 
and s 2 (q 2 ) is tested, then the 4-parameter (S, T, U and 5b(m 2 z )) fit to all electroweak data 
is performed by assuming SM running of the charge form factors. Finally, constraints 
on m t and m H are discussed in the SM fit. The total x 2 °f the SM is parametrized in 
terms of mt, m H , a s and 5 a . In section 6, consequences of the new precision measurement 
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of the left-right asymmetry [31] and the impact of a top quark mass measurement are 
considered. Section 7 summarizes our observations. Details of the theoretical formulae 
used are collected in the appendices. In appendix A, we give all the SM radiative correction 
terms completely at one- loop level, and partly at two-loop level for 0(aa s ) terms. They 
are classified into three parts, the propagator corrections, the vertex corrections and 
the box corrections. In appendix B, we discuss the renormalization group improvement 
of the charge form factors and hadronic contributions to the gauge boson propagators. 
Appendix C gives the complete analytical formulae for the S, T, U parameters and the 
ZbiJjL vertex form factor 5b(m z ) in the SM. Here all the known two-loop level corrections 
are included. We also give convenient approximations to the exact formulae. Appendix 



D provides explicit expressions for the A, B, C, D functions |32] that are used to express 
all the one-loop correction factors. 



2 Basic formalism 

2.1 S-matrix elements, weak boson masses, and charge form factors 

All the precision experiments sensitive to electroweak physics at one-loop level so far 
are concerned with processes involving external fermions, that is, leptons or quarks (ex- 
cluding top quarks), whose masses can safely be neglected in the correction terms as 
compared to the weak boson masses. There are the Z boson properties as measured at 
LEP and SLC, the neutral current (NC) processes at low energies (^C rn z ), the measure- 
ments of charged current (CC) processes at low energies and those of the W mass at pp 
colliders. The relevant observables in these processes are then expressed in terms of the 
S'-matrix elements of four external fermions which form a scalar product of two chirality 
conserving currents. All the information on electroweak physics is contained in the scalar 
amplitudes which multiply these current- current products. 

For example, consider the S'-matrix element responsible for the generic 4-fermion NC 
process ij — * ij (or any one of its crossed channels). This includes e + e~ — > ff as well as 
v ii1 ~^ v yfl- The matrix element has the form 

Tn = MijJ, ■ Jj, (2.1) 

where Jf and Jj 1 denote currents without coupling factors, that is, Jf = i[)f^^P a ipf for 
i — fa, where P a — (1 + a , y^)/2 with a = ±1 are the chiral projectors.^] All radiative 
1 We use the chirality index a = +1 for right-handedness and a = — 1 for left-handedness 
throughout the paper; e.g., P+i = P + = Pr and P_i = P_ = Pl for the chiral projectors, 
f + = Jr and /_ = Jl for chiral fermions. 
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effects interfering with the tree-level SM amplitudes can be cast into the above form as 
long as terms of order m 2 /m 2 z in the one-loop amplitudes are neglected (mj denoting the 
external fermion mass). The one-loop corrections then appear in the scalar amplitudes 
Mij which depend on flavor and chirality of the currents and on the invariant momentum 
transfers s and t of the process. 

In neutral current amplitudes, the photonic corrections attached only to the external 
fermion lines are U(l) gauge invariant by themselves |J. Therefore, finite and gauge 
invariant amplitudes can be obtained by excluding all the external photonic corrections. 
We find the following closed form for the generic neutral current amplitude of (|2.1|) 
at one-loop order (see details in Appendix A) : 



M. 



NC 
ij 



QiQj 



e 2 (s) + e 2 (T\ + T{)(s)-te 2 A^(s) 



(Qih j )^- + {hiQ j ) 



s — m| + is^f-8(s) 



x{(/ 3 « - Q l s 2 ){h 3 - Q 3 s 2 )[f z (s) + g 2 z (T\ + r{)(s) - ig 2 z A zz (s 

+(hi - QiS 2 )g 2 z [l 3j (c 2 T j 2 + Ti)(s) + T{(s) - Qjtfis) - s 2 + tA, z (s)) 
+(I 3j - Q 3 s 2 )g 2 \l 3t (c 2 r 2 + r*)(s) + r*(s) - Qi(s 2 (s) -s 2 + iA lZ (s)) } 



+ Bg c {s,t). (2.2) 

Here s is the momentum transfer of the current Jf , and t is the momentum transfer 
between the fermions % and j. The hatted couplings e = gs = g z sc and all the ultraviolet 
singular loop functions are renormalized in the MS scheme, and hence they depend either 
implicitly or explicitly on the unit-of-mass \i. Three of the four charge form factors of 
Table 1, e 2 (s), s 2 (s) and g z (s), appear in the NC amplitudes : 



Re 11^(9 



9l 



1 + ^ReTV T Z (q 2 

s ' ' 

1-Ren5(g 2 ) 



(2.3a) 
(2.3b) 
(2.3c) 



Imaginary parts of the propagator correction factors denoted by A 77 (s), 
A zz (s) are defined as follows : 



A lZ (s) and 



A 77 (g 2 ) = ImTT^g 2 ), 
A, z ( q 2 ) = sclmTP T Z ^ q 2 ), 

A zz ( q 2 )=lmUZ(Q 2 )- lm[l 



■zz , 



(2.4a) 
(2.4b) 

(2.4c) 
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The vertex functions r{ a (s) and the box functions Bf a f (s, t) are process specific. The SM 



contributions to all the two-, three-, and four-point functions in eq. (|2.2f) are calculated 
in the 't Hooft-Feynman gauge. Their explicit forms are found in appendix A. 

The residues of the 7- and Z-poles are separately //-independent and gauge invariant, 
and therefore physical observables. For q 2 = 0, the vanishing of the vertex functions 

lf(0) = 0, lf(0) = 0, (2.5) 

is ensured for all f a by the Abelian and non-Abelian parts of the Ward identities, re- 
spectively. The universal residue of the photon pole gives the square of the unit electric 
charge e 2 (0) = 47ia. 

Likewise, the charged current (CC) process ij — > i'f can be expressed by 



r cc 



M 1 

13 s-m 2 w + is^ti(s 



1 EEfl? A 9w( s )+9 2 \^(a)+T i *{8)-iA ww {8)]}+B™(8,t), (2.6) 



with an appropriate CKM factor Vn'V^, accounting for quark family mixing. The W 
propagator corrections appear in the charge form factor g\y{s) and in the imaginary part 
A ww (s): 



qW) = 9 2 1 - RelC^) , (2.7) 



A W ,,( 9 2 ) = Imll-V) " lmIlT 2 imw) ■ (2.8) 

Factorization of the external photonic corrections does not hold for the charged current 
processes, and hence all the one- loop correction terms are included in eq. ( |2.6|) . Ex- 
plicit forms of the SM contributions to the propagator function U^ W (q 2 ) are found in 
Appendix A. 

The gauge boson two-point functions appearing in eqs. ( |2.3| ), Q2.4Q , (|2.7|) and (|2.8| ) are 
defined as follows : 

n-tf) = nS^|!W), (2 . 9) 

q ~ In v 

where m v is the physical mass of the gauge boson V (that is, m w , m z and m 7 with 
m 7 = 0) and the subscript T stands for the transverse part of the vacuum polarization 
tensor U^q), 

n£?(g) = (-^ + ^)n^(g 2 ) + (g 2 ) . (2.10) 

Contributions from the longitudinal part of the gauge boson propagators are consistently 
neglected in the one- loop corrections, because they give terms of order m 2 /my (V = Z 
or W) in the weak amplitudes. 
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The gauge boson propagators are calculated in the 't Hooft-Feynman gauge and the 
so-called pinch terms |2|, of the vertex functions arising from diagrams with the 



weak boson self-couplings are included in the overlined functions H^ B (q 



-,2\ 





= n?V) - 


n?V) 


= n? (a 2 ) - 


nf(g 2 ) 


= nf V) - 


IT T (g 2 ) 


= I^ W {q 2 ) - 



e 



2 



2 q B (q ;%,%), (2.11a) 
^(? 2 -f)5o(9 2 ;%,%), (2.11b) 



2 s4 



(q -m z )B (q ;m W} m w ) } (2.11c) 



4vr 2 

# ' J2 ™2 Ua2d ^„2._ _ \ , »2n /„2, 



4vr 2 



(g -m w )[c B (q ;m w ,m z ) + s B (q ; m w , m 7 )].(2.11d) 



Here -Bo is a Passarino-Veltman function |32| in the notation of appendix D. The overlines 
on the vertex functions T^s) in eqs. ( |2.2| ) and ( |2.6| ) and V" 1 (s) in eq. (|2.6|) indicate the 



subtraction of the pinch term associated with this prescription (note, the pinch terms in 
eqs. ( [2. 11|) have a negative sign in our convention). 

The absorption of the above q 2 dependent propagator-like parts of the vertex functions 
into the effective charges improves over the usual method of absorbing the relevant 
vertex term at zero momentum transfer [H] in two ways. One is that the remaining vertex 
parts do no longer give rise to large logarithms of the type ln(— q 2 /m 2 v ) at \q 2 \ 3> m^, 
and hence the effective charges are useful in making the improved Born approximation 
even at very high energies. The second is that the effective charges are now gauge 
invariant and hence their properties can be discussed independently of the other 

process specific corrections of the same order. Most importantly, we can obtain explicitly 
renormalization group invariant relations between the MS couplings and the effective 
charges 

1 1 



e 2 (q 2 ) e 2 (/i) 



1 + ReTT^g 2 ) , (2.12a) 



s 2 (g 2 ) = s 2 (») + ^^-ReTT T ;(g 2 ), (2.12b) 

within the 't Hooft-Feynman gauge of the electroweak theory. This enables us to discuss 
the renormalization group improvement of the above two effective charg whole, 
that is, without separating the contributions from the SM fermions and the rest. The 
trajectories of all the MS couplings (e = gs = g z sc) are completely fixed by the above two 
equations at one-loop level, which can be used to study quantitatively the heavy particle 
threshold corrections in GUT theories ]2T| . 

In the analysis presented here the MS couplings act as the expansion parameters of the 
perturbation series, since we find them the most convenient when studying consequences 
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of various theoretical models beyond the SM. Their usefulness in the SM analysis has been 
emphasized in ref . , and they are often used in the analysis of new physics contributions 
to the precision experiments ||. However, it is not convenient to use the MS couplings at a 
specific unit-of-mass (/x) scale, such as fi = m z , when dealing with a theory with particles 
much heavier than the weak bosons because of the appearance of large logarithms of their 
masses. Hence, we adopt the following renormalization conditions 

e 2 = e 2 (m|), s 2 = s 2 (m 2 z ), (2.13) 

consistently for all processes studied. The above conditions renormalize all the loga- 
rithms of large masses with the help of the renormalization group identities ( 2.12j ) at 



\q 2 \ < 0(m 2 z ). Note that the running of e 2 (q 2 ) and s 2 (q 2 ) at low energies arises from 
the QED x QCD interactions PE| , and hence the ratio e 2 (q 2 )/s 2 (q 2 ) is not an appropriate 
expansion parameter of the weak corrections at \q 2 \ <C m 2 z . Note further that, apart from 
details concerning the higher order terms, the effective charges e 2 (q 2 ) and s 2 (q 2 ) ( 2.12| ) are 
the same as the real parts of the corresponding star-scheme f|| charges, e 2 (q 2 ) and s 2 (g 2 ), 
respectively. More details on the treatment of the renormalization group improvement 
and the hadronic contributions to the charge form factors are given in appendix B. 

Since we adopt the LEP convention |§ regarding mass and width (m v and IV) for both 



Z and W, the Breit-Wigner propagator factors in eqs. (|2.2| ) and fl2.6|) have the running 



width factor, and the imaginary parts fl2.4cj) and ( |2.8|) have the associated subtraction 



terms. These masses and widths can also be defined in terms of the more conventional 



pole masses and widths |37[] , denoted by and T v , p , as follows 



m v = m 2 Vp + T 2 Vp , (2.14a) 



v 



T v , p Jl + (T v Jm v , p y. (2.14b) 



The Breit-Wigner propagator function with the fixed width and that with the running 
width are then related by the exact relation [BS 



1 l+iT v /m x 



m Vp + im VtP T V)P s-m v + isT v /m v 



(2.15) 



The imaginary part of the numerator A zz (q 2 ) ( |2.4c| ) and A ww (q 2 ) ( |2.8| ) are arranged 
such that the imaginary parts of the full amplitudes vanish exactly at zero momentum 
transfer: Ayv r (0) = 0. The theta function 9(s) (0(s) = 1 for s > and 9(s) = for 
s < 0) in the running width factor of eqs. ( |2.2|) and ( |2.6|) then ensures the reality of the 
amplitudes at s < 0. It should be noted that the imaginary part Ayy(q 2 ) vanishes at 
q 2 = my at one-loop level, if all the contributing particle masses can be neglected. As long 
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as the relations ( [2. 14 ) and ( p. 15 ) are respected, physical consequences for observables near 
the W- or Z-poles remain unchanged. When constraining the electroweak parameters, 
however, we often refer to the weak currents at zero momentum transfer. The masses 



in the LEP convention are more appropriate to use in this case ||38|| , since they absorb 
reducible higher order contributions from the W and Z widths. 



2.2 Vertex and box corrections 

In this subsection, the vertex and box corrections are numerically estimated in the 
SM, while their explicit forms are given in appendix A. First the neutral current (NC) 
amplitudes near the Z-pole and at low energies (|g 2 | <C m 2 z ) are discussed, then the 
charged current (CC) amplitudes in the zero momentum transfer limit. Except for the 
Zb L b L vertex, all the vertex and box corrections are assumed to be dominated by these 
SM contributions in the following analysis. 

Four types of vertex form factors appear in the NC amplitudes (|3g). T{ and T J 2 
appear both in the 7// and Zff vertices, while and T{ appear only in the Zff 
vertices: 



T 7/ V) = -i{Q,[l + r{(g 2 ) ] + h f T f 2 (q 2 )} 



(2.16a) 



9z {(/a/ " Q f s 2 ) [ 1 + T{(g 2 ) ] + h f [c 2 Y f 2 {q 2 ) + r{(g 2 ) ] + r{(g 2 )} . (2.16b) 



The SM contribution to the vertex form factors that are non-vanishing at one-loop order 
are r{ a (g 2 ), T^iq 2 ) and T{ L (q 2 ). They can be expressed by 



T{*(q 2 ) 
Tt(q 2 ) 

rf(g 2 ; 



(9r 



V 47T 



1 IZ 



(s 2 



( 9l 



Zff, 2 



V 47T 

/' 



r{ z (q 2 



E 



9l 



Wff 



Atx 



1 IW 



(<1 2 ) 



9l 



Wff 



47T 



E 



9l 



Wff 



An 



with the gauge-boson-fermion coupling convention 

•jf f = eQ fl g z R " = g z {-QfS 2 ) 

9z (hf - Qf s 2 ) , 



iff 
9l 



zff 
9l 



Wff 

9l 



9 

V2 



v ff 



(2.17a) 
(2.17b) 

(2.17c) 
(2.17d) 



:2.18) 
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Explicit forms of the functions T\ v (q 2 ), T 2W (q 2 ) and Tj nW (q 2 ) are given in eqs. ( |A. 18| ) 
( |A.20| ) in appendix A. 2. External fermion self-energy corrections are included in the 
functions r{ z (g 2 ) and T{ w {q 2 ). For right-handed fermions T^g 2 ) = holds, since only 
those diagrams with W exchange contribute to the vertex function T2 at one-loop order. 
The vertex functions T^s) are found to be proportional to the square of the fermion mass 
inside the loop, and are non-vanishing only for f a = hi in the SM, within our approxima- 
tion of using diagonal KM matrix elements and neglecting terms of order (rrib/m z ) 2 a. For 
large m t [m 2 3> |s|), the SM contribution to T^s) is proportional to m^/m^ [p9| , |40[1 . The 
functions T{ a (s) can, in general, be present, but happen to vanish for all f a in the SM; they 
are, however, found to be non-vanishing in some extended models such as the minimal 
SUSY-SM. These analytic expressions agree with the known results of refs. 



The numerical values of the vertex form factors T{(q 2 ), T 2 (q 2 ) and T^{q 2 ) at q 



m 



are given in Table 0. All the numerical results presented in this section and in the following 
sections are obtained by setting 



4vr/e 2 = 128.72 
s 2 = 0.2312 



(2.19a) 
(2.19b) 



with e 



gs 



g z sc in the one-loop correction terms. They are fixed by using the 
renormalization conditions ( |2.13D and the SM predictions for a(m z ) and s 2 (m 2 z ) at m t = 
150 GeV, m H = 100 GeV, a s (m z ) = 0.120 and S a = 0. We emphasize that we do not 
change the numerical values of eq. (|2.19|) when discussing experimental constraints on the 
charge form factors e 2 (m|) and s 2 (m|). All our predictions for the Z parameters can be 
reproduced simply by using the numerical values listed in Table || and eq. (|2.19|) , together 
with the imaginary parts of the gauge boson propagator corrections 



« S ( m z)Ms 


0.11 0.12 0.13 


A 77 (m|) 
Azz(m 2 z ) 


0.01726 0.01760 0.01763 0.01766 
0.00248 0.00257 0.00257 0.00258 
0.00005 0.00003 0.00003 0.00003 



(2.20) 



which are obtained by using the perturbative order aa s approximations of appendix A 
with the effective quark masses of eqs. (p.25| ) and ( |B.26| ). It is worth noting that the 
2 We note the following misprints in ref. |40| . In the last line of eq. (2.7), the factor 
l/(m 2 - M 2 ) 2 should read l/(m 2 - M 2 ). In the first line of eq. (2.8), the term 4g 4 M 2 
should read 4g 2 M 4 , and in the last line of eq. (2.9), the term m/2M 2 should read m 2 /2M 2 . 
Our vertex functions T J 1W , T 2W and T J mW are then related to their functions p, A and H 

r! ft 

by the identities : T[ w (q 2 ) = p (-q 2 , m 2 w , m 2 f ,), T 2W (q 2 ) = [p + A] (-q 2 ,m 2 v ,m 2 ,) + 



2 [B (q 2 ;m w ,m w ) - B o (0; m w , m w )], and T f ^ w (q 



2W \i . 

2, __ _ 2 
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real part of the vertex corrections (Table [|) and the imaginary part A 77 (m|) interfere 
with the leading Z-pole amplitude: the latter contribution has been subtracted in the Z 
parameters [26|, whereas the former contributions modify the scattering amplitudes by as 
much as 0.5%, and hence they can contribute to the cross sections at the 1% level. 

Note further that the vertex correction without the pinch term subtraction || 41 

(q 2 ) function by 

(g 2 ) = It (g 2 ) - f-R e [B (q 2 ; W, W) - B (0; W, W)) (2.21) 



T f 2 L (q 2 ) is related to the v{ L( 1 



in the t' Hooft-Feynman gauge. The difference is universal (/-independent) and we find 
'm|) — 0.00134. The vertex corrections are slightly larger in magnitude 



rf < 



4) 



after subtraction of the pinch term. 

It is convenient to introduce the following special form factor 



5b(s) = T\ L (s) + c 2 T u 2 L (s) + T b 3 L (s) + higher order terms 



which is treated also as a free parameter in our fit at s 



(2.22) 



m 2 z to deal with the strong m t 
dependence of the Zb L b L vertex (see also ref. P2"| , |T4] |). In this way, the importance of the 
Zb L b L vertex correction can be assessed independently of the specific SM mechanism 
and also the data analysis is kept separate from the evaluation of 5b in a specific model. In 



the SM, the parameter 8b can be evaluated by including 0(a s m 2 ) [f|6j and 0(mf) |47|,[4q] 
two-loop corrections of the SM, which are given explicitly in appendix C.4: see eq. ( C.54 ). 

At low energies, light fermion masses may not be neglected compared to the momentum 
transfer q 2 . In the limit of | q 2 \ /m 2 z <C 1 and m 2 /m 2 z <1 but at fixed m 2 /q 2 , the vertex 
functions reduce to 



J- i 



1Z 



r^(g 2 



(Q 2 



rn 



w 



rn 



w 



Jz(q 2 ;f) + 
Jw(q 2 ;f) 
Jw(q 2 ;f) 



rrtr 



O 



rrir 



(2.23a) 
(2.23b) 
(2.23c) 



The functions Jz(q 2 ', f), Jw(l 2 ', f) and Jw(q 2 ', f) have the same form as the fermionic con- 
tribution to the neutral gauge boson vacuum polarization functions: see eq. (|A.27| ). The 



form factor T^iq 2 ) is often called the neutrino charge radius term [^9|. The subtraction 



of the pinch term makes it gauge invariant ||34|| . 

For the NC process f a (Pi)f'f3(P2) —> faiVzlfpiVi)) as we U as f° r its crossed channels, 
the box correction terms in eq. ( |2.2|) can be expressed as 

1 



16vr 2 



n Zff Zf'f 



[h(u, s; m z , 0) - I 2 (t, s; m z , 0)] 



(2.24a) 



16 



16tt 2 



B°t a (s,t) 



16vr 2 



Wff" Wff" 

9l 9l 



9 Z a ff 9 Z - f J 



x 



+I 1 (u,s;m w ,mf>») 
-I 2 (t,s;m w ,m f ,„) 



for hfhf < 



for 



[I 2 (u, s; m z , 0) - h(t, s; m z , 0)] 



(2.24b) 
(2.24c) 



where s = (p± — P3) 2 , t = (p± — p^) 2 and u = (pi + P2) 2 are the Mandelstam variables 
satisfying s + t + u = 0. In the second term of ( |2.24b| ), /" and /'" are the weak isospin 
partners of / and /', respectively, where all external and internal fermion masses except 
for mfin are neglected : the upper term (-Z3/-/3// < 0) should be taken for (/, /') = (£,u), 
(v,£), and (is, d), whereas the lower term (I^fl^f, > 0) for (/,/') = (£, d), and (v, u). 
The explicit form of the box functions l\ and I 2 are given in eq. (|A.30|) of appendix A. 3. 
These analytic expressions agree with the known results of refs. [^0|-|44]| . It is worth noting 
here that the box contributions to the helicity amplitudes have the above simple current 
product form only when the external fermion masses can be neglected. 

The numerical values of the box functions Bij(s,t) for the process e + e _ — > ff are 



given in Table f| for s = — 2 1 



rrir 



They contribute negligibly to the Z parameters, 



because they do not interfere with the dominant Z-pole amplitudes being almost purely 
imaginary near the pole. The imaginary parts appear in the box functions only above the 
V^-pair production threshold. 

The box contributions are found to be non-negligible in some low energy NC processes. 
In the s = t = u = limit, one finds 



7i(0,0;m v ,0) = 
J 2 (0,0;m y ,0) 



m 



2 ' 



mi 



(2.25a) 
(2.25b) 



The WW box contributions to the processes with the I\ function, that is, the low energy 
u—£, v-d and e-u scattering processes are found to be significant. 

Precise values of the charged current matrix elements are needed only at low energies. 
The muon decay constant is given by 

9 2 w(0)+9% 



G) 



AV2m 2 



(2.26) 



where the factor 5, 



G 



9 



2 r 



'G 



1 + 



4 s 2 



- 1 In — 



0.0055 



(2.27) 



denotes the sum of the vertex and the box contributions in the SM. Its numerical value 
above is obtained for the couplings of eq. (|2.19|) . The identity ( |2.26|) gives the physical W 
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mass in terms of g^(0), once the S G value is known for a given model. The overline here 
again indicates the removal of the pinch terms with the consequence that the numerical 



value is significantly (about 20%) smaller than the standard factor [50 

*2 r , 7 



>G 



9 



8tt 2 



3 - 



4s s 



in I 

cr. 



0.0068 



(2.28) 



which was obtained simply by subtracting the singular vertex function at zero momentum 
transfer. The difference 

^2 

s2 



5 G - S G 



9 
4vr 2 



B (0; W, W) - c 2 B (0; W, Z) - s 2 B (0; W, 7) 



[2.29) 



is the pinch term contribution [|34}| . Note that the sum of the propagator and the ver- 
tex/box corrections is scheme-independent and that the correction term 5 G of eq. ( 2.27 ) 
should be used together with the charge form factor <7jy(0) which contains the associated 
pinch term. 



2.3 Constraints due to a, Gp and m z 

Among the electroweak observables the three quantities a, Gp and m z have been 
measured with outstanding precision, namely Aa/a ~ 5 x 10 -8 , AGp/Gp ~ 2 x 10 -5 f25 



and Am z /m z w 8 x 10 -5 ||26|| . For this reason a, Gf and m z are chosen as our basic 
electroweak parameters and treated as constants in the analysis (see Tables 1 and 2). 

On the other hand, the tree-level properties of the gauge boson propagators are fixed 
completely by three parameters, the two gauge couplings g and g' for the SU(2)l and 
U(1)y gauge groups, respectively, and one vacuum expectation value v = (y2 Gf)^ 1 ^ 2 ~ 
246 GeV, in models where the electroweak symmetry breaking sector has the custodial 
SU(2) symmetry ||5l| . Consequently, the four charge form factors e 2 (g 2 ), s 2 (q 2 ), gz{<l 2 ), 
and giviq 2 ) are completely determined by finite quantum corrections in this class of models 
when expressed in terms of the three constants a, Gf and m z . 

In this subsection, the prescription for calculating all charge form factors in terms of 
(a, Gf, tti z ) is given explicitly in an arbitrary model with the broken SU(2)l x U(1)y 
gauge symmetry. Their numerical predictions are given in the SM. 

The form factor e 2 (q 2 ) = 4na(q 2 ) is fixed by the following identity 

^- 1 = 4vrRek> 2 )-iT?^(C 



;2.30) 



which gives the renormalization group improved running a(q 2 ) as explained in appendix B. 
Here II T 1 
eq. Q. 



Here n T (q 2 ) is the 77 propagator function without the overall coupling factor e 2 ||: see 
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In principle, the effective coupling a(m 2 z ) can be calculated from the observed a value 
by using the above identity. In practice, however, the right-hand side suffers from non- 
perturbative QCD corrections to the light quark contributions. We make use of the 
dispersion analyses [27-29| to estimate the hadronic contributions to the running of a(q 2 ) 
and s 2 (q 2 ) form factors at < \q 2 \ < m z . Details can be found in appendix B. 

In order to take account of uncertainty in the hadronic contribution and also possible 
new physics contributions, the parameter 8 a is introduced as an external parameter in 
the analysis : 

1 



a(m|) 



= 128.72 + 5 



which can be expressed by 



u a ~ 



Shad + 4vrRe 



TT H (m 



n?>) 



NewPhysics 



(2.31) 



(2.32) 



for m t = 150-200 GeV as stated in eqs. (|B . 3 2| ) and (|B.30| ) of Appendix B. Here 5had = 
0±0.10 ( |B.22| ) is the present estimate |28| for the hadronic contribution. The parameter 
S a being treated as an external parameter serves also to assess future improvements in low 
energy e + e~ hadroproduction experiments as well as possible new physics contributions. 

The remaining three charge form factors can be fixed by introducing the three radiative 
parameters S, T and U that are defined by the following identities : 



= 1 — aT . 



m 2 w 0z(O) 
47r s 2 (m|)c 2 (m|) 



47T 



772, 



a(m% 



m 



s 



s + u 



(2.33a) 
(2.33b) 
(2.33c) 



9w( Q ) a { m z) 

The parameters S, T and U can be calculated perturbatively in any models from the 
gauge boson propagator functions of eq. ( |A.1| ) by 



S 

T 

U 



16 7rRe 

a 

16vrRe 



rf 3< 3 i 2 \ 



f33 



n T ,z(o) 



Re 



nJ 3 (o)-n"(o) 



4 3 z (o) - (o) 



(2.34a) 

(2.34b) 
(2.34c) 



For models without custodial SU(2) symmetry, the T parameter is sensitive to the ultra- 
violet cut-off, and hence is un-calculable from (a, Gf, rn z ) alone. In this case it should 
be regarded as the fourth basic parameter of the theory. 
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Our definitions ( p. 33 ) of the three parameters S, T, U are inspired by the pioneer- 



ing work of Peskin and Takeuchi Q. Our definition, in contrast to theirs, includes all 
radiative effects from both SM and new physics contributions. The original parameters, 
denoted below by the index PT, are approximately related to ours by subtracting the SM 
contributions evaluated at m t = 150 GeV and m H = 1000 GeV : 

(2.35a) 
(2.35b) 
(2.35c) 

provided the scale of new physics is much larger than m z . The expressions ( |2.34j ) agree 
with the modified S, T, U parameters of ref. [j4j| . The same form of the definitions 
without the pinch terms (in the 't Hooft Feynman gauge) have been used in some earlier 
works HIM H- 

Explicit forms of the SM contributions to the S, T, U parameters are given in ap- 
pendix C, together with the SM contribution to the Zb L b L form factor 5b{m 2 z ). All 
the known two-loop corrections of order aa s [[46"], [54H5"6"1| and order [47], 48], |57|, [58| are 



SpT 


« S- 


SsM(m = 


150 GeV, m H = 


1000 GeV) , 


Tp T 


sa T - 


T S M(m t = 


150 GeV, m H = 


1000 GeV) , 


[7pT 




Usuim = 


150 GeV,m^ = 


1000 GeV) , 



included. The recently found |5(|. ^9[ small two-loop corrections of order m 2 H are ne- 
glected. For practical reasons we adopt the perturbative order aa s [|46],|54|-^(| corrections 
at a s = ot s {m z )y^ in calculating all the parameters S, T, U and 5b(m z ). The reader 
can therefore unambiguously reproduce our results. The effects due to non-perturbative 
threshold corrections |6"0" H6"2]| should be evaluated separately, and one can obtain more 
precise predictions of the SM from our formulae by adjusting the effective top-quark mass 
to produce the same S, T, U, and 5b(m z ) values. It should be noted that at present the 
uncertainty in the SM contribution to the T parameter is such that m t can be predicted 
with a few GeV uncertainty for a given T value [62]. Fig. |T] shows the SM contribu- 



tions to the S, T, U and 5b{m z ) parameters as functions of m t for m H =1-1000 GeV at 
a s (m z )=0.12. It is worth noting that the T and Sb{m z ) parameters are proportional to 
m 2 for large mt (m 2 3> vr? z \ the parameters U and 4(m|) are almost independent of m H , 
the T parameter decreases with increasing m H , and the S parameter becomes negative 



for small m H . 



Once the S, T, U parameters are calculated in a given model, the three charge form 
factors can be predicted as follows : 



1 1 + 8 G - a T 



g 2 z (0) 4V2G F m 2 z 



(2.36a) 
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g { m z) 
e 2 (m|) 



16 7T 



(5 + 17) . 



(2.36c) 



The expression ( |2.36a| ) follows from eqs. (|2.33a|) and ( 2.26 ) up to terms of order a 2 . Its 



explicit form takes account of the reducible order m\ corrections |p3|l , and it makes clear 
that the combination 



5 G — aT 



(2.37) 



determines the neutral current charge form factor g z {q 2 ) in terms of Gpw? z . In fact, 
the pinch term contribution to T in eq. ( |2.34b[ ) and the one removed from the vertex 
contribution in 5 G ( |2.27| ) cancel in the combination 5 G — aT. 

It is clear from eqs. ( |2.36|) that <?f (0) is fixed by 5 G — aT, s 2 (m 2 z ) by gf (0), «(m|) 
and 5*, and #jy(0) by s 2 {m 2 z ), a{m 2 z ) and S + U. It is instructive to express these form 
factors approximately as linear combinations of the parameters S, T, U and S a : 



^|(0) = 0.5456 + 0.0040T, 

s 2 (m 2 z ) = 0.2334 + 0.00365 - 0.0024T 
g^(0) = °- 4183 - 0.00305 + 0.0044T + 0.0035f/ + 0.00145 o 



0.00265 n 



(2.38a) 
(2.38b) 
(2.38c) 



Expressed in this form, it becomes obvious that essentially (?|(0) measures T, s 2 (m 2 z ) 
measures S — 0.7T, and <7vr(0) measures T + 0.8U — 0.7S, if the SM values of S G and 
8 a are assumed. Here the coefficients are obtained by setting 5 G = 0.0055. Results for 
arbitrary 5 G are obtained by the replacement : 



T —> T + 



0.0055 - 5, 



G 



a 



(2.39) 



Note that the combination 5 G — aT vanishes in the SM (5 G ~ 0.0055) for T w 0.75. Fig. ^] 
shows that this cancellation occurs at around m t pa 175 GeV. The SM predictions for 
the neutral current experiments can then be reproduced rather accurately by using the 
'tree-level' predictions with 5 G — aT = S = in eqs. ( |2.36|) , since the SM contribution to 
S is rather small. This should not, however, be interpreted as absence of any quantum 
corrections |M] (that is, 5 G = T = 0), but rather as evidence for the large quantum 
correction aT pa 0.0055 within the SM (see also section 5.3). 

Finally, the running of the remaining three charge form factors are calculated by 



s 2 {q 2 



[m. 



m 2 ) si(o) 



Re 



Re 



n r, 7 W ) - n T, 



(2.40a) 



U? <z (q 2 ) - 4>) 



2 s 2 Re 



Tl 6 T Q z (q 2 ) - l4>) 
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+s 4 Re 
Re 



TKV)-nS(o) 



n 11 <~ 2 

L1 TW 



(q 2 ) - <w (o) 



^2.40b) 
(2.40c) 



1 1 

Equation (|2.40a|) is the solution of the RG equation (see appendix B), and hence is valid 
at arbitrary q 2 . At \q 2 \ < m%, the parametrizations of the dispersive analysis [[27-29 



arc 



used for the light quark contribution. Equations ( j2.40b ) and ( |2.40c ) are valid perturbative 
expressions provided \q 2 \<zO(m 2 z ). At very high energies (|g 2 | 3> Tn z ), the more elaborate 
expressions ( p.38|) - ([B.41J) should be used to estimate accurately the charge form factors 
g 2 z (q 2 ) and g 2 w {q 2 ). 

Fig. H displays the four charge form factors l/d(q 2 ), s 2 (q 2 ), g z (q 2 ) and glv^q 2 ) as 
functions of 



for both time-like (q 2 > 0) and space-like (q 2 < 0) momenta. They 
are obtained in the SM for several m t and m H values, namely m t = 100, 150, 200 GeV, 
and m H = 100, 1000 GeV. The trajectories are fixed such that the known values of 
the three basic parameters (a, Gp, and m z ) are reproduced for 5 G = 0.0055, 5 a = 
and a s {m z ) = 0.12. The running of the form factors a(q 2 ) and s 2 (q 2 ) at \q 2 \ m 2 z 
is due to the QEDxQCD quantum effects [P5[ , and its detailed treatment is given in 
appendix B. The threshold effects are clearly seen in the time-like trajectories. Light 
hadron threshold effects do not show up since we adopt the dispersion integral fit of the 
hadronic contributions to the vacuum polarizations in the space-like region fl27|-p9f1 also 
for their contribution in the time-like region. The running of the g z (q 2 ) and g\y{q 2 ) form 
factors freezes at \q 2 \ m 2 z . It is clearly seen that the weak boson threshold effects are 
significant for all the charge form factors in the time-like region^. 
In section 4, the charge form factors (|2T 



and (|2.7|) are determined from the three 
sectors of the electroweak precision experiments under the assumption that there are no 
new physics contributions to the vertex and box corrections, except for allowing the Zb L b L 
vertex to take on arbitrary values. 



3 Predictions of electroweak observables 

In this section, all electroweak observables are expressed in terms of the helicity ampli- 
tudes of eqs. ( |2.2| ) and ( |2.6| ), together with the external QED and QCD correction factors. 
3 Note that the charge form factor g^iq 2 ) suffers from an infrared singularity at q 2 = 



m 



due to the opening of the W + 7 threshold on the pole ||65|| . The charged current 



cross section near the iy-pole may be expressed more conveniently in terms of gj^/(0), or 
Gpmly. 
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The predictions are restricted to the models respecting SU(2) L x U(l)y gauge symmetry 
with spontaneous breakdown to U(1)em an d presented as functions of the charge form 
factors s 2 (0), <?f (0), ^(0), s 2 (m|), g|(m|), and the vertex form factor 5 b (m|). It is 
assumed that the remaining vertex and box correction are dominated by the SM contri- 
butions. 



3.1 Z boson parameters 

The following observables on the Z-pole (s = m z ) are used in the fit 

Tz, cr°, Re, Apg, P T , v4 L r, A F ' B) A F ^, R b . 



(3.1) 



Since the Z mass m z is measured very accurately, the value m z = 91.187GeV is treated 
as a constant in the fits. The contributions from the SM box corrections are very small 
on the Z-pole (see Table ||), thus the cos ^-dependence of the box correction factors is 
neglected. 

The total cross section for the process e + e~ — > ff is given by 



a/ = a{e e 
s ( 
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♦//) 

M e L [ + M e L s R 2 + 



2\ C 



fv 



+ 



Mil ~ M LR 
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Kl - M RR 



r*f 
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2\ C fA 



at(s) 
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(3.2) 



for unpolarized beams, where the last term proportional to a(s)/ir accounts for the final 
state QED correction. Here and in the following 



<%(S = (Pe- +Pe + )\t= ( Pe - ~ Pf Y) (3.3) 

denote the NC amplitudes of eq. (|2.2j ). The factors C q y and C q A for quarks contain the fi- 



nal state QCD corrections for the vector [66] and axial vector current |)7],|68[ contributions, 



respectively, together with the finite mass corrections of the final state fermions pG 



a 



qV 



2 

i 12 " 



+ a + 1.409 a 2 - 12.767 a s 



K (3 - /%) 

2 

m 2 AJ~S) /no 



X 



a + 8.736 a* + 45.146 a A 



C qA = 3< /?; + a + 1.409a^- 12.767a d -6 



=F<r 



/(mt)+6 !2M)(3 + .n4)l) 
s v m z 7 J \) 



s 

2 -i 



-a + 14.286 a' 



(3.4a) 



(3.4b) 
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with 



= „.( 5 ) 



MS 



7T 



A, 



4^2(^5) 



(3.5) 
(3.6) 



where rh q (y/s) denotes the MS running quark mass at \i — y/s. The masses of the 
three lightest quarks (u, d, s) are neglected, while the bottom and charm quark running 
masses, rhb{\fs) and m c (y / i), are obtained from rh b {m b ) and m c (m c ) by the two-loop 
renormalization group equations: 



rh b (y/s) _ rn c {^/s) 
m b (m b ) m c (m b ) 

rh c (m b ) 



a^(x/i) 



m b ) 



m c (m c ) 



a^(m b ) 



7 (5) 

~J5T 



"o 

J 4 ) 



■fo^ + fofV^y/i) 
.&S 5) +6f ) a( 5 )(m 6 ) 

& ( 4 ) +& ( 4 ) a (4) K) - 



»(5) „(5) 

jj Zq 

~J5) ,(5) 



1 



"0 



(4) (4) 
Ti 7n 



4 4) 



J 4 ) 



(3.7a) 
(3.7b) 



where 70, 71 



( n/ )_33-2n/ ( n/ ) _ 153 - 19n/ 



h y f 



7^ 



6 



Ti 



for 

( n/ ) 303 — 10 n/ 



12 



36 



(3.8a) 
(3.8b) 



are the coefficients of the /5-function and the anomalous mass dimension in the effective 
n^-flavor QCD. The running coupling a^ A \jj) of the effective n/ = 4 theory is calculated 
from a given a^ 5 ^(m z ) = OL s (m z )y^/n by solving the three-loop QCD renormalization 



group equation with the two-loop matching condition [72 



a^(m b )=a^(m b ) + ^-U 5 \m b f 

36 L 



(3.9) 



at /i = m b . The relation between the MS quark mass rh q (m q ) and the physical mass m q 
is given in ref. [73] as 



m q [m q ) 



l + ^ a (n f ) {mq)+Kq[a (n f){mq)] 2] 



(3.10) 



with K b 12 A (rif = 5) and K c w 13.3 (n/ = 4), for bottom and charm quarks. The 
following table summarizes the running quark masses, m b (fi) and m c (fi), for Oi s {m z )y^ = 



0.11, 0.12, 0.13, m b 



4.7 ± 0.2 GeV, m c 



1.4 ± 0.2 GeV and m z = 91.187 GeV (the 
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difference mj — m c is fixed to 3.3 GeV [74] in evaluating m c (/x)) 



"4 m z)Ms 


0.11 0.12 0.13 


m c (GeV) 


1.40 ±0.20 1.40 ±0.20 1.40 ±0.20 


rh c (m c ) 
rfi c (m b ) 
rh c (m z ) 


1.13 ±0.18 1.03 ±0.19 0.86 ±0.20 
0.90 ±0.17 0.76 ±0.17 0.56 ±0.17 
0.65 ±0.13 0.53 ±0.12 0.37 ±0.12 


m b (GeV) 


4.70 ±0.20 4.70 ±0.20 4.70 ±0.20 


m b (m b ) 
m b (m z ) 


4.17 ±0.19 4.06 ±0.18 3.92 ±0.18 
3.05 ±0.16 2.83 ±0.15 2.59 ±0.15 



(3.11) 



The function f(m t ) in the 0(a 2 ) axial part of eq. ( |3.4b| ) is given by |56 



f(m t ) = 2 In 



37 
12 



28 
81 



rrir 



- 0.5767 



rrir 



+ 0.7873 



rrir 



2m t 



m t L'Z 81 \2m t J \2m t 
The minus sign should be taken in front of f(m t ) in eq. ( |3.4b|) for u, c quarks, 
plus sign for d, s, b quarks. These formulae are sufficient to calculate the factors 
C q a as functions of a s (m z ), m b and m c . For charged leptons, the correspondin 
are 

&(3 - PI) 



, (3.12) 

and the 
C q v and 
g factors 



with 



ft- 



Ami, 



(3.13a) 
(3.13b) 

(3.14) 



The effect of the charged lepton masses is negligible except for the r lepton. 

Near the Z-pole, s ~ m|, the cross sections are sensitive to the total Z width, T z , 



and hence it should be evaluated at two-loop level |]41|)^4],|f5|. The Z width is calculated 
in a similar way as the total cross section case (|3.2j ): 



rrir 
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(3.15) 
(3.16) 



by using the Z — > f a f a decay amplitudes 

T(Z -> fj- a ) 



Here e z is the normalized Z wave function, J^ a are the currents of eq. 
amplitudes M[ can be expressed by 

M[ = (h fa - Q f s 2 )\g z (m 2 z )+g z ReTt(m 2 z ) 



(3.17) 
H), and the scalar 



+g z Re h fa (c 2 T f 2 a (m 2 z ) + T^(m|)) - Q f (s 2 (m z ) - s 



(3.18) 
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It is straightforward to evaluate the partial and total widths from the above formulae, 
once the three form factors g z (Tn z ), * 2 (m|), 5b{m 2 z ) and a s (m z ) are given. Fig. |3] shows 
the predicted T z (GeV) in the plane of s 2 (m|) and g z {m z ) for a s {m z ) = 0.11, 0.12, 0.13 
and 6 b (m 2 z )=0 (a), -0.01 (b) and -0.02 (c). In the SM, 5 b (m z )& - 0.003 holds for all m t 
(see Fig. |), 5 b (m 2 z ) = -0.01 (-0.02) for m t « 175 (270) GeV. Q It is clearly seen from 
the figure that T z increases with growing a s and 5 b , and that it remains roughly constant 
when a s increases by 0.01 and, simultaneously, 5 b decreases by about 0.006. The net effect 
is a strong anti-correlation between the fitted a s and 5 b values (see section 4.1). 

In the SM, all the form factors are calculable in terms of m t and m H . In table |5| the SM 
predictions are shown for the partial and the total Z widths for several m t and m H values, 
for a s {m z ) = 0.12, 5 a = 0, 5 G = 0.0055 and (m b , m c ) = (4.7, 1.4) GeV. The numerical 
values turn out to be larger by about 1/5000 than those quoted in ref. Uncertainties 
in our predictions are estimated as follows : (i) Change of m b and m c by 0.2 GeV affects T b 
by less than 0.2 MeV (£1/2000 of r b ) and T c by about 0.03 MeV; (ii) Setting g z = g z {m 2 z ) 
and s 2 = s 2 (m 2 z ) in the amplitudes (|3.18| ) affects the total width by about 0.2 MeV for 
the m t and m H values of Table § (iii) If the imaginary parts in the amplitudes ( [3. 18 



arc 



also included, the total width increases by about 0.01 MeV; (iv) QCD higher order effects 
may affect the hadronic widths at the level of a A s ~ 1/4000; (v) The present uncertainty 
in l/a(m|), 8 a = ±0.1 affects s 2 (m|) by ^0.00026 ( pT38b|) , and hence the Z width by 
about ±0.65 MeV (~ 1/3000 of T z ). These uncertainties are still an order of magnitude 
smaller than the actual experimental error of A(r^) = 7 MeV |26| (A(T Z )/T Z ~ 0.003). 
Note that we adopt the perturbative order aa s [^6| , [54] -|56]| corrections at a s = a s (m z )j^ 



in calculating all the SM predictions, since it allows the reader to reproduce our results 
straightforwardly. The effects of non-perturbative threshold corrections |60-p2| may be 



4 The m r dependences of the electroweak Z boson observables are not completely ab- 
sorbed into the three form factors, g z (m 2 z ), s 2 (m 2 z ) and S b (m 2 z ). Mild m t -dependences 
remain in the two-loop QCD correction factor f(m t ) of eq. ( |3.12 ) and in the Zb L b L 
vertex function T\ L {m 2 z ) (Table |). When 8 b ( m 2 z ) is allowed to vary in the fit, these 
residual m r dependent terms are determined by using the SM m r dependence of the 
S b (m 2 z ) form factor (see Fig. |I]), which can be inverted approximately as m 4 (GeV) = 



21.77^-10 4 5 6 (m|) -9.9 - 31.2 valid in the region 75 GeV< m t <400 GeV. We set 
m t =75 GeV for 6 b (m 2 z ) > -0.0036 and m t =400 GeV for 8 b {m%) < -0.0405. With this 
prescription the parameter 5 b (m 2 z ) covers the full m r dependences of the vertex corrections 
within the SM, while it still allows S b (m z ) to measure large new physics contributions to 
the Zb L b L vertex because of the relatively mild mrdependences of the f{m t ) and T\ L (m 2 z ) 
factors. 
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accounted for by adjusting the effective top-quark mass to produce the same T parameter 
value. 

Once the Z width, T z , is determined the formula ( |3.2|) gives the total cross section 
for the process e + e~ — > // at all energies, up to the cos ^-dependence of the box form 
factors which can be safely neglected near the Z-pole. At LEP, the on-pole cross sec- 
tions a® are obtained after subtracting the 7-exchange contribution to the amplitudes. 
Because of this subtraction, we cannot simply compare 07 (m|) of eq. ( |3.2|) with the cor- 
responding published measurement. In fact, the subtraction procedure is not completely 



model- independent and the following two cases are examined : (i) In the amplitudes (2.2) 
only those terms multiplying the Z propagator factor are retained; (ii) From the full 
amplitude ( j2.2| ) the 7-exchange amplitude QiQj[e 2 (m 2 z ) — ie 2 A T7 (m|)]/s is subtracted. 
The above two prescriptions differ by contributions from the 7 vertex corrections and 
the box corrections, but the numerical predictions for a® are found to differ by at most 
0.0003 nb and are thus negligibly small compared to the actual experimental error of 
A(cr°) =0.14 nb. The pole amplitudes (i), the term with the Z propagator factor in 
eq. ( |2.2|) , are used below when confronting the theoretical predictions with the LEP/SLC 
experiments. 

It must be pointed out here that the quantities quoted as cr° by the LEP electroweak 
working group |26| are not the peak cross sections as obtained above, but that they are 
rather defined by the following identities^ 

12n T e T f 

z 



oJ(LEP) = — (3.19) 
1 mi 1 1 



This quantity does not agree with the pole cross section <r° as calculated above, but agrees 
rather accurately with the modified expression: 



o$(LEP)»o$-(l + 53£>). (3.20) 



For example, the SM predictions for m t =175 GeV, m H =100 GeV, a s (m z ) = 0.12 and 
5 We thank T. Mori for pointing out our misunderstanding of cr° affecting the earlier 
version of the present work. The notation of the LEP electroweak working group is 
misleading, since ref. |2(| does not explicitly state that their <r° value is not the peak 
cross section. In order to avoid any ambiguity it would be better to call this quantity 
{12it /m 2 z )T e T} l /Y 2 z and explain precisely from which experimental quantities it is cal- 
culated. It is also desirable to publish the total hadronic cross sections at yfs = m z 
without subtracting the 7-exchange contributions, since the full total cross sections can 
be calculated unambiguously. 
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5 a = give : 



/ 


°°t 




a° f (LEP) 


I — e, [A 


1.995 nb 


1.998 nb 


1.997 nb 


h 


41.399 nb 


41.476 nb 


41.463 nb 


b 


8.928 nb 


8.945 nb 


8.942 nb 



(3.21) 



The right-hand side of ( |3.20| ) reproduces the LEP definition ( |3.19| ) with an accuracy 
of 1/3000, while the peak cross sections a® as obtained from eq. ( |3.2| ) with the Z-pole 
part of the amplitudes are off by about 1/1000 to 1/600. The former uncertainty of 
about 1/3000 is typically on the order of the higher order corrections, while the latter 
difference, especially the difference between <r° and a°(LEP), shows up clearly in the fit 
CIS db SI gnificant shift in the fitted 5 b (m 2 z ) and a s values. 
Figs. |, |, | show 

a°(LEP) = £ a°(LEP), (3.22) 

f=u,d,s,c,b 

R £ = a°(LEP) / a°(LEP) = T h /T e , (3.23) 
R b = cr°(LEP) / cr°(LEP) = T b /T h , (3.24) 

respectively, in the plane of s 2 (m|) and 5 b (m z ) for a s (m z ) = 0.11, 0.12 and 0.13. All 
the three quantities turn out to be almost completely independent of gf (m|), as the 
predictions at g z {m z ) = 0.55 (solid lines) and g|(m|) = 0.57 (dashed lines) are almost 
degenerate. Fig. |] shows that <r° is sensitive to both ot s and 5 b , but an increase of a s (m z ) 
by 0.01 can be compensated by a simultaneous decrease of 5 b by about 0.006, just as for 
Y z . Fig. H] shows that the ratio Ri is only sensitive to a linear combination of s 2 (m 2 z ) 
and 5 b {m 2 z ). At fixed s 2 (m|), the correlated change of ot s and S b leaving T z and <j° 
unchanged, keeps also the Re value roughly unchanged. The reason for this behaviour is 
in the fact that the cv s -dependences of the three Z-resonance observables, Tz, cr° and Ri, 
are solely contained in just the quantity T h which depends on a s and 5 b approximately in 
the combination 

S b (m 2 z ) + 0.6a s (m z ). (3.25) 

Hence, in order to get a s independently of S b , the measurement sensitive to another 
combination is required. For instance, Fig. |6| shows that the ratio R b does measure 
S b (m 2 z ) rather independently of a s and s 2 (m 2 z ). An accurate measurement of R b offers 
the key to disentangle a s and S b (see also section 4.1). 

The asymmetries on the Z-pole provide the measurement of the universal parameter 
s 2 (m|) almost independently of g z (m 2 z ) and S b (m 2 z ) and with little or no dependence on 
the QCD coupling a s . 
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The forward-backward (FB) asymmetry is given by 
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(3.26) 



for leptons, and 
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(3.27) 



for quarks (q = b,c). Here, the physical heavy quark masses m q are used in the factor 
(3 q = Jl — Amy/nrg. The QCD corrections for the FB asymmetries |TB[ have not been 
included in eq. ( gg7p . The reported asymmetries from LEP A^(LEP) and A^(LEP) 
have been corrected for these effects assuming a linear a s -dependence and a s = 0.12. 
Therefore, we estimate the LEP asymmetries for a given value of a s (m z ) by using the 
following simple formula : 



A°'1(LEP)=A' 



0,9 
FB" 



l + k A 



7T 



0.12 



(3.28) 



7T 



with k& = 0.75 P6 |. The uncertainty in the coefficient Ak A = ±0.25 affects the above a s 
dependence by less than 1/1000 in the range 0.11 < a s (m z ) < 0.13. The QCD correction 
depends on details of the final charm and bottom quark tagging procedure, and hence 
it is desirable to have the a s -dependence of the corrected asymmetry value from each 
experiment. 

The r polarization asymmetry is defined by the ratio of the left- and right-handed r 
pair cross sections : 
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(3.29) 
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By neglecting the r mass one finds 
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(3.30) 



Likewise, the left-right beam polarization asymmetry is defined by 



.4 



LR 



E 



(3.31) 
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where the cross sections for completely polarized beam are expressed in terms of the 
helicity amplitudes by 



a f = <j(e r .e 



L C R 
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The cross section for the electron beam polarization P e is then 

Of(Pe. 



1-P P r 1+P. 
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We comment here that the factorization identities 
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(Pr) 



do not hold exactly even in our Z-pole approximation to the amplitudes 
do not factorize into Z production and Z decay amplitudes at s — m 2 r 



instance for the SM predictions at m t 
6 n = 0: 



175 GeV, m 



H 



100 GeV, a s (mr 



(3.32a) 



(3.32b) 



(3.33) 

(3.34) 
(3.35) 

), since they 
We find for 
= 0.12 and 



s 2 (m|) 


0.23040 


-4lr 

~Pr 


0.14801 
0.14802 


^FB 

Hp,) 2 


0.01667 
0.01643 



(3.36) 



The identity ( [3.34j ) holds rather accurately, but the identity (|3.35|) is violated by a factor 
of 1.4%. This is mainly because of the subtle cancellation among the squared amplitudes 
of eq. ( |3.26|) rendering the asymmetry Ap B sensitive to our detailed treatment of the order 
a 2 effects such as the treatment of the imaginary part and the choice of the couplings g\ 
and s 2 . 

In Fig. 0, all asymmetry parameters on the Z-pole are plotted as functions of s 2 (m|). 
For each asymmetry, the contributions from both the 7-pole and Z-pole terms are ex- 
amined using the following helicity amplitudes : (i) The full helicity amplitudes ( |2.2| ) 
including the 7 and Z exchange as well as the box contributions, (ii) The helicity am- 
plitudes obtained from the full amplitudes (|2.2| ) by subtracting the real and imaginary 



parts of the 7 exchange contribution QiQj[e 2 (m 2 z ) — ie 2 A 77 (m|)]/s. (iii) The helicity 
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amplitude retaining only the Z-pole term, the term multiplying the Z propagator factor 
in eq. Q2.2Q . (iv) The helicity amplitude in the improved Born approximation (IBA) of 
the Z-exchange amplitudes : 



9 2 z {m 2 z ) 



/ e A = »ZK ZJ L 3e a Ve V Z)\ [3fp V/ V ZJj 
V "PJlBA S - m 2 +is Iz_ 



on the Z-pole s = m 2 z . In Fig. |7|, the predictions of (i) are denoted by 'Full', (ii) by 
'Full — 7', (hi) by l Z only', and (iv) by 'IBA'. The prescriptions (ii) and (hi) give almost 
identical predictions, and we adopt (hi) in the fit. It is worth noting that the subtraction 
of the 7-exchange amplitudes affects the asymmetry Ap B significantly, but not the other 
asymmetries. Note particularly that the IBA gives consistently larger asymmetries by 
as much as 10% for Ap B , and by about 5% for the rest. Hence, the 'process-dependent' 
effective sin 2 9w factor determined from each asymmetry by making use of the IBA-like 
formula (|3.37|) differs significantly from the process-independent universal form factor 



s 2 (m|). We find approximately, 

sin 2 9$(A°/ B ) « s 2 (m|) + 0.0009 , (3.38a) 

sin 2 6$(A° LR ) « s 2 (m|) + 0.0010 , (3.38b) 

sin 2 &w( A ¥b) ~ s 2 ( m z) + 0-0010 , (3.38c) 

sin 2 C(^fb) ~ s 2 {m\) + 0.0009 . (3.38d) 

A related study is found in ref . |77| . 

In the SM, all the form factors gz{ m z)i ^ 2 ( m l) an d ^>( m i) are calculable as functions 
of rrit and m H (see appendix C for details). The main uncertainty in these calculations 
appears in the parameter 5 a (|2.31|) which parametrizes the uncertainty in the hadronic 
vacuum polarization contribution to l/a(m|). Hence, all Z parameters can be predicted 
accurately in the SM as functions of four parameters: mt, m Hi a s an d S a . 

Figs. |H shows the mt-dependence of all Z parameters for three m H values 60 GeV 
(dashed lines), 300 GeV (solid lines) and 1000 GeV (dash-dotted lines), at a s {m z ) = 0.11, 
0.12, 0.13 and 5 a = (l/a(m|) = 128.72). Shown by horizontal lines are the experimental 
data from LEP ]2E| and SLC [[H]] (see sections 4 and 6). The m t - dependence is sizeable for 
all the observables. In Rb and cr°, the m r dependence comes mainly from the Zb L b L form 
factor 5b(m z ), and hence these parameters have little sensitivity to m H (see Fig. H). The 
m^-dependences of all asymmetry parameters including P T come from the form factor 
s 2 (m|). Re receives mrdependences from both 5b(m 2 z ) and s 2 (m 2 z ). Finally, the total 
Z width is the only quantity sensitive to the form factor g z (m 2 z ). In conclusion, the 
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m t -dependence of V z is a combined effect of all three form factors g|(m|), s 2 (m 2 z ) and 

Likewise, Fig. || shows the ot s dependences of the hadronic Z parameters for the three 
m t values 100 GeV (dashed lines), 150 GeV (solid lines) and 200 GeV (dash-dotted lines), 
all at m H = 100 GeV and 5 a = 0. It can be seen that the ratio Rf, and the asymmetries Apg 
and A F b are almost independent of a s . Tz and Re grow linearly with a s because of the final 
state QCD correction factor ( |3.4|) . cr° decreases with increasing a s , since it is proportional 
to the factor Y h /Y 2 z . The ratio Ri exhibits the strongest dependence to a s . As emphasized 
above, however, the a s -dependences of all Z observables are approximately proportional 
to a common factor Sb(m z ) + 0.6a s (m z ), and hence either an accurate determination of 
Sb{m 2 z ) (via Rb) or else the assumption of SM dominance to the form factor <5 fe (m|) is 
crucial for the extraction of a s (m z ) from these experiments. 



3.2 Low energy neutral current experiments 



The data of four types of low energy neutral current experiments are analysed : 
neutrino- nucleon scattering (z/^-g), neutrino-electron scattering (z/^-e), atomic parity vi- 
olation (APV), and polarized electron-deuteron scattering (e-D). Theoretical predictions 
are given for all model-independent parameters |24],[7|,|79| characterizing the electroweak 
low energy neutral current experiments. They are the effective v^-q coupling factors [fTSfl 



2 2 ir2 r2 

9l> 9r, 5 l> o r , 

for the v^-q scattering experiments, the effective neutral current parameters |79] 



for the v^-e scattering experiments, the weak charge of nuclei [^0] 

Q W (A,Z), 



for parity violation in atoms, and the effective neutral current couplings p3 



2C*i u — Cid, 2C*2 U — C; 



2d 



(3.39) 



(3.40) 



(3.41) 



(3.42) 



for the e-D polarization asymmetry. Definitions of these model-independent parameters 
are given below and re-expressed in terms of the helicity amplitudes of eq. ( |2.2| ). 

In this subsection terms of order a ■ (q 2 /m 2 v ) are neglected, while keeping terms of 
order q 2 /m 2 v and a ■ (m 2 /q 2 ). The generic amplitude for the process ij — > ij follows then 
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from eq. 

M S C = i>{(QiQi) [*V) + e 2 ri(g 2 ) + e 2 T{(q 2 )} 

+ (Qrh J )e 2 T J 2 ( q 2 ) + (Q J I 3l )e 2 r 2 (q 2 )} 

+-r i -r{(^ - Q l s 2 ) (/ 3j - Q,^i(0) 

-{hj-Qi?)Qi& \s 2 {q 2 )-s 2 ]} 

+B» c (0,0) + O(e 2 X r ). (3.43) 

All electroweak observables of the low energy neutral current sector are calculated by 
using the above approximation. Contributions from the neglected terms are completely 
negligible. The numerical predictions for all observables (|3.39|) -( |3.42|) depend on just 



the two universal charge form factors s 2 (0) and g§(0), since the running of the charge 
form factors l/a(q 2 ) — 1/a and s 2 (q 2 ) / a(q 2 ) — s 2 (0)/a at low energies \q 2 \ -C m 2 z are 
governed completely by known physics only and are hence accurately calculable (see 
appendix B). Although the expression ( |3.43|) with the MS coupling normalization ( |2.19|) 
is used in all numerical calculations presented below, we often quote below a slightly more 
compact expression that is obtained from eq. (|3.43| ) by dropping the terms proportional 
to [s 2 (q 2 ) —s 2 ] and replacing s 2 by s 2 (q 2 ) in the term multiplying the Z propagator factor. 
This is a valid approximation to eq. ( |3.43| ) differing only by terms of order c/z[s 2 (q 2 ) — s 2 ] 2 . 

3.2.1 Neutral currents in — q scattering 

The neutral current data from the v-q scattering experiments can be conveniently 



parametrized by the four model-independent parameters [75 



Sfa = < + € . ( 3 - 44a ) 

8l = ul-dl, (3.44b) 

for a = L or R. The effective chiral couplings q a (— u l^l^ u r^r) can be directly 
expressed in terms of the helicity amplitudes of eq. ( |3.43| ) by 

L " (q = u,d; a = L,R) , (3.45) 



with the notation = Mgg. The amplitudes (ggg) can then be written in compact 



form 

M^ q = .. . 

2^ q t 2 



M L r = \Q q e 2 ^# + I \h qa ~ Q q s 2 (t)} + « (0, 0) , (3.46) 
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and the charged current factor is approximated by 

[l + 4.c.] 1/2 



a 



G f 



1 + 



(-t)aa ' 



(3.47) 



The QED correction factor 5 C c . is accounted for (following Sirlin and Marciano J8T[] ) by, 



a 

71 



In 



2 (-t)c.c. 



+ 2 



0.017 



for H) c . c . = 20GeV 2 . 



Note that the leading logarithm approach of ref. |36| gives 

,2 



- i n m w 

7T (—t)c, 

0.013 



for (-t) c . c . = 20GeV 2 



(3.48a) 
(3.48b) 

(3.49a) 
(3.49b) 



for the above correction factor. In our numerical calculation we adopt the factor (|3.48b| ). 
The fyrj charge radius factor Y v 2 (t)/t and the box form factors B^(0, 0) in the amplitude 
( p.46j ) can easily be read off from the generic expressions in appendix A : 



r u 2 (t) 



TOyy 167T 2 



Jw(t] m^) , 



with 



J w (t; Trip) = 4F 3 (t;m M ,m^) - -mm 



2 

W 



19 



m 



IV 



t m. 



??? 



2 ' 
11' 



(3.50) 

(3.51a) 
(3.51b) 



from eq. ( |A.27b| ) and 

B^r (0,0) 
Bit '(0,0) 
5^(0,0) 



+ 



3& 



64f 2 m; 64vr 2 m|V2 3 



9 3 / 



3 & 



16 7r 2 m 2 



+ 



647r 2 m| 

3j| 
647r 2 m| 

3j| 
64 vr 2 m| v ;', 



3 / ' 

1 1 ^ 
-- + -s 4 

2 3 



(3.52a) 
(3.52b) 
(3.52c) 
(3.52d) 



from eq. (A. 35). These expressions are sufficient to evaluate the helicity amplitudes ( |3.46| ) 
as functions of s 2 (0) and gf (0), for the MS coupling normalization of eq. Q2.19Q . We set 
m w = 80.24 GeV and m z = 91.187 GeV in all numerical calculations. 
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At (-t) n . c . = 20GeV 2 , 



7a~2 
16 7T 

and the g Q 's are approximated as 



-20 GeV 2 ; m 



i ••"in 



-0.0037, 



(3.53) 



q a ^ 0.9923 p h q<x -Q q s 2 (t) 



+0.0031 (q a = u L ) 

+0.0026 \q a = u R ) 

-0.0074 (g a = d L ) 

-0.0012 (q a = d R ) 



(3.54) 



for the MS coupling normalization of eq. (|2.19|) . Here the universal p parameter is defined 
by 

„_ Sl(0) i ^ fl|(0) 



5 G - aT 0.54864 



4x/2G F m! 1 

The relation between the form factor gf (0) and the T parameter is seen in eq. (|2.36 
The running of s 2 (t) is estimated as 

a(t = -20GeV 2A 



(3.55) 



-20 GeV 



a 

^2 



[s 2 (0) - 0.0097] 



1.0295 s 2 (0) -0.0100. 



(3.56) 



The approximations (|3.54j) - (|3.56|) are found to give excellent numerical predictions for all 
q^y as functions of the two charge form factors, s 2 (0) and gf(0). 

The major effects of radiative corrections can be made transparent by parametrizing 
the model-independent coupling factors of eq. ( |3.45| ) in terms of the effective couplings 
p vq and s 2 vq of ref. 

2 

2 3^' 

s 2 " 
2 ' 3 vq > 

2 



/_1 1 



Pvq\ o 8 uq) 



+ A Ui 
+ A di 
+ A M/ 
+ A, 



(3.57a) 
(3.57b) 
(3.57c) 
(3.57d) 



The extra terms A Qa are fixed such that they do not interfere with the leading terms in 
the most accurately measured quantities, that is, g\ and g 2 R . One finds 



A, 



A,. 



if 

8tt 2 
9l 



>TC Z 



2 A, 



s 2 

+ y 


) «/3L 


s 2 




~ 3 




si 




8tt 2 


5 ' 



(3.58a) 
(3.58b) 
(3.58c) 
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with 



h j3L 



'9 3 a2 8 „ 4 
2c 2 V8~2 S + 15 S 



(3.59) 



The radiatively corrected amplitudes can then be expressed approximately in terms of the 
effective strengths { p V q of the neutral current and the effective weak mixing factor 's 2 ^' in 
the v^-q scattering process. In terms of the two universal charge form factors (0) and 
s 2 (t) they are given by 



1 + 



(-t)c 



Pvq 



[l+5c.c.] 1/2 1 + ^ 

e 2 - 

4 q {t) = s (f) + Y^ J w(t;m^) - 

The box factors in eq. (|3.60|) are obtained from eq. (|3.52|) : 

5 15 
2 



8tt 2 

X2 



2 "2 7 ' 
7rCr 



(3.60a) 
(3.60b) 



9% 



9% 



--,2 



1Qtt 2 c 2 

Z2 



-s 2 - 



-s H s 

5 9 

— s H s 

10 9 



0.0074, 
a 0.0018 



(3.61a) 
(3.61b) 



4 

5 61 

87r 2 c 2 ' a7 ~ 16^6*1 2 ~ 20' 
where m^/m-l is replaced by c 2 in order to reproduce the expressions in ref. [52]. With 
the estimates ( |3.53| ) and ( |3.56| ), we find 

0.9923 p + 0.0074, (3.62a) 

1.0295 s 2 (0) -0.0155. (3.62b) 



Pvq 
°vq ~ 



These equations are useful in understanding qualitatively the effect of the v^-q scattering 
experiments off isoscalar targets, but we find that they give slightly inaccurate approxi- 
mations to the quantities q a ( |3.45| ). 

In the following table, we compare the numerical predictions for the basic quantities q a 
and the model-independent parameters of eq. (ft.44|) by using the exact matrix elements 
( |3.45| ) and by using the approximation ( p37|) , for g 2 (0) = 0.5492 and s 2 (-20GeV 2 ) = 



0.2359 (the SM predictions for m t = 175 GeV and m H = 100 GeV) : 





eq. (3.45) 


eq. (|3.57|) 




exact 


approx. 




0.3435 


0.3343 


U R 


-0.1537 


-0.1537 


d L 


-0.4260 


-0.4336 




0.0769 


0.0769 


si 


0.2995 


0.2998 


9l 


0.0295 


0.0295 




-0.0634 


-0.0763 


K 


0.0177 


0.0177 



(3.63) 
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It is clearly seen that the formulae ( 3.571 ), although reproducing u L and d L rather poorly, 
give, as expected, an excellent approximation for the most precisely measured parameter 
g\. They give, however, a rather poor approximation for the parameter 8\ being off by 
20%, which is unsatisfactory in view of the experimental uncertainty (see section 4.2.1). 

Fig. 1C illustrates the relation between the model-independent parameters (g 2 , g^) 
and the two universal form factors (s 2 (0), <?f (0)). The present data ]78| (see section 4) 
constrain the 2-dimensional parameter space to the ellipse drawn in the same figure. The 
dashed line is the p = 1 (S G — aT = 0) curve : <?f (0) = Ay/2G F m% = 0.5486. The thinness 
of the ellipse in the (gf,, g^) plane implies a strong correlation between s 2 (0) and gf (0). 
It is worth noting that the effective charge s 2 (0) derived from v^-q scattering experiments 
at q 2 « —20 GeV 2 is larger than the process-dependent effective mixing factor s 2 q by as 



much as 0.01: see eq. ( |3.62b| ). 

3.2.2 Neutral currents in — e scattering 

The total cross section for the processes u^e — > v^e and v^e — > v^e in terms of the 
helicity amplitudes M^ e and M^ are given by 



a 



a 



< - ^ ^ dz 



An Jo 



M, 



LL 



+ V--Z) 

2 



Mi 



LL 



+ 



Mlt 



An Jo 

where the variable z is related to the momentum transfer t by 

(a -m 2 e ) 2 



z = -t/t 



max j "max 



2m e E u , 



(3.64a) 
(3.64b) 

(3.65) 



with the approximation s = {pv+p e ) 2 ~ 2m e E u . The amplitudes in eq. ( p. 64 ) are obtained 
from eq. ( |3.43|) 



LL 2 t 2 

M v L t = --e 2 S^ 
LR 2 t 



gl(Q) 

t — 171% 



+ 5lY(0,0), 



+ ^)^ + «(0,0) 



(3.66a) 
(3.66b) 



where the charge radius factor T 2 (t)/t is given by eq. ( |3.51| ) and the box form factors 
B%,B% by 



5 L T(0,0) 
fl& e (0,0) 



9' 



lQir 2 mf v 167r 2 m^ v 2 



1\2 



lQn 2 m 2 z ^2 



(3.67a) 
(3.67b) 
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see eq. ( A.35Q . It is then straightforward to express the cross sections ( [3.64 ) in terms of 
the universal charge form factors s 2 (0) and gf (0). Our results (|3.66|) and ( p. 671 ) agree 
with ref. [83]. 

As in the case of the u^-q scattering analysis it is useful to introduce the process- 
dependent effective couplings p ue and s 2 e p4| : 



M, 



LL 



2V2G 



FPv 



M, 



LR 



2V2G 



F Pi 



From eqs. (|3.66| ) and ( |3.68| ) one finds 



Pv 



P + 



16tt 2 



4 

e 2 



Qs 4 



Sue(t) =S (t) + J^ J w(t; Trip) - 



16tt 2 



17 



-s 2 + 6s 4 



(3.68a) 
(3.68b) 

(3.69a) 
(3.69b) 



by neglecting higher order terms and by setting m 2 v /m 2 z = c 2 . Here p and Jw(t;m M ) 
are given by ( |3.55| ) and fl3.51| ), respectively. The cross sections can then be expressed in 
terms of the model-independent parameters p ue and s 2 e by 



a 

a 
~E~u 



V 



2m e G 2 F 2 
~Z. Pvt 

7T 

2m e G 2 F 2 

7T 



dz 



\-sUt)] 2 + (l-z) 2 [s 2 ue (t) 



dzi (1-z) 



4e(t) + 4e(t) 



where t = -2m e E v z (|3~65|) . For E v = 25.7 GeV (CHARM-II ||), we find 
In this momentum region the running of s 2 (t) is negligible: 



At) 



1.0072 s 2 (0) -0.0018 {t = -m\ 



1.0080 s 2 (0) - 0.0020 (t = -2 m 2 ) 

Also the i/-charge radius factor J w (t;m^) has little t-dependence: 

-0.0061 (t = 0) 
~J w (t,m^ = { -0.0060 (t = -m 2 ) 
-0.0059 (i = -2m 2 ) 



(3.70a) 
(3.70b) 

(3.71) 

(3.72) 



e 2 



16 7T 



(3.73) 



Thus, the t-dependence of the effective mixing factor s 2 ue (t) ( |3.69b[ ) is negligibly small. 
From eqs. (|3~69| ), ( gjj ) and (|3~73D follows 



4e(0) 



Pve « p + 0.0121 , 
2N ~ 1.0072s 2 (0) -0.0103. 



(3.74a) 
(3.74b) 
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In the limit of negligible t-dependence of s 2 e , eq. (|3.70|) becomes 

1 





_ 2m e G\ 
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Pie 


a pe 


2m e G 2 F 






71 


Pie 



- -sD 2 + -(sV 2 



2 

1 / 1 



3 ^2 



(3.75a) 
(3.75b) 



with s 2 ve = s^e(O). This is the form entering the analysis of ref. | 79fl : they combined 



the three experiments [8|] and expressed the result in terms of the model-independent 
parameters p ue and s 2 e ( |3.40| ). In our analysis the above parametrization ( |3.75|) is used to 



reproduce the combined measured cross sections from the fit [79| in terms of p ue and s ve . 
These cross sections are then analysed in our framework by using the defining equation 



Fig. |TT| illustrates the constraint by the data similarly to Fig. [10]. The approximation 
( |3. 74j ) is found to reproduce our results accurately. The dashed line denotes the curve 
p = 1 (5 G — aT = 0). The ratio of the u^e and v^e cross sections is measured accurately, 
and hence the form factor s 2 (0) is constrained fairly independently of gf (0) from the u^—e 
scattering experiments. 

3.2.3 Neutral currents in e—q interactions 



The effective Lagrangian of the parity-violating e-q interaction 
G 



2q faT^e ' ^ql^l^q 



(3.76) 



can be rewritten in terms of left- and right-handed currents as follows : 

£pv = -% E [d, (Jr - Jt) ■ (4 + Jl) + Ciq (J e R + Ji) ■ (4 - J 



V2 „ 

The effective couplings Ci q , Ci q expressed in terms of the helicity amplitudes read : 

1 



(3.77) 



C2q 



°2 9 "T Ly 2q 



2y/2G F 

1 



2 v / 2Gi 



MH ~ M r % + M e A - M e A 



LR 



RL 



RR 



, Mi) 



(3.78a) 
(3.78b) 



Here and C^j) denote the sum of the contributions from the photonic correction to 



the axial vector Zee vertex and the Z7 box correction |85[ , which are not included in our 
helicity amplitudes (|3.43j ). They are found in refs. fl85l,|8li : 



.(7) 


9 2 Z e 2 


1? - 


167T 2 






<(7) _ 


9 2 ze 2 


2q - 


16tt 2 



{-2(/ 39 - 2 Q q s 2 ) + QhqQq (1 - 4 s 2 ) (in ^§ + -) } , (3.79a) 

{-2 hq Q 2 q (I -AS 2 ) +6 Qq(l 3q - 2 Qq § 2 ) (in ^§ + |) } . 



(3.79b) 
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M 2 « (— £) « 1.5 GeV 2 is used in the analysis of the SLAC e-D scattering experiments [87 
By inserting eq. ( [3.43D into Cf l q and defined above, one finds 



2V2G F ■ = e -{-Q q ) \2{Tl L - r? a ) + 1?! (t) + 



t — m\ 



4m 



■3q 



2Q q s\t) 



Bnl + Bl\ 



B eq 
n RR J 



2V2 g f -c% = ?- (-Q q )2(rr -n 



(*) + 



(3.80a) 
Sl(0) -(-i + 2,^))(^ 



t-m z ^ 2 



LL 



R e< ? 



Re? 



R eg 



(3.80b) 



By adopting the SM predictions for the vertex and box form factors of appendix A, 
the model-independent parameters Ci q of the low energy effective Lagrangian ( |3.76| ) are 
readily evaluated as functions of s 2 (t) and g§(0). More explicitly, one finds 



4y/2G F m 2 z C^ 



t/m 



'3g 



2Q q s\t) 



+ 



9% 



167T 2 
F 



^ 167T 2 ^ 
1 + (1 



;i-4s 2 )j z + 2(j w - 

S2 



4s 



2\2 



+ 2c 2 (g 



c 



u 



2"(g = d) 



(3.81a) 



1 - t/m| 



1 - 4s 2 (t) 



16vr 2 



2 J 3g Q 9 (2 I 3<z - 4 Q g s 2 )J z + 2Q q J w -A I 3q J w ) 



9% 



16vr 2 



5(i 

2 v 



4 s 2 



(h q -Q q s 2 Y + (Q q s 2 y +2c 2 (q 



u 



(3.81b) 



where the factors Jz = Jz(t;m e ), J w = J w (t;m Ue ) and Jw = Jw(t', m v e ) are given in 
appendix A. The sum of ( |3.81| ) and ( |3.79|) agrees with ref. [|3~5" |. 
At t — —1.5 GeV 2 , s 2 (t) is calculated from s 2 (0) as 



s 2 (t) w 1.0183 s 2 (0) - 0.0058 



and the numerical values for J z , J w and Jw are Jz ~ —6.97, Jw ~ —6.80, Jvk p 
The non-universal (vertex and box) corrections for Ci q are estimated numerically as 



(3.82) 
-7.69. 



1M1 
'lul 



IPliiJlBA 



+ 0.0061 + 0.0007, 



-0.0011 + 0.0009, 



(3.83a) 
(3.83b) 
(3.83c) 
(3.83d) 



KIiba + 0-0082 + 0.0048, 
[CmW " 0.0070 + 0.0043, 

where the second terms in the r.h.s. denote the electroweak vertex/box corrections for 

J iq 



C^f, and the last terms denote the external photonic corrections, C-j\ The improved 
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Born expressions [C'* f ] IBA , s can be expressed by 



[C 



Mi 

lgJlBA 



1 - t/m% 
P 



2gJlBA 



1 - t/m 2 z 

with p = g 2 z (0)/(AV2G F m 2 z ) as in eq. (gjp . 



h q -2Q q s\t) 



l-As 2 (t) 



(3.84a) 
(3.84b) 



In the polarized eD experiment only the combinations 1C\ U — C\& and 2C*2 U — Cid 
[24| are well measured. A model-independent determination of these two combinations 
is performed in section 4. Fig. [12] shows the relation between the model-independent 
parameters i2C\ u — Cid, 2C 2u — C 2 d) and the two universal parameters (s 2 (0), (0)), 
together with the l-a contour of the result of the analysis obtained in section 4 from the 
experimental data ||87|| . Note that the vertex and box corrections (especially the WW 
box contributions) in eqs. ( |3.80| ) are important in these combinations yielding : 



2Ci u — C 
2C*2m — C' 



id 



2d 



[2Cfl - C^]iba + 0.0134 + 0.0005 , 
[2C 2 ^ - Oba + 0.0234 + 0.0052 . 



(3.85a) 
(3.85b) 



As before, the second terms denote the vertex/box corrections in while the last terms 



denote contributions from Cj q . The majority of the non-universal contributions above 
come from the WW box diagram. Since the typical contribution of the improved Born 
approximation to these factors are [2C£[ - C^] 1BA » 0.7089 and - C^] mA « 0.0751 
for g 2 z (0) = 0.5492 and s 2 (-1.5GeV 2 ) = 0.2375 (the SM predictions for m t = 175 GeV and 
m H = 100 GeV), the non-propagator correction terms are appreciable in these observables. 

In the case of atomic parity violation the momentum transfer is so small that the 
matrix elements for nucleons should be calculated. Marciano and Sirlin |8H] introduced 
effective couplings C\ p and C ln for nucleons, which may be separated as in eq. ( |3.78| ) 



C 



f~1 flM _j_ /"ft. 



(3.86a) 
(3.86b) 



Here Cfl and are the contributions from the neutral current amplitudes Q3.43 ), which 
can be expressed in terms of and Cfj by 



U lp — Z U IU "+* U ld 5 



M 



or more explicitly, 
(aV2G f m 2 z ) C 



M 
Ip 



C* 



-2s 2 (0) 



(3.87a) 
(3.87b) 
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(4V2G F ml)C?l = gU0) 
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16tt 2 1 16 
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+ 
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9% [9 



16vr 2 1 16 



1 -4s 
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2\2 



1 _ fl + (1 - « 2 ) 2 
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2 , m 2 e 1 
In 



3 ~~ ml 9 ' 



(3.88a) 
(3.88b) 

(3.89a) 
(3.89b) 



which can be obtained from eqs. flPTD by taking the g 2 -> limit. c£ } and are the 
contributions from the photonic correction to the axial vector Zee vertex and the Z7 box 
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4V2G F ml Ctf = @£ {-(1 - 4 S 2 ) 2 + 5(1-4. 



4V2G F m|^ = f^ 

167T Z 



l_4s 2 ) 2 + 4(l-4S 2 ) AT+(ei)fl 



(3.90a) 
(3.90b) 



The last terms on the right-hand sides of eqs. ( |3.90| ) denote the 'jZ-box corrections which 
are sensitive to the nucleon structure. The constants K, (£i)# and have been 



estimated in ref. [531 to be 



K = 9.6 ± 1 , {if B = 2.55 , (Ob = 1-74 . (3.91) 
By estimating numerically the vertex/box corrections in eqs. Q3.88 ) and ( |3.90|) , we find 



Ci p 



IBA 



L 1 



IBA 



+ 0.0107 + 0.0027, 
0.0038 + 0.0023, 



(3.92a) 
(3.92b) 



where the second terms denote the weak vertex/box corrections for and C^, and the 
last terms denote the photonic corrections, C$ and . The improved Born approxi- 



mations [Cj!f] IBA and [C^] IBA are given simply by 



t^+pllBA 



Ml 

In J IBA 



;-2s 2 (0) 
1 

2 



Their typical numerical values are found to be [Cfp 



I IBA 



(3.93a) 
(3.93b) 



0.0223 and [C^U, « -0.5005, 



Ira J IBA 

for <?|(0) = 0.5492 and s 2 (0) = 0.2389 (the SM predictions for m t = 175 GeV and 
100 GeV). Note that the non- universal corrections are important especially for 



m 



H 
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Ci p , where the effect comes mainly from the WW box contribution, or the term with the 
factor |c 2 in eq. Q3.88a| ). 

The weak charge Qw{A, Z) of an atom is given in terms of C\ p an C\ n by 

Qw(A,Z) = 2(A - Z)C ln + 2ZC lp , (3.94) 

which in the case of cesium is 

QwCllCs) = 156 C ln + 110 C lp . (3.95) 

Numerically they are estimated as 

Qwr(Cs) « #f(0) [ -41.92 - 400.99s 2 (0) ] + 1.77 + 0.65 . (3.96) 

where the first term comes from the IB A approximation to CM, and ( 3.93|) , the 
second term comes from the electroweak vertex/box contributions to them, and the last 
term from the external photonic corrections of eq. ( |3.90| ). It is clear from the above result 
that the vertex and box corrections should be carefully taken account of in extracting the 
electroweak parameters from the Qw measurements. 

In Fig. [13], the parameter Qw(Cs) of eq. ( |3.95| ) is shown as a function of the two 
universal parameters (s 2 (0), g|(0)) in the range 0.20 < s 2 (0) < 0.26 and 0.52 < #1(0) < 
0.57 together with the l-a contour of the data |S0[ (dashed lines). The horizontal straight 
dashed line denotes the line p — 1 (5 G — aT = 0). It is worth noting that the correlation 
in the s 2 (0) and #|(0) form factors obtained from the cesium weak charge Qvk(133, 55) in 
the figure is opposite to that obtained from the v^-q scattering experiments (see Fig. |T~0D . 
The cesium Qw measurement implies an ant i- correlation between s 2 (0) and #|(0), or p. 
This is opposite to the trend observed for the constraints from the v^-q experiments. For 
further discussion, see section 4. 

3.3 Charged current experiments 

In the charged current sector we consider two precision experiments: the muon lifetime 



25 and the W boson mass measurements 25. BS 



From the matrix element (|2.6j ) one finds for the muon decay constant 



AV2 



in 



vv 



where the factor 5 G denotes the sum of the vertex and the box contributions. It has been 



calculated within the SM in ref. |50 



g 2 r / 1 \ 1 



^-s^r-vi^-v 111 ^' (3 - 98a) 

^ 0.0055, (3.98b) 
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where the pinch term [34] has been subtracted as explained in section 2: see eqs. ( |2.27|) 
and (|2~28| ). 



The expression (|3.97|) enables one to predict the physical W mass in terms of the 
charge form factor gjy(O). Numerically, one finds : 

■*-*g^ (3-) 
15155.9 ^(0) + 46.7 6 ° - ° M55 + 35.21 GeV 2 . (3.99b) 

Once the numerical value of 5 G factor is known, the measurement of the m w mass deter- 
mines directly the charge form factor g^(0). 

The form factor g^(0) can be calculated in terms of the S, T and U parameters in 
the SU(2) L x U(l)y models. Insertion of the expansion (|2.38c|) leads to 

m w (GeV) = 79.840 - 0.291 S + 0.417T + 0.332 U-0. 136 5 a , (3.100) 

in excellent agreement with eq. ( |3.99| ) for 5 G = 0.0055. The prediction for a different S G 
value follows from the above expression by simply making the substitution ( |2.39| ). 

Fig. 113 shows the SM predictions for m w in the plane of m t and m H , for 5 a = and 



5 G = 0.0055. In the 0(aa s ) corrections to the SM contributions to the S, T, U parameters 
a s (m z ) is set to 0.12. Changing a s (m z ) by ±0.01 affects the prediction of m w by about 
=F0.004 GeV. The mean and standard deviation of the present m w measurement (see 
section 4.3) are indicated by dashed lines. 

Note that among the electroweak observables examined in this paper, only m w is 
sensitive to the U parameter. Hence, when performing a general fit to the S, T, U param- 
eters, the mean ((U)) and standard deviation (AC/) of the U parameter are determined 
solely by the mean {(m w )) and standard deviation (Am w ) of m w : 

(U) = [ (m w {GeV)) - 79.840 + 0.291 (S) - 0.417 (T) + 0.136 S a ]/0.332 , (3.101a) 
AU w Am u/ (GeV)/0.332. (3.101b) 

Here (S) and (T) denote the best-fit values from other experiments. The present exper- 
imental error of Am w = 0.16 GeV induces AU = 0.48, while Am w = 0.05 GeV, the 
precision anticipated in future LEP200 experiments, would give AU = 0.15. The full 
error AU should be slightly larger than the above estimate, since S and T were fixed and 
set at their best values in deriving ( |3.101b| ). 
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4 Experimental data and the electroweak parameters 



Based on the formalism introduced in the previous sections the values for the form 
factors are inferred from fits to the data of electroweak precision experiments : g 2 z{ m z)i 
s 2 (m|), 5b{m 2 z ) from the LEP/SLC experiments on the Z-pole, g z {0), s 2 (0) from the low 
energy neutral current experiments at q 2 ~ 0, and g\y{0) from the W mass measurements 
at pp colliders. 

4.1 Z boson parameters 

The analysis is based on the data from the LEP and SLC experiments published up to 
the year 1993 |26|,|89],|)0|]. Discussions of the recent update from LEP |91|] and the precision 
measurement of the left-right asymmetry at SLC [[31] are postponed to section 6. 

The Z line-shape parameters resulting from a combined fit performed by the LEP 
electroweak group |26[| are : 

/ 1 -0.157 
1 



m z (GeV) 
rz(GeV) 
a°(nb) 

Ri = oVo°t 

^FB 



Pc 



91.187 ±0.007 
2.489 ±0.007 
41.56 ±0.14 
20.763 ± 0.049 
0.0158 ±0.0018 \ 

The other electroweak data used in our fit are p6[,|S£ 

P-r = -0.139 ±0.014, 
0.10 ±0.044 
= 0.099 ± 0.006, 
= 0.075 ± 0.015, 
R b = a b /a° h = 0.2203 ± 0.0027 



0.017 0.012 0.075 \ 

-0.070 0.003 0.006 

1 0.137 0.003 

1 0.008 



(4.1) 



4° 

^LR 



(SLD [90]), 



(LEP + SLD). 



(4.2a) 
(4.2b) 
(4.2c) 
(4.2d) 
(4.2e) 



Definitions of all the above observables and their theoretical expressions have been given 
in section 3.1. 

The Z mass, m z = 91.187 GeV, is treated as an input parameter neglecting its error. 
This is justified because of the smallness of the experimental uncertainty and correlations. 
For the fits to be described below a few general conditions are anticipated : (a) only three 
neutrinos (N u = 3) contribute to the invisible width of Z, (b) the perturbative QCD 
corrections with the finite quark mass effects are taken as given explicitly in section 3.1, 

(c) the vertex and box corrections are calculated in the SM and given in Table |^ and |], 

(d) the Zb L b L vertex is taken into account by the quantity 5i,(m 2 z ), which is treated in 
the fit as a free parameter just as the universal parameters gz{"m 2 z ) and s 



m 
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Various methods to determine the QCD coupling constant have led to consistent results 
with a typical uncertainty of Aa s (mz) ~ 0.01. However, this is far from making it precise 
enough to be used as a fixed input parameter, since the fitted electroweak parameters are 
found to be rather sensitive to the assumed value of a s {m z ): see, for instance, eq. (|4.3|) 
below. For this reason, and also for the convenience of GUT studies, a s = Oi s {m z )y^ is 
treated throughout our fits as an external input parameter and, consequently, the best-fit 
values of the fit parameters and the minimum x 2 are always presented as functions of a s . 
Once a precise determination of a s from independent data is available, it is straightforward 
to get the correspondingly adjusted best-fit values without repeating the fit. It is also 
easy to infer from our results the quantitative consequences of a particular GUT model 
predicting the relationship between a s and sin 2 w {m z )y^. 

The overall fit to all Z parameters listed above gives the following result : 



g 2 (m 2 z ) = 0.5542 - 0.00030 ±0.0017 (l 0.14 -0.36 \ 

1 0.20 

V 1 / 



o.oi 

s 2 (m|) = 0.2313 + 0.00008 ±0.0007 p, 

5 b (m 2 z ) =-0.0061 - 0.00430 ± 0.0034 



(4.3a) 



9 /a s - 0.1029\ 2 . 

where the errors and correlations are nearly independent of a s . The above parametrization 
for the a s dependences of the mean values and Xmin are accurate interpolations of our 
fit results (at the level of 1%) in the range 0.09 < a s < 0.15. The bottom and charm 
quark masses were set to m;, = 4.7 GeV and m c = 1.4 GeV. A shift of the bottom mass 
by ±0.2 GeV implies only the fitted Sb{m 2 z ) value to be displaced by ±0.0002, which is 
negligibly small compared to its error (±0.0034). Similarly, shifting the charm quark mass 
by ±0.2 GeV does not affect the above results, as expected. In particular, in the favored 
range 0.11^a s (m z )<^0.13 the quality of the fit is good, e.g. Xmin = 4.6 at ot s = 0.12 for 
9 — 3 = 6 degrees of freedom. 

Results of the fit (|4.3|) are displayed in Fig. [1^ by the 1-cr allowed contours in the three 
projections, (s 2 (m|), g z (m 2 z )), (s 2 (m|), Sb(m 2 z )), and (<5{,(m|), g z (m 2 z )). The contours are 
shown for three representative a s -values : a s =0. 11 (dashed), 0.12 (solid) and 0.13 (dash- 
dotted). Also shown by the lattices are the SM predictions for m t = (100-200) GeV and 
m H =(50-1000) GeV. In these SM predictions, all known two-loop corrections of the 
0(mf) and 0(aa s ) level [[46| , p4] - |58| , |63]| are included, as explained in detail in appendix 
C. Hence, the predictions depend weakly on a s due to the 0(aa s ) corrections, as well 
as on 5 G and 5 a which are needed to predict the charge form factors from the known 
(a,GF,m z ) values. The SM predictions in Fig. [15| are calculated for a s =0.12, <f G =0.0055 
and S a =0. Changing a s by ±0.01 has little effect, but changing 5 a by ±0.10 leads to 
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a shift in the SM predictions for s 2 (m 2 z ) by +0.00026, which is as large as 40% of its 
uncertainty: see eq. ( [4.31) . 



It is clearly seen from the figure and eq. ([4.3|) that the weak mixing form factor 
s 2 (m|) is determined almost independently of a s , while the Z coupling strength g 2 z {Tn 2 z ) 
is anti-correlated with the assumed a s value as a reflection of its sensitivity to the total Z 
width. This anti-correlation leads in the SM to a preference of larger mt for smaller a s , 
since g\ (m%) grows with m t (see Fig. ||). Furthermore, the fitted 5b(m 2 z ) value depends 
strongly on the a s value assumed. The minimum of \ 2 is reached at a s = 0.1029 in 
eq. ( [4. 3 b| ), a value slightly outside the range 0.11^o; s (m z )^0.13 expected from various 



QCD analyses [[3(J . 

It is instructive to elucidate the properties of the fit to the Z parameters in three 
steps. First, the relatively small sensitivity of the parameter s 2 (m 2 z ) to a s can be under- 
stood easily, since it is derived essentially from the asymmetry parameters being either 
completely or nearly insensitive to QCD corrections. Indeed, the fitted values of s 2 (m 2 z ) 
as determined from each asymmetry measurement (see also Fig. |7|) turn out to be : 



s 2 (m 2 z ) 


= 0.2309 




±0.0010 


(from A%), 


(4.4a) 


s 2 (m 2 z ) 


= 0.2316 




±0.0018 


(from P T ), 


(4.4b) 


Ami) 


= 0.2365 




±0.0055 


(from< R ), 


(4.4c) 


s 2 (m 2 z ) 


= 0.2313 + 0.00004 


a s -0.12 
0.01 


±0.0011 


(from Ape), 


(4.4d) 


s 2 (m 2 z ) 


= 0.2302 + 0.00004 


a a -0.12 
0.01 


± 0.0035 


(from A%), 


(4.4e) 



almost independent of g z (m 2 z ) and Sb(m 2 z ). Note that although the quark (q = b,c) 
forward-backward asymmetries have mild a^-dependences due to the perturbative QCD 
corrections ||76|| , they still can be neglected compared to the experimental uncertainties. 
From the above asymmetry data alone one finds 

s 2 (m 2 z ) = 0.2312 ±0.0009 (from A%, P T , A° LR ), (4.5a) 



-2/ 2 

s (rriry 



0.2312 + 0.00002^^P ± 0.0007 (from A%, P T , A£ R , A ^, A%£). (4.5b) 



The precision of the above determination of s 2 (m 2 z ) from the asymmetry data alone is 
almost as good as that of the global fit to all the Z parameters. These asymmetry 
measurements are particularly important for GUT studies, since the parameter s 2 (m 2 z ) is 



directly related to the unifying coupling s 2 (/j) = sin w (fi)y^ via eq. (|2.12[) . 

Next, the best-fit value s 2 (m|) pa 0.2313 is taken to probe the sensitivity of the 
remaining four observables to the parameters gz{m 2 z ) and 5b(m 2 z ). As explained in sec- 



tion 3.1, three of the remaining four observables, T z , a\ and Re, are sensitive to the a 
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value assumed, but only through the combination <5 b (m|) + 0.6a s (|3.25|) , or equivalently 



a s + 1.65fe(m|). Y z is also sensitive to g%ijn 2 z ). Hence, a 2-parameter fit to the above 
three observables for s 2 (m%) = 0.2313 leads to : 



g\{m\) _ = 0.5547 ±0.0017 

a s + \M h {ml) = 0.106 ±0.007 



p corr = -0.46 . (4.6) 



The above result is found to be insensitive to the a s value in the range 0.10 < a s < 
0.14. The above result for g z (m 2 z ) is consistent with the global fit ( |4.3|) , as may be 
verified by evaluating g z (m 2 z ) at the minimum of \ 2 i a s — 0.1029). The anti-correlation 
above reflects the fact that V z remains unaltered, while increasing g z {m 2 z ) and decreasing 
a s + 1.65fe(m|) simultaneously. 

Only one Z observable is now left, namely Rb- In section 3.1 Rb was found to be 
sensitive to the parameter lj,(m|) alone. A 1-parameter fit to Rf, yields : 

5 b (m 2 z ) = 0.0012 ± 0.0068 , (4.7) 

keeping the other parameters fixed at s 2 {m 2 z ) = 0.2313, g 2 z {m 2 z ) = 0.5542 and a s = 0.12. 
However, this fit is insensitive to variations around the values of the fixed parameters. 
Note, the SM predicts a negative value of 8b(m 2 z ) for large m t (see Fig. [!]). Thus, there 
is poor agreement with the expected large m t behavior of the Zb L b L vertex correction 
from the present Rf, measurement alone. Since the parameter a s enters the fit only in the 
combination a s + 1.65b(m 2 z ), the fitted <5f,(m|) can be interpreted as a constraint of a s . 
From eqs. (|4.6j ) and (|4.7|) follows 



a s = 0.104 ±0.013. (4.8) 
The fit to the Z shape parameters with both 5b and a s left free yields : 

Pcorr = -0.85. (4.9) 



8 h (m%) = 0.0014 ±0.0070 
a s (m z ) = 0.103 ±0.013 



The large errors and the strong anti-correlation among them show that it makes little 
sense to extract a s model-independently from the electroweak experiments on the Z-pole, 
as also noted in ref. p2| . The low best-fit value of a s reflects essentially the actual value 
of R b , which is larger than the SM prediction in the range 150 GeV < m t < 200 GeV 
(see Fig. U). It is therefore necessary to assume the SM contributions to 5b(m 2 z ), and to a 
lesser extent those to g z {m 2 z ), in order to extract a s from the electroweak Z parameters. 
The result of such an analysis is given in section 5.4, where consequences of the minimal 
SM are studied. 
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Finally, we present the result of 1-parameter fits to four observables, T z , cr°, Re and 
R b , respectively, in terms of the parameter 5 b (m 2 z ), for various values of gz(m%), s 2 {m 2 z ) 
and a s . Here, we neglect correlations in the errors and find: 

8 b (m 2 z ) = -0.0068 - 0.0084 ^ ( "^ ) ~° 7 5542 + 0.00 20 s " 2(m ^J 7 2313 - 0.0061 ± 0.0077 

(from Tz) , (4.10a) 

5 b (m 2 z ) = -0.0210 + 0.0000 ^ (m ^°' 5542 + 0.00 04 §2(m ^~ ° 7 2313 - 0.00 63 ± 0.0168 

(from a° h ), (4.10b) 

5 b (m 2 z ) = -0.0078 + 0.0000 ^"f^ 5542 + 0.0011 ^f^ff^ - 0.0061 ± 0.0044 

(from R e ) , (4.10c) 

5 b (m%) = 0.0012 + 0.0000 ^^"y 5542 - 0.0001 s " 2(m ^~ ° 7 2313 - 0.0001 ± 0.0068 

(fromi^). (4.10d) 

The above fits clearly confirm quantitatively our observations that T z , a\ and Re measure 
the combination 5 b + 0.6 a s (|3.25|) , that T z is also sensitive to g z (m 2 z ), and that R b is 
sensitive only to 5 b {m 2 z ). At present the data Tz, 0"° and Re favor a negative 5 b (m 2 z ) 
value consistent with the SM prediction for 150 GeV < m t < 200 GeV, while R b data 
gives a S b (m z ) value consistent with zero, at a s m 0.12. The combination of all the above 
measurements together with all the asymmetry data, and properly accounting for the 
correlations in the errors, yields 

8 b (m 2 z ) = -0.0062 - 0.0014 ^^"y 5542 + 0.00 09 ' 2(m ^~ 07 2313 - 0.00 46 ± 0.0031, 

(4.11) 



in accordance with the result ( f4.3| ). Note that the coefficient in front of a s in eq. ( 4.11|) 



is smaller than 0.6 in the combination (|3.25|) as a consequence of including the additional 
information due to R b . 

4.2 Low energy neutral current experiments 

The two universal parameters s 2 (0) and g z {0) can be extracted from four types of 
low energy neutral current experiments : the neutrino- nucleon scattering (u^-q), the 
neutrino-electron scattering (u^—e), atomic parity violation (APV) and the polarized 
electron-deuteron scattering experiments (e-D). Effects due to small, but finite, momen- 
tum transfer in these processes are accounted for by assuming the running of these form 
factors to be governed by the SM particles only (see Fig. ^]), which, at low energies, is 
an excellent approximation. Vertex and box corrections are calculated by assuming that 
they are dominated by the SM contributions. For details of the theoretical predictions, 
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see section 3.2. For each sector, first a model-independent parametrization of the data is 
given, and then the fit result in the (s 2 (0), gf (0)) plane. 



4.2.1 Neutral currents in v^—q scattering 



For the v^-q data, the results of the analysis of ref. |7§] are adopted. In terms of the 
model-independent parameters (fl , £,fl l fl,5£,<^), the following fit has been obtained: 



9l 
9% 
Si 

si 



0.2982 - 0.0058(m c - 1.5) ] ± 0.0028 ± 0.0029, 
0.0309 - 0.0053(m c - 1.5) ] ± 0.0034 ± 0.0028, 
-0.0588 - 0.0025(m c - 1.5) ] ± 0.0233 ± 0.0042, 
0.0206 + 0.0010(m c - 1.5) ] ± 0.0155 ± 0.0039, 



(4.12) 



where the former and the latter errors denote the experimental and the parametrization 
errors. The correlation matrices for the two types of uncertainties also quoted in ref. 
are respectively 



Z 1 



r corr 



-0.751 
1 



-0.100 
0.064 
1 



V 



0.118 ^ 
0.097 
-0.436 
1 / 



n (p ar ) 

r corr 



0.606 ^ 
-0.677 
-0.712 

1 / 
(4.13) 

The fitted parameters depend on the assumed value of the charm quark mass (m c in GeV 
units) PB[ entering the slow-rescaling formula [9H] for the charged current cross sections. 



( 1 -0.914 
1 

V 



-0.975 
0.945 
1 



The data [78[ constrain the charm quark mass to 



m c = 1.54 ±0.33 GeV. 



(4.14) 



After summing the experimental and the parametrization errors in quadrature, and inte- 
grating out the m c dependence of the above parametrization under the constraint ( f4.14| ) 
the new model-independent parametrization of the v^-q data gives : 



9l 
9l 
Si 

s% 



0.2980 ± 0.0044 
0.0307 ±0.0047 
-0.0589 ± 0.0237 
0.0206 ±0.0160 



\ 



-0.559 -0.163 
1 0.156 
1 



0.162 ^ 
-0.037 
-0.447 
1 / 



(4.15) 



which properly accounts for the uncertainty in m c . The parametrization Q4.15 ) serves as 
input to our analysis. 

By using the theoretical formulae ( 3.44 ) and ( |3.45|) of section 3.2.1 the data ( 4.15| ) can 
be confronted with the predictions in terms of s 2 (0) and (?|(0). Corrections due to small, 
but finite, momentum transfer are evaluated at 



(-t)n.c. = H) c.c. = 20 GeV 2 , 



(4.16) 
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in eqs. (|3.47| ) and (|3.51|) and in the running of s 2 (t): see eq. ( |3.56|) . The fit result is : 

gUo) = o.55oo±8:88S 
,- 2 (o) = o.242ot8:SS 



p corr = 0.916, (4.17a) 

= u.242U_5;5ii2 j 
xLn = 0-13, (4.17b) 



Asymmetric errors are quoted. The non-gaussian behaviour of the x 2 function reflects the 
non-linear transformation between the charge form factors (<?f(0), s 2 (0)) and the model- 
independent parameters (g 2 Ll g\, 8\, 5^), as seen in Fig. [Hj. The strong positive correlation 



between the fitted values of s (0) and gf (0) is a consequence of eq. (|4.15|) : the precisely 



measured combination g\ ± g 2 R in ( |4.15| ) dominates the total neutral current cross section 



off isoscalar targets. The 1-a contour of the above fit is shown in Fig. [16[ It can be 
reproduced rather accurately by the following parametrization: 

g 2 (0) = 0.5497 ±0.0080 1 

s 2 (0) = 0.2413 ±0.0136 J Pcorr ~ U " J °' 1 8 J 

xLn = 0-13, (4.18b) 



which serves merely for estimating the constraints from the v^— q experiments. We stress 
that all the quantitative analyses in the following sections are performed by fitting directly 
to the original parametrization of the data ( |4.15|) . 



4.2.2 Neutral currents in v^—e scattering 



The u^-e data from the three experiments: CHARM, BNL E374 and CHARM-II |84 



have been summarized in ref. J79] in terms of the model-independent parameters s ve and 

Pve- 

n,.,~(nYt e = 1.007 ±0.028 1 

p corr = 0.09. (4.19) 



Pve = {p)lt = 1-007 ±0.028 

s 2 ue = (sin 2 9 W ) V £ = 0.233 ± 0.008 



As explained in detail in section 3.2.2, first the total cross section a ve and a ve is recon- 
structed by using the formula (|3.75| ), and then the fit is performed by using the theoretical 
expressions ( j3.64j ). The reconstructed cross sections are found to be 



a I/e /^(10- 42 cm 2 /GeV) = 1.56 ±0.10 
a P 7E p (10- 42 cm 2 /GeV) = 1.36 ± 0.09 



p corr = 0.51 , (4.20) 



and the 2-parameter fit to the above data gives 

Pcorr = 0.09. (4.21) 



g|(0) = 0.5459 ± 0.0154 
s 2 (0) = 0.2416 ± 0.0079 



The same result follows if we use the approximation (|3.74|) directly to fit the parametriza- 



tion ([4.19|) . Here Xmin = 0, since the number of degrees is 2 — 2 = 0. The result is 
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shown in Fig. 16 by the l-a contour. The weak mixing form factor s 2 (0) is measured 
more accurately in the u^—e experiments than in the v^-q experiments, whereas for g\ (0) 
it is the other way around. 



4.2.3 Atomic parity violation 

As for the APV experiments the result of the analysis 
transitions in the cesium atom (A, Z) = (133, 55) are used : 

Qvk(135,55) = -71.04 ± 1.81. 



on the parity violating 



(4.22) 



The quoted uncertainty is the quadratic sum of experimental and theoretical errors. After 
correcting for the vertex and box corrections [^5j as explained in detail in section 3.2.3, 
one finds 

s 2 (0) = 0.2294 - 0.6178 [f z (0) - 0.5486] ± 0.0082. (4.23) 

Here the value g|(0)=0.5486 stands for the prediction at p—1 or T = S G /a. The result is 
shown in Fig. [16] by l-a contours. As anticipated in the previous section, the correlation 
between the fitted s 2 (0) and gf (0) values is opposite to that from v—q fit. As a conse- 
quence, the constraints on both s 2 (0) and g§(0) are improved significantly by combining 
the two types of experiments. 

4.2.4 Polarization asymmetry in e— D scattering 



Finally, for the SLAC eD polarization asymmetry experiment |87[ a model-independent 
fit is performed to the original data by using the two combinations, 2C\ U — and 
2C*2« — Cid of the coefficients of the effective parity violating e-q neutral current operators 
24j| : see eq. fl3.76|) . In the quark parton model with the valence quark approximation the 



observed polarization asymmetry is expressed in terms of the above parameters by 

6 Gf f,_ „ „ , „ „ , l-(l-y) 2 



A 
~Q 2 



5v^e 2 (-Q 2 ) 



(2 Ciu — Cid) + (2 C 2u — CW) 



(4.24) 



which depends on the scaling variable y, but not on x. The mild Q 2 dependence due 
to the running of the effective QED charge e 2 (— Q 2 ) is accounted for. There have been 



extensive studies [^,|96|, which show that the above approximation is in fact valid on more 
general grounds, but that it may suffer from higher-twist contributions. We therefore 
perform a new model-independent fit to the original data [S7|, and obtain quantitatively 
the theoretical uncertainty in the fitted parameters. 

By taking account of the sea-quark contributions and finite R = <tl/ '<Jt HU; as well as 
possible higher twist contributions f96| , |97l ) the above simple expression for the asymmetry 
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Q4.24Q is modified as follows 
A 6G 



with 

c = 1.345 - e s /5 . (4.27) 

Here the x-dependent parameters e u and e s denote the relative contribution of the sea 
u-quark and that of s and s quarks, respectively, which are parametrized by 

(1-x) 4 

e u = e ' , (4.28) 
\/x 



The uncertainty in the factor e above is estimated to be 

e = 0.1 ±0.03. (4.30) 
The effects of introducing sea-quark contributions in the fit is shown in Fig. |T^(a). As 



found in ref. [p5j, the effect is very small along the tree level SM prediction as shown 
by the straight line in the figure. Some representative values of sin 2 6 W in the SM are 
denoted by blobs. The longitudinal to transverse cross section ratio R = ctl/ct is allowed 
to vary within the rather conservative limits 

R = 0.2 ±0.2. (4.31) 

The effect of introducing the R parameter alone is shown in Fig. 0(b) and the result 
turns out to be insensitive to its uncertainty, especially along the tree-level SM trajectory, 
confirming the earlier observation of ref. [^] . Finally, the parameter 5 in the factors b and 
c parametrizes the higher twist effects as expected in the MIT bag model P7| . Taking 
as the magnitude of the uncertainty the largest value of the MIT bag model estimate of 
ref. Ml yields 



5= (1.58 ±1.58) x 10~ 3 . (4.32) 

The effects of introducing the S parameter alone are shown in Fig. |TT|(c). As in the case 
of the sea-quark contributions (Fig. |l7Ka)), the effect is negligibly small along the line 
representing the tree-level SM prediction. Note that the higher-twist effects are found to 



53 



be rather model dependent |98| . The MIT-Bag model estimates [|)7| adopted here lead 
to quite small corrections, as in the neutrino scattering off isoscalar targets [pSR . Further 
study on the higher twist effects may be needed to achieve precision measurements of the 
electroweak parameters in these reactions. 

After allowing for all of the above uncertainties, one finds 

2C lu -C ld = +0.94 ± 0.26 \ _ . . 

2C 2u -C 2d = -0.66 ±1.23 J Pcorr ~ 1 j 

with Xmin = 9-95 for 11 data points, that is, a good fit. The above result is shown 
in Fig. |H|(d). Because of the strong correlation, only a linear combination of the two 
coupling factors is measured well. 

By using the theoretical formulae ( |3.78| ), a fit is made to the data ( [4.33| ) in terms of 
the two parameters s 2 (0) and g\ (0). In order to fix the g 2 -dependent factors (Q 2 = — q 2 ) 
ri(-Q 2 ), r 2 (-<5 2 ) and s 2 (-Q 2 ) in the amplitudes we choose (Q 2 ) = 1.5 GeV 2 . Note, 
however, that Q 2 -dependence of each data point |TJJ and that of the QED running coupling 
e 2 (— Q 2 ) in eq. ( |4.25|) have been respected in the model-independent fit (|4.33|) . The result 
is : 

s 2 (0) = 0.2273 + 0.3067 [#f(0) - 0.5486 ] ± 0.0092, (4.34a) 
xLn = 0-46 - 1.77 [£§(0) - 0.5486 ] , (4.34b) 



and shown in Fig. [16]. Note that the parametrization ( f4.34b| ) is valid only in the vicinity 



of the SM predictions gf(0) ~ 0.55 (but is valid in the whole region of Fig. [16]), and that 
the global Xmin ls zer ° ; since the two parameter parametrization ( |4. 33| ) is adopted as the 
original data of our fit. 

4.2.5 Summary of low energy neutral current experiments 

In this section the fits to the electroweak observables in the four low energy neutral 



current experiments are summarized. The fit results are illustrated in Fig. [16] by 1- 
a allowed regions in the (s 2 (0), (?|(0)) plane. Since all four pieces of information are 
consistent with each other, a combined fit can been performed : 



Pcorr = 0.53, (4.35a) 



f z (0) = 0.5462 ±0.0036 ' 
s 2 (0) = 0.2353 ±0.0044 _ 

X 2 min = 2.22. (4.35b) 



The fit with 7 = 9 — 2 degrees of freedom is good and its result is shown in Fig. [16] by the 
ellipse with the thick 1-cx contour. 
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It is sometimes useful to analyse the neutral current sector with and without inclusion 
of the neutrino data, since in some models they receive different new physics contributions. 
To this end the fit is done separately for u^—q and u^—e experiments : 



Pcorr = 0.75, (4.36a) 



f z (0) = 0.5496 ±0.0068 ' 
s 2 (0) = 0.2414 ±0.0047 

xLn = 0-19. (4.36b) 



The fit for the APV and eD experiments gives : 



-0.62, (4.37a) 



£|(0) = 0.5510 ±0.0165 1 
s 2 (0) = 0.2280 ±0.0088 J Pcorr 

Xl in = 0.46. (4.37b) 



These two fits are again consistent and their combination reproduces, of course, the above 
global fit flPSD . 



4.3 Charged current experiments 

The W mass measurements have been updated recently by the CDF and DO collab- 



orations. By combining the most recent measurements |88|] and the previous result of 
PDG |25|] one obtains 

m w = 80.24 ±0.16 GeV. (4.38) 
Note that in this analysis the W mass definition follows the LEP convention ||, as opposed 



to the pole mass definition: see eq. ( |2.14 ). The pole mass should be smaller by about 



0.03 GeV. The difference is still negligibly small as compared to the error of 0.16 GeV. It 
is worth noting that the W propagator with running width factor gives a more accurate 
description of the scattering amplitudes when no imaginary parts are introduced outside 
the propagator factor. 

The electroweak parameter ^(0) is now obtained by combining the m w measurement 
with the fi life-time parameter Gp ( |3.99| ) : we find 

9 , , 6 r - 0.0055 
f w (fi) = 0.4225 - 0.0031 ± 0.0017 , (4.39) 

a 

where 6 G = 0.0055 is the SM estimate for the process specific correction to the \i life-time: 
see eq. fl3.98p . No other experiment in the charged current sector is accurate enough to 
provide adequate information for our electroweak analysis. Precise measurements of the 
W shape parameters [|100|1 would improve our knowledge in this sector considerably. 
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5 Systematic analysis 

In this section, first the g 2 -dependence of the two charge form factors g z (q 2 ) and 
s 2 (q 2 ) is examined between q 2 = and q 2 = m 2 z . Next a combined fit in terms of the 
S, T and U parameters is made assuming the g 2 -dependence of these charge form factors 
to be governed by the SM. Finally, only the SM particles are assumed to contribute to 
the radiative effects and the preferred range of the two mass parameters m t and m H is 
searched for. Also the a s and 5 a dependences of the fits are discussed in detail. 

5.1 Summary of all experimental constraints on the electroweak parameters 

The information on all electroweak precision data has been represented in the previ- 
ous sections in terms of the charge form factor values (see eqs. (|4.3| ), (|4.35|) and ( |4.39| )) 



and is, for convenience, collected in Table || In addition, the fine structure constant a 
determining the charge form factor e 2 (0) = 4na (see Tables |I] and 0) has been used as 
an input parameter. In calculating x 2 the model-independent parametrizations of the 
original data are used as inputs for the fit : eqs. (|Q| )-(|4^) for the Z parameters (sec- 
tion 4.1), eq. (|4.15 ) for the v^-q scattering experiments (section 4.2.1), eq. ( 4.19|) for the 



v, 



t -e scattering experiments (section 4.2.2), eq. (|4.22|) for the atomic parity violation ex- 



periments (section 4.2.3), eq. ( f4.33|) for the e-D polarization asymmetry measurements 
(section 4.2.3), and eq. ( |4.38| ) for the W mass measurements (section 4.3). The \ 2 fits 
in each of the various sectors look all fine and it is concluded that the whole body of 
data is consistent with the assumption of the SU(2)l x U(1)y universality and the SM 
dominance of the vertex and box corrections. 

5.2 Testing the running of the charge form factors 

If there are new particles coupled to the weak gauge bosons with masses near or 
below m w and m z , their signal can be identified as an anomalous running of the charge 



form factors [11, 12]. In principle, the running of all four charge form factors provides 



us with information on new physics contributions via eq. ( |2.30| ) for l/a(q 2 ), eq. ( |2.40a| ) 



for s 2 (q 2 )/a(q 2 ), eq. ([2.40 bp for l/g z (q 2 ) and via eq. ( |2.4Ucj ) for \jg 2 v {q 2 ). At present, 
only two of the four form factors, s 2 (q 2 ) and cjz{q 2 ), have been determined with sufficient 
accuracy at two different energy scales, q 2 = and m 2 z . 
The results collected in Table |] yield : 

47T All 



9 2 z«) 9 2 z(0) 



-0.33 + 1.2(a s - 0.12) ±0.17 



772/ 



s"(m" 7 ) s 2 (0) 

_/ 9 f ; — = -2.47+1.1 (a s - 0.12 ±0.62 



p corr = -0.49. (5.1) 
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In the absence of a precise value for a(m|) the SM prediction a(m|)sM 
more generally, 5 a = l/a{m 2 z ) — 128.72 = 0) is used above. 

Fig. [18| illustrates SM running of the charge form factor g z (q 2 ), 

4:71 47X 1 

4 



1/128.72 (or, 



9z( m z) 



3f(o) 



Sz(m 2 z ) - S z (0) 



(5.2) 



as a function of m H , together with the experimental constraint ( p.l[ ). The g 2 -dependent 
Sz function is defined in terms of the gauge boson two-point functions in eq. ( |B.41a| ) of 
appendix B. The difference (|5.2| ) takes the form (|2.40b|) of section 2. The m t dependence 
of the SM prediction is very small compared to the experimental error for m t > 100 GeV. 
The SM is consistent with the data as long as the Higgs boson mass is not too small. 
Note that the 1-cx constraint on m H , m H > 2.9 GeV (67%CL), is obtained merely by 
comparing the Z boson coupling strengths at q 2 =0 and q 2 = m 2 z . These values are, 
however, obtained by neglecting the Z — > Hff contribution to T z , and are anyway 
excluded by direct searches at LEP {m H > 63 GeV) ||1 1 1 . 

The Higgs boson does not contribute to the running of the other neutral current form 
factors, l/d(q 2 ) and s 2 (q 2 )/a(q 2 ). They are affected by loops of charged particles only, 
and, for instance, the top quark contributions to the running of these form factors are 
parametrized in appendix B, in eqs. ( |B.27] ) and ( |B.28| ): 

1 



mm 



128.71 + £ had + 0.024 1 + 5 



?:2, 



ffl 



Z)SM 

s 2 (0) 



ayrrir 



a 



SM 



-3.09+^ + 0.009(1 + 5 



7i 
71 



100 GeV \ 2 

m t ) 
100GeV \ 2 
m t J 



(5.3a) 
(5.3b) 



The mt dependences of these runnings are very small for m^lOO GeV. 

The running may be appreciable, if there is a charged fermion with mass near to 
half the Z mass [12|. The case of a light wino, the fermionic partner of the W in the 
super symmetric SM, is shown in Fig. [H| (a) 47r/gf(m|) — 47r/g!(0) , (b) s 2 (m 2 z ) / a{m 2 z ) — 



s 2 (0)/a , and (c) 5 a = l/a(m 2 z ) — 128.72 . The singularity at m v 



m z /2 of the charge 



form factor 4u / g 2 z (m 2 z ) in (a) reflects [^] the deviation of the Z line-shape from the 
Breit-Wigner form assumed both in the experimental fit and the corresponding theoretical 
formulae, and is unphysical. The 1-cx bound on the wino mass, m wino > 46.1 GeV, as read 
off from Fig. ^ is unrealistic, since the threshold 2m wino = 92.2 GeV is less than a half 
width away from the Z-pole. In order to derive constraints on particles very near to 
the threshold, one should look for deviations of the Z line shape from the simple Breit- 
Wigner form [[H], ^j5j. When calculating the predictions for (b) and (c) the hadronic 
vacuum polarization contribution to the running of these form factors is set to <5h a d — 0, 
while the present estimate [^8j is <5h a d = + 0.1 (|B.22| ). Wino of masses around 50 GeV 
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may shift 5 a = l/a(m 2 z ) — 128.72 from its canonical value 5 a = by about 0.1, which is 
of the same order as the present uncertainty in the SM prediction. 

It is clearly seen from Fig. 18 and from eq. (5J5) that the results ( |5.1|) are consistent 
with the SM predictions in the range m H > 60 GeV and m t > 100 GeV. The study of 
the two examples, a very light Higgs boson and a supersymmetric wino, demonstrates 
that more accurate values of s 2 (0) and g z (0) are needed to detect effects of new physics 
through the running of the charge form factors. Accurate measurements of the charge 
form factor l/a(q 2 ) at \q 2 \ ~ m| should also provide independent information. 

Fig. 2C shows the above results in the (s 2 (m 2 z ), g 2 z{ m2 z)) plane, where the Z parameter 
fit ('LEP+SLC') is taken from Fig. [15| for a s = 0.12, and the combined low energy fit 



of Fig. [II] has been rescaled to the m z scale by assuming SM running of the two charge 
form factors, s 2 (q 2 ) and g z (q 2 )- The combined low energy neutral current data (see 
Table ^| and Fig. [IB]) are displayed for various choices of m t and m H in order to put in 
evidence their small, but finite, effects on the running of these form factors. The four 
contours are obtained for m t =100, 200 GeV and m H =60, 1000 GeV. At m t = 175 GeV 
and m H = 100 GeV, the fit ( |4.35a| ) for the low energy neutral current data can be re- 
parametrized as 

g\(m\) = 0.5533 ±0.0037 ' 
s 2 (ml) = 0.2266 ±0.0047 



Pcorr = 0.53. (5.4) 



It is seen from the figure that the low energy neutral current fit and the Z parameter fit 
in terms of the charge form factors s 2 (q 2 ) and g z {q 2 ) are in accordance with the running 
of these form factors as predicted by the SM. 

The thick solid contour marks the result of the fit to all neutral current experiments 
as summarized in Table || assuming the SM for the running of g z (q 2 ) and s 2 (q 2 ) : 



g 2 (m 2 z ) = 0.5544 - 0.00023^^ ±0.0015 
s 2 (m|) = 0.2312 + 0.00008^P ± 0.0007 p c 
5 b (m 2 z ) = -0.0064 - 0.00437^^12 ± 0.0034 



o.oi 

m., - 0.1024\ 2 



/ 1 0.16 

1 0.20 | , (5.5a) 




*U = 4.67+1 0.0127 J' (5 ' 5b) 
In the global fit the uncertainty due to m t and m H in the running of the form factors is 
negligible in the range m t =100-200 GeV and m i? =60-1000 GeV. The Xmin vam e of 6.6 
for a s = 0.12 is acceptable for 15 (=18 — 3) degrees of freedom. In conclusion, there is no 
indication of new particles with mass ^m z in the running of the charge form factors. 

Note that the errors in (|5.1[ ) are dominated by those of the low energy experiments. 
Further improvements in the low energy precision experiments are required to detect a 
signal of relatively light new particles, should they exist, through anomalous running of 
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the charge form factors. In comparing the global fit of Fig. with the individual fit to 
the low energy NC data in Fig. 16, the fit from the v^-q and v^—e experiments (|4.36|) 
are remarkably consistent with the Z parameter fit of Table || whereas the fit of the e—q 



sector ( 4.37 ) based on the APV and the e-D asymmetry measurements are about 1.5 



standard deviations away. For m t = 175 GeV and m H = 100 GeV, the fit ( 4.36a|) for u^—q 
and Vn-e scattering is re-parametrized as 



g%{m 2 z ) = 0.5568 ± 0.0048 
s 2 (m 2 z ) 



Pc 



0.75, 



0.2331 ±0.0072 # 

while the fit (|4.37a|) for the APV and polarized e-D experiments gives 



(5.6) 



g\{m\) = 0.5583 ±0.0170 



s 2 (m 2 z ) 



0.2188 ±0.0093 



Pc 



-0.62. 



(5.7) 



Further studies of polarization asymmetries in the e—q sector as well as studies of the 
neutral current processes at TRISTAN energies might be potentially rewarding. 



5.3 Testing the 3 parameter universality 

Once the g 2 -dependence of the charge form factors is assumed to be governed by 
SM physics alone, all radiative effects to the gauge bosons depend on three universal 
parameters : S, T, U . They include the SM radiative effects as well as new physics 
contributions, as opposed to the original definitions of ref. [|J. While the charge form 
factors gz(m%), s 2 (m 2 z ), g^r(0) can be directly confronted with experiments, the S,T,U 
parameter fit suffers from uncertainty in the QED effective coupling a(m 2 z ), the reason 
being the fact that the charge form factors are determined by the S, T, U parameters under 
the (a, Gf, m z ) constraints (see discussion in section 2.3). The magnitude of a(m|) is 
controlled by the external parameter S a = l/a(m 2 z ) — 128.72. 

A 4-parameter fit yields : 



S = -0.35 
T = 0.39 
U = 0.41 



-0 016 a *-°- 12 
o.oi 



-0.058 ■ 



-0.12 



0.01 

±0.058 
0.0064 -0.0043 



±0.067 ^±0.33 
-0.004 ^ ±0.36 



/ 1 



±0.024 ^±0.54 

±0.0034 



0.83 
1 



-0.18 -0.12 \ 
-0.40 -0.32 



1 



0.20 
1 



(5.8a) 



from the result of the global fit summarized in Table || The best-fit values of S, T, U 
and 5b are weakly dependent upon ot s and 5 a as quoted explicitly in eq. (|5.8a|) . The 
minimum of x 2 turns out to be practically independent of 5 a . We therefore add to the fit 
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the independent knowledge 5 a = 0.0 ±0.1 [|^] leading then to : 

X - = 467 + I 0.0127 ) + loloj • (5 ' 8b) 

The correlation between S and T is strong, since they are constrained by the precisely 
measured weak mixing form factor s 2 (m|) via eq. ( |2.38b]) . 

The above results are shown in Fig. |21]by l-er contours as projections onto the (S, T), 
(S, U), and (U, T) planes. The contours are drawn for three a s values, a s =0.11 (dashed 
lines), 0.12 (solid lines) and 0.13 (dash-dotted lines), and for m t = 150,200GeV and 
m H = 100, 1000 GeV in the running of the charge form factors s 2 (q 2 ) and g 2 z {q 2 ) between 
q 2 = and q 2 = m|. The fit results depend slightly on m t and m H in the above range. The 
numerical values of eq. ( p.8| ) are obtained for m t = 175 GeV and m H = 100 GeV. The SM 
predictions of appendix C are drawn in Fig. [H] by lattices in the region m t =100-200 GeV 
and m^=50-1000 GeV. 

The fitted T parameter depends only slightly on a s , when the parameter 5b is allowed 
to vary freely within the experimental constraints. If we fix 5b by a theoretical model, 
then the T parameter should have stronger a s dependence due to the correlation —0.31 
between the errors of T and 5b (see section 6.3 for more discussions). The S parameter 
depends on 5 a . The fitted S value is shifted by about 0.07 (that is, 20% of its present 
uncertainty of 0.33) for [5 q ]sm ~ <5had = ± 0.10. 

The parameters S, T and U measure electroweak radiative effects in the gauge boson 
propagators. The fit ( |5.8|) shows that the data favor negative S and positive T at a s = 
0.12 and 5 a = 0. The point S = T = U = 5b = 0, which represent no electroweak 
radiative effects in the gauge boson propagators nor in the Zb L b L vertex, is about 4.5 
standard deviations away from the minimal for a s = 0.12 and 5 a = 0. However, if 
in addition the electroweak radiative effects are dropped in the muon decay by setting 
5 G = in eq. ( |2.26| ), then according to the substitution rule ( |2.39| ) the 'no-radiative 
effects' point becomes T = 0.0055/a = 0.75, S = U = 5b = in the fit (J5.8|), which is 
only 2.6 standard deviations away from the minimal. Although this result still assumes 
the SM radiative corrections for the remaining vertex/box corrections, it is essentially 
the mechanism that led the authors of ref. |64j] to state that there had not yet been an 
evidence for genuine electroweak radiative effects. Our analysis makes it clear that it 
is more natural to interpret significant radiative effects in the T parameter which are 
approximately cancelled by the effect of the radiative effect 5 G in the prediction of the 
electroweak observables. 

The resulting x min of eq. (|5.8b| ) agrees nearly with that of eq. ( E>.5b|) . The effective 
number of degrees of freedom is in both cases 15, namely 19 — 4 respectively 18 — 3. The fit 
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to the NC data contains actually only three parameters, S, T and 5b(m 2 z ), corresponding 
to the charge form factors s 2 (m|), g z {rn z ) an d Sb(m z ) in the global fit. The present fit 
depends in addition upon U, when the charged current data (and thus the forth form 
factor, (7^/(0)) are included. 

5.4 Testing the Minimal Standard Model 

In the minimal SM, all the parameters g|(m|), s 2 (m|), <?f(0), s 2 (0), ^^(0) and 
Sb(m 2 z ) depend uniquely upon the two mass parameters m t and m H . Consequently, the 
results of the fits summarized in Table |6| are constraining m t and m H . We should repeat 
here that the SM contributions from the top-bottom doublet to the form factors are 
calculated by using the simple 0(aa s ) two-loop formula p4|-p6[1. Non-perturbative tt 
threshold effects [|60| -p2| will affect these corrections and the predicted m t value will shift 



upwards by as much as a few GeV ]B2] from the effect in the T parameter. Our approach 
separates clearly the data analysis in terms of the generic form factors and the analysis 
of the SM contributions to these form factors. Uncertainties in the latter process can 
hence be studied separately. In fact if the SM mt-dependence of the fit is dictated by 
the mt-dependence of the T parameter alone, then the sole effect of the non-perturbative 
threshold corrections can be expressed as a rescaling of the m t parameter in the following 
analysis. 

Fig. ^ shows the result of the global SM fit to all electroweak data in the (m H , m t ) 
plane ||102| , [103|1 for three representative a s values. The "x" indicate the minimum of x 2 j 
7.4, 6.6, 10.3 for a s = 0.11, 0.12, 0.13, respectively, the inner contours correspond to 1-a, 
the outer to x 2 — Xmin + 4-61 (that is, 90% CL). Dashed lines show the best m t values 
for a given m H . Note the positive correlation between the preferred values of mt and 
m H , which is found to be independent of the assumed a s value. On the other hand, the 
preferred range of m H depends rather sensitively on a s . For the cases a s {m z ) = 0.11 and 
0.12 smaller m H values are preferred, whereas for a s (m z ) = 0.13 larger m H is slightly 
favored. If the lower bound for m H , m H > 63 GeV at 95% CL measured by the LEP 
experiments | |101| |, is imposed, m t below 100 GeV is clearly disfavored for all a s , in 



agreement with the directly established lower top mass limit [|104| , |105| . 

The x 2 function in the global fit to all electroweak data can be represented in terms of 
the four parameters m t , m H , ot s {m z ) and 5 a together with the constraint S a = 0.0±0.1 
by: 

2 / m t — (m t ) \ 2 2 
X SM (m t , m H , oc s , 6 a ) = — ^ J + X H { m H, a s, <U , (5.9a) 
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where 



and 



, rriTT , , itltt (ol s — 0.12\ / 5 a \ 

(m t ) = 145.2 + 12.5 In — ^ + 0.9 n 2 — ^ - 1.9 — - 4.6 —2- 5.9b 

v 7 ioo ioo v o.oi / vo.ioy ' v ; 

Am, = 14.6 - 0.23 In ^ - fo.38 - 0.05 In m " 150 , (5.9c) 

100 v 100' 10 v ; 



2 , , , fi11 , /5 Q -0.31\ 2 /a, - 0.1173 + 0.005 5 Q \ 2 

0.0136 J 100 V 0.0700 J 100 V0.10 



a s - 0.1244 + 0.025 5 a \ , m H /a, - 0.1322 \ 2 m H . / 5 Q 



(5.9d) 

Here mt and m H are measured in GeV. This parametrization reproduces the exact \ 2 
within a few % accuracy in the range 100 GeV < m t < 250 GeV, 60 GeV < m H < 
1000 GeV and 0.10 < a s (m z ) < 0.13. The best-fit value of m t for a given set of m H , a s 
and S a is readily obtained from eq. ( |5.9b|) with its approximate error of (|5.9c| ), mutatis 
mutandis for m H . Due to the quadratic form it is easy to get the a s or 5 a independent 
results. Also additional constraints on the external parameters a s and S a , such as those 
from their improved measurements, can be discussed without difficulty. As explained in 
section 4.1, the SM does not fit well the ratio Rb- If we remove from our global fit the 
data on Rb, we find that the best-fit m t value above becomes larger by 3.9 GeV, almost 
independent of m H and a s , and the Xmin decreases by 2.4. 

Fig. |23] displays the overall \ 2 of the SM fit, XsM' as function of m t for m H = 60, 300, 
1000 GeV and a s (m z ) = 0.11, 0.12, 0.13. Also the uncertainty due to 5 a is shown for 
three cases, 5 a = —0.1 (a), (b), +0.1 (c). The results of the parametrization eq. (|5.9|) 
is shown by the dotted line. It is remarkable to see that the present knowledge of 5 a 
to ±0.10 affects the best-fit value of m t by about 5 GeV, while the uncertainty in a s 
of ±0.01 affects it by about 2 GeV. This observation emphasizes the importance of the 
asymmetry measurements for the prediction of m t through s 2 (m z ), where the dependence 
on S a in the SM is not negligible: see eq. ( |2.38b|) . On the other hand, the a s -dependence 



of the fitted m t comes from the constraint due to the Z total width, T z , which in turn is 
sensitive to m t mainly through g z {m 2 z ). We come back to this point in the next section 



when discussing the new left-right asymmetry measurement ||31|| . 

In Fig. [3] the overall \ 2 is plotted as functions of m H for m t =120, 140, 160, 180, 



200 GeV and a s (m z ) = 0.11, 0.12, 0.13 setting 5 a = 0. The dotted lines indicate our 
approximation Xsm °f (|5.9|). Obviously, the best-fit value of m H depends very sensitively 
on the m t and a s values. A small value of the Higgs mass is favored for m t < 140 GeV, 
values of a few hundred GeV for m t around 160 GeV and large values for m t > 180 GeV. 
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The preference of lighter m H is more pronounced for small a s , while heavier m H for larger 
a s . However, the m H dependence of x 2 is very mild and meaningful upper bound on m H 
can only be obtained for small a s and small m t . The upper and lower bounds on m H 
will be discussed more quantitatively in section 6 after inclusion of the new left-right 



asymmetry data [31 



For given m t and m H the QCD coupling a s (m z ) may be extracted within the SM 
from the electroweak data alone with the result : 



a s = (a,) ± 0.0060 , 
(a s ) = 0.1165 



m 



H 



0.00085 (— ) 2 + 0.00031 (In 

U00' V 100 



+ 2.6 - 0.0006 



6 n 



0.10 



(5.10a) 
(5.10b) 



where m t and m H are measured in GeV. The above parametrization reproduces well 
the a s dependence of the \ 2 function (|5T9|) in the range 100 GeV < m t < 200 GeV and 



60 GeV < m H < 1000 GeV. The error on a s determined from the electroweak data is 
found to be approximately 0.0060, almost independently of the assumed m t , m H and S a , 
while the mean value (a s ) is slightly sensitive to them; 



(a s ) 



(150, 60) GeV 
(175, 60) GeV 
(200, 60) GeV 
(150, 1000) GeV 
(175, 1000) GeV 
(200, 1000) GeV 



(5.11) 



0.1159 for (m t , m H ) 

0.1153 for (m t , m H ) 

0.1145 for (m t , w, H ) 

0.1220 for (mt, m H ) 

0.1214 for (m t , m H ) 

0.1206 for (m t , m H ) 

for 5 a = 0. There is a tendency in the SM fit to prefer larger a s for larger m H . 

Furthermore, if all radiative effects are assumed to be dominated by the SM con- 
tributions, the present electroweak data have some sensitivity to the parameter 5 a = 
l/a(m|) - 128.72. By excluding the last term in eq. ( |5~9dT ), (5 a /0.1) 2 the elec- 

troweak data alone provide the constraint : 



5 a 



(6 a ) ±0.24 
0.010 



m t — 150 , m 

0.139 — ~~rr. + 0.246 In 



,/ a s -0.12 
11 - 0.112 - 



(5.12a) 
(5.12b) 



10 100 0.01 

where m t and m H are measured in GeV. The above parametrization is valid in the range 
120 GeV <m t < 200 GeV, 60 GeV < m H < 1000 GeV and 0.11 < a s < 0.13. For some 
representative m t and m H values the exact evaluation of the x 2 function leads to : 



<*«> 



-0.09 for 

-0.45 for 

-0.87 for 

0.59 for 

0.25 for 

-0.12 for 



(m t , m H ) 
(m t , m H ) 
(m t , m H ) 
(m t , m H ) 
(m t , m H ) 
(m t , m H ) 



(150, 60) GeV 
(175, 60) GeV 
(200, 60) GeV 
(150, 1000) GeV 
(175, 1000) GeV 
(200, 1000) GeV 



(5.13) 
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for a s =0.12. The above fit is consistent with the direct measurement [<5 a ]sM ~ <^had 



± 0.10 when m t and m H are in the preferred range in Fig. g^. This confirms the 
importance of the direct <5h a d measurement in constraining the model parameters from 
the electroweak precision measurements. 

Considering the Xmin P er degree of freedom (see parametrization ( |5.9|) and Figs. f22H2"4]) 
the SM predictions provide a good description of the data over a still wide range of mt 
and m H for the values of a s and 6 a in the ranges: 0.11^o; s ^0.13 and — 0.1^5 a ^0.1. In 
conclusion, the analysis of the present precision experiments does not show a signal of 
new physics beyond the SM. 

6 Discussion 

In this section, the consequences of the update of LEP data, the new precision mea- 
surement of the left-right asymmetry at SLC |H| and the impact of a direct top mass 
measurement are considered. Finally, the predictions of all electroweak observables within 
the SM are discussed. 

6.1 Update of LEP data 



Recently the LEP Electroweak Working Group has published a report |91 



summa- 



rizing the combination of preliminary LEP data for the 1994 La Thuile and Moriond 
conferences. During 1993 the four LEP experiments have performed a high precision scan 
roughly 1.8 GeV above and below the Z resonance and within 200 MeV of m z . The new 
Z shape parameters agree with the ones quoted in section 4.1 within one standard devia- 
tion. The Z mass moved to 91.1895 ± 0.0044 GeV with improved uncertainty. Changing 
of the 'constant' m z from 91.187 GeV to 91.1895 GeV does not lead to noticeable effects 
in the analysis. The total Z width increased to 2.4969 ± 0.0038 GeV with considerably 
reduced uncertainty, also the forward-backward lepton asymmetry increased to 0.0170 ± 
0.0016. Other parameters, <r°, R^, Rb, have changed very little. The correlations in the 
Z line-shape parameter fit have become slightly smaller. 

For the time being no attempt has been made to incorporate the updated values, since 
the analyses of the 1993 data are still preliminary. 

6.2 The new left-right asymmetry data at SLC 

As emphasized in sections 3 and 4, the left-right asymmetry as well as the other asym- 
metry measurements at LEP have the advantage of determining the universal parameter 
s 2 (m 2 z ) almost independently of the other form factors, gz( m %) and 5b(m 2 z ), and almost 
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unaffected by uncertainty in a s . Since the parameter s 2 (m|) is directly related to the MS 
coupling s 2 (/x), these asymmetry measurements are particularly important for the GUT 
studies. 



The new measurement of the left-right asymmetry [31 

0.1656 ±0.0076, 



implies 



4° 



[m. 



0.2282 ±0.0010. 



(6.1) 



(6.2) 



This value is 2.5 standard deviations smaller than ( |4.5b| ) . Excluding the possibility of a 
shift caused by a systematic effect this measurement may be considered as a statistical 
fluctuation and then be combined with the other asymmetry data on the Z-pole, that is, 



the lepton (e, /i, r) forward-backward asymmetry | 26fl , the r polarization asymmetry |2 
the left-right asymmetry |!T| and the quark (b,c) forward-backward asymmetries |^6 



as 



well as with the old left-right asymmetry data from SLD [90|. The result is 

= 0.2302 ± 0.0005 . 



[m, 



(6.3) 



The new average (denoted by "ALL" ) is shown in Fig. ^5| together with the individual 
contributions.^ Note that s 2 (m|) derived from the r forward-backward asymmetry is 
as small as ( |6.2j ) from the new left-right asymmetry. Although the inclusion of the new 
left-right asymmetry lowers the s 2 (m^) fit value by about 1.5 standard deviations, the 
quality of the fit (x 2 = 6.6 for 5 degrees of freedom) does not indicate an inconsistency 
with the other data, as may be seen also from the histogram of the distribution in the 
figure. 

With the proviso of excluding a shift due to systematic error sources we include the 



data (|6.1|) into our global analysis, and discuss its effect by comparing the results with 
those obtained in section 4. The 3 parameter fit to the Z parameters only gives 



9z(™ 2 z) 



0.5538 - 0.00031 ^ 
0.2303 



-0.12 



0.01 

0.00006 a ^ 012 



±0.0017 



8b(m%) 



±0.0005 p c 
0.0071 - 0.00432 ± 0.0035 



o.oi 



0.01 



/ 1 



V 



0.11 
1 



-0.37 \ 
0.16 

1 / 



v 2 . 

A mm 



5.78 ± 



a. 



0.1000\ 2 



(6.4a) 



(6.4b) 



0.0127 

The result is shown in the (s 2 (m|), gzim 2 ?)) plane by the thick lines of Fig. |26](a) for 
three values a s = 0.11,0.12,0.13 along with the old fits (thin lines) copied from Fig. |l|. 



In Fig. ^ and in the following analysis, we use the combined result of ]90| and 



as the data for A£ R : A^ r = 0.1637 ± 0.0075 gives s s 



0.2284 ±0.0010. 
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The SM prediction for S a = is also shown in the range 100 GeV < m t < 240 GeV 
and lGeV< m H < 1000 GeV. It can be seen that the new ^4lr measurement by itself 
implies large m t (m t ^200 GeV) for m H > 50 GeV. The combined fit, however, favors 
m t ~ 180 GeV for m H ~ 100 GeV. The remaining two parameters g\{jv?z) an d ^( m l) 
are less affected. The Xmin P er degree of freedom is 8.4/6 for a s = 0.12, which is fine. 

Next, the 4-parameter fit in terms of S, T, U and 5b is performed analogously to the 
one in section 5.3. Combining the above result with eq. ( 4.35Q from the low energy neutral 
current experiments and eq. ( }4.39|) from the W mass measurements leads to 



S = -0.67 -0.024^^ +0.066^+0.30 /l 0.87 -0.25 -0.19 \ 

1 -0.42 -0.33 
1 0.19 



T= 0.30 -0.060 ^=g^ -0.004^+0.36 
U= 0.24 +0.053 +0.024^+0.54 Pcorr 

§ b = -0.0074 -0.0044 ±0.0034 



/ 



1 

(6.5a) 

ocn , f a.- 0.0998 ^ ■ ( S <* V (a. en 

X - = 8 ' 6 ° + { 0.0126 ) + lolo ) ' (6 ' 5b) 
where mt = 175 GeV and m H = 100 GeV are used to calculate SM running of the form 

factors between q 2 = and q 2 = m z . Fig. |26|(b) shows the 1-a contours in the (S, T) 

plane for the three values a s = 0.11,0.12,0.13 for 5 a = 0. The old fits ( p.8| ) are also 

shown by thin lines. The results are insensitive to the above (mt, m H ) values assumed in 

the running of the charge form factors in the region m t > 100 GeV and m H > 50 GeV, 

although they are considerably modified for m H <^50 GeV (see Fig. pj|) . It is worth noting 

this qualitative difference between the fit to s 2 (m 2 z ) and g 2 z (m 2 z ) and that to S and T. 

As a matter of fact, the experiments on the Z resonance are far more precise than those 

from the low energy neutral current experiments implying that the global fit to all the 

electroweak measurements in the neutral current sector measures essentially s 2 (m 2 z ) and 

g z (m 2 z ). In the SM the two charge form factors can be calculated for arbitrary m t and 

m H , as shown in the figure for m t =100-240 GeV and m^=l-1000 GeV. On the other 

hand, in our definition, the T parameter determines g z {0) rather than g z {m 2 z ). Hence, 

only if the running of the g 2 z {q 2 ) between q 2 = and q 2 = m\ is small, can we make 

the global fit to the S, T parameters. For this reason we restrict the SM predictions to 

the region m H =50-1000 GeV in the (S, T) figure. It is remarkable that the electroweak 

data including the new left-right asymmetry measurement clearly favor negative S, thus 

putting severe constraints on technicolor models Note that in the (S, T) plane only 

the S parameter is strongly affected by the new Alr data, while the T parameter is 

constrained, independent of the S parameter, by g z {m 2 z ) from V z . 

Next, the impact of the left-right asymmetry measurement on the SM fit is discussed 

using all electroweak data. Fig. [27] shows the results of the SM fit in the (m t , m H ) plane 
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for at a (m z ) — 0.11, 0.12, 0.13, and for S a = —0.1 (a), (b), and +0.1 (c). The contours 
of x 2 = Xmm + 1 an d X 2 = Xmin + 4.61 are shown by thick lines. The minima of \ 2 m 
the figure are marked by crosses : 12.1, 11.4, 15.7 for a s = 0.11, 0.12, 0.13, respectively, 
for 5 a = 0. The 1-a contour for each a s value is now clearly outside the physical region 
allowed by the direct Higgs searches at LEP (m H > 63 GeV, denoted by "LEP limit" 
in the figures), although the m H dependence of the \ 2 is very mild for a s >0.12. The 
result favoring a light Higgs boson reflects the fact that the new left-right asymmetry 
measurement shifts the S parameter to negative values. 

Finally, the status of the SM fit is studied in detail as in section 5. To this end the 
representation of the \ 2 °f the SM fit including the new left-right asymmetry data is 
obtained (analogous to section 5.4) : 

(TTlt — (.TTlt) \ ^ 
— ^ J + Xh(%i a s, 6 a ) , (6-6a) 

where 

, rrifT , 9 itltt /a, — 0.12\ / 5„ \ 

(m t ) = 162.2 + 12.6 In — ^ + 0.8 n 2 — ^ - 1.2 — - 4.8 —2- , 6.6b 

v 1 ioo ioo v o.oi J vo.ioy ' v ; 

, win / TriTT\ m t — 175 

Am t = 12.0 - 0.09 In — £ - (0.31 - 0.05 In — ^ ) — , 6.6c 

100 V 100^ 10 v ; 

and 

2, rx nKR , f5 a -0.75\ 2 fa s - 0.1164 + 0.005 5, 

XH{ m H, a s,6 a ) = 9.56+ — — — — + 



0.39 J \ 0.0060 
a s - 0.1365 + 0.030 5 a \ m H fa s - 0.1255 \ 2 m H ^ 

0.0144 ) n 100 ~ V 0.0639 J ^ 100 + vOlO 

(6.6d) 

Fig. |28] (in analogy to the previous results of Fig. |23|) shows the total \ 2 °f the SM 



fit as functions of m t for m H =60, 300, 1000 GeV and a s (m z ) = 0.11, 0.12, 0.13. The 
uncertainty 5 a is shown for three cases : 6 a = —0.1 (a), (b), +0.1 (c). The dotted lines 
are obtained by the approximate formulae ( |6.6|) . It is obvious from Fig. |28] and Fig. or 



from eq. (|6.6b| ) and eq. (|5.9b|) , that the best-fit value of m t is shifted by about +17 GeV 



for given m H , a s and 6 a values. Here again the uncertainty of 5 a is important for the top 
mass prediction, as observed from ( |6.6b|) and Fig. |28] : 5 a = ±0.1 causes a shift +5 GeV in 



the best-fit value (m t ). The a s -dependence of the (m t ) values is considerably weakened. 

Fig. |29] (in correspondence to Fig. |24| in the previous fit) shows the total x 2 of the SM 
fit including the new left-right asymmetry data [[3TJ, as functions of m H for mt =100- 
200 GeV. Three a s cases are displayed; a s {m z ) = 0.11 (a), 0.12 (b), and 0.13 (c), all 
for 5 a = 0. The dotted lines show our approximation (|6.6|) , valid only in the 'physical' 
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region 63 GeV < m H < 1000 GeV. As seen, the best-fit value of m H is as low as 10 GeV 
for m t <^150 GeV, while it increases with m t for m t ;>150 GeV. This trend can also be 
appreciated from the global fit of Fig. p6j(a) in the (s 2 (m|), g%{m 2 z )) plane. 



6.3 The impact of the top mass measurement 

The top quark searches of the two collaborations CDF and DO at the Tevatron entered 



in their decisive phase ||105| , |106|1 . The range of values for the top quark mass coming out 
of the fits to the electroweak precision data is within reach for direct observation in the 



detectors at the Tevatron. In view of the recent publication by the CDF collaboration [ |106j| 
it is instructive to examine the impact of the constraint 



m t 



17 A ± 16 GeV. 



(6.7) 



First, the m t - dependence of the global fit to the electroweak data in terms of the 
charge form factors s 2 (m|) and gf(?7i|) is considered, now assuming SM dominance to 
the 5fe(m|) form factor. Using the Z parameters including the new Alr measurement [31] 



one obtains 



0.55430 - 0.00109 - 0,12 a0023mt i6 174 



s 2 (m 2 
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0.19, (6.8a) 



(6.8b) 



V 60 J ' V 0.0068 
which is a good approximation in the region 150 GeV < m t < 200 GeV. Here the errors 
and the correlations are almost independent of the mt value. The fit to all electroweak 
data gives 

S =-0.62-0.097 
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63 J \ 0.0067 
The appearance of essentially the same combination 

m t - 174 
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0.10 



a, - 0.12 - 0.0022- 



16 



(6.9a) 
(6.9b) 

(6.10) 



in eqs. fl6.8|) and (|6.9|) is the expected consequence of the strong correlation between 
5b{m 2 z ) and a s as discussed in detail in section 4. 
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Next, the above constraint on the top quark mass (|6.7|) is imposed on the x function 
of the SM fit in the previous subsection. The result displayed in Fig. 31 shows the 



improvement over Fig. 27. Now, light Higgs boson masses are moderately favored, as a 



consequence of the constraint (6/7) being somewhat larger than the best-fit value of m t 
obtained by freely fitting the two parameters, m t and m H without the m H constraint from 
LEP. 

It is instructive to anticipate the impact a precise measurement of the top mass would 
have in the context of the present electroweak data. The top quark mass is expected to be 
measured eventually with an uncertainty of about 5 GeV at Tevatron by the end of this 



decade [107|, which may be improved to about 3 GeV at an upgraded Tevatron ||108|| . The 
uncertainty is expected to be reduced by an order of magnitude to a few hundred MeV 
at next linear e + e~ colliders | 109 |. The top mass acts then like an external parameter 
and the only remaining free parameter is the Higgs mass. Fig. shows the 95% CL 
constraints for three values of a s {m z ) = 0.11,0.12,0.13, and for 5 a = 0. For small mt 
values, rather strict upper bounds on m H are found. On the other hand no strict upper 
bound is obtained for m t ^180 GeV. In the region 160 GeV < m t < 190 GeV, the upper 
bound on m H at the 95% CL is approximately expressed as 

1.20 + 1.12^^ for a s = 0.11 
ln^<<| 1.55 + 1.25^P for a s = 0.12 , (6.11) 



100 



16 

1.95 + 1.45^=^ for a s = 0.13 



16 

where m t and m H are measured in GeV. The upper bound is lower for smaller m t . Since 
these bounds are very sensitive to the m t value as well as the assumed a s value, further 
accurate measurements of m t are needed to obtain more stringent limits on m H . Never- 
theless, it is remarkable that the constraint on the top quark mass ( |6.7|) would favor a 
relatively light Higgs boson, m H = O(100 GeV), which may exist in the minimal SUSY- 
SM. Also, the direct m H bound from LEP [JlOTJ ™h > 63 GeV (95%CL) implies that the 
top quark should be heavier than about 145 GeV, since otherwise the Higgs boson should 
have been discovered at the 95%CL. This lower m t bound changes by about =f5 GeV for 
5 a = ±0.1. 

One comment is in order. Though our approximate formulae of the \ 2 f° r the SM 
fit, (|6T6|), reproduce the exact result within about 1% accuracy in the Higgs mass range 
63 GeV < m H < 1000 GeV as seen Figs. [ZB|-[Z9], one should not use them in finding the 



confidence levels of m H for small m t , since the neighborhood of the minimum of the x 2 
is outside the above range, where the exact \ 2 an d the approximate formulae are fairly 
different as seen from Fig. BO. 
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6.4 Summary of the data and the SM fit 

Table |^ collects the complete list of all input data (except for a, Gf and m z ) and the 
corresponding minimal SM predictions for several sets of (m t ,m H ,a s ) values. The total 
X 2 of each sector is also given in the table. The correlations between the errors (given 
in the text) are properly taken into account. The numbers demonstrate that the present 
electroweak experiments are well described by the SM, perhaps except for a combination of 
a light top and a heavy Higgs, see the case (mt, m H ) = (150, 1000)GeV in the last column 
of the table. Its total x 2 a t oc s = 0.12 is 30.22 for 19 data points, whose x 2_ P r °bability 
corresponds to 95%. In Table [F] also the results of two approximations are listed. The 
'no-EW column is obtained by dropping all electroweak corrections to the two-point 
functions (S = T = U = 0) as well as vertex/box corrections (6 G = 5b = r, = By = 0), 
while retaining the QED running of the charge form factors a(q 2 ) and s 2 (q 2 )/a(q 2 ) due to 
light particles (excluding the W and t contributions). The 'IB A' column shows the result 
of the improved Born approximation, where all the gauge boson propagator corrections 
are retained and hence all the four charge form factors are kept exact, but all vertex/box 
corrections (5b — A = B^ = 0) dropped, except for 5 G in the // decay. 

It is amazing to note that the 'no-EW hypothesis is, from a statistical point of view, 
not completely unacceptable. The comparison between the 'no-EW and the 'IBA' hy- 
pothesis is surprising, since in the 'IBA' prediction all the most important electroweak 
corrections are supposed to be contained, including the dominant m 2 corrections in the T 
parameter. It is even more striking, if 5 G in IBA is set to (this may be called a genuine 
IBA), to obtain s 2 (m 2 z ) = 0.2286 for m t = 175 GeV and the total x 2 jumps nearly to 100. 
The measurement of the Z parameters are equally well described by the 'no-EW and the 
full calculation for m t = 175 GeV. This confirms the observation of ref. |6~3] , |110| | that there 
is no evidence of the genuine electroweak correction in the present electroweak precision 
experiments. As explained in sections 2.3 and 5.3, this is because of the accidental can- 
cellation between the propagator corrections and the remaining vertex/box corrections. 
The no-EW calculation for all the asymmetries on the Z-pole give almost the same values 
with the predictions of the exact calculation for m t = 175 GeV and m H = 100 GeV. As 
discussed in section 4.1, R b also gives a large contribution to x 2 m the full calculation. 
For a large top quark mass, the Zb L b L vertex from factor 5b(m 2 z ) decreases (see Fig. [I]), 
and hence it gives smaller Rb. For this reason the present data of Rb agree better with 
the no-EW and the IBA calculations, where 5b(m 2 z ) is set to 0. 

The most significant differences between the no-EW prediction and the full SM pre- 
dictions in (m t = 175 GeV , m H = 100 GeV) column appear actually in the predictions 
for the low energy u^-q scattering and the atomic parity violation experiments. When 
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evaluating the no-EW and IBA predictions, all the external photonic corrections and the 
tree-level propagator effects are retained, as explained in section 3.2. The difference be- 
tween the full SM predictions and the no-EW or IBA predictions is mainly caused by the 
absence of the WW box contribution in the latter. 

Another significant difference appears in the predictions for m w , where the no-EW 
prediction (79.95 GeV) is much smaller than the observed value, 80.24 ± 0.16 GeV. This 



observation has also been made in refs. ||110H112|] . In contrast to the low energy neutral 
current experiments above, the difference here is due to S and U contributing to m w 
proportional to -0.2945 + 0.332C/ (c.f. eq. (|3~T00|) ). For instance, the full SM for m t = 
175 GeV and m H = 100 GeV predicts S = —0.2323 and U = 0.3577, which implies for 
m w a shift by 0.19 GeV corresponding to more than one standard deviation. 

Finally, Fig. shows separately for each sector the \ 2 of the SM fit as functions of 
m H for m t = 100 — 200 GeV. In all sectors, the preferred Higgs mass range is strongly cor- 
related with the assumed top mass. For m t =170-180 GeV, a light Higgs boson is favored 
by the Z parameter measurements and by the low energy neutral current experiments, 
while the data of m w alone prefer a rather heavy Higgs boson. Although the overall trend 
of the total x 2 shown in Fig. |29| is dominated by the contribution from the Z parameter 
measurements, also the W mass measurement plays an important role for some m t , m H 
ranges. For instance, a relatively light Higgs boson (m H <^100 GeV) appears incompatible 
with a heavy top quark (m t ~ 200 GeV) by the m w measurement alone. 



7 Conclusions 

A novel method to confront electroweak data with theory at the quantum level has 
been proposed and a comprehensive analysis has been carried out. The electroweak ob- 
servables were first expressed in terms of model-independent parameters, which in turn 
were expressed in terms of S'-matrix elements of processes with four light fermions and 
factorized into the short-distance part and the part related with the external QED/QCD 
corrections for neutral current processes. Only two quantities, the Fermi coupling constant 
Gp and the W mass are considered for charged current processes. Since all electroweak 
observables were expressed in terms of helicity amplitudes, they can be evaluated in an 
arbitrary model on and off the Z resonance. Our formalism is hence useful to study effects 
of tree-level deviations from the SM, arising, for instance, from an additional Z boson. 
After careful evaluation of the external QED / QCD corrections, the theoretical predictions 
were confronted with experiment in three steps of increasing theoretical stringency. First, 
in the class of theories respecting the electroweak gauge group SU(2)l X U(l)y broken 
spontaneously to U(1)em the radiative effects were classified into process-independent and 
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process-dependent ones. Apart from the Zb L b L vertex, all vertex and box corrections were 
assumed to be given by the SM, while new physics contributions were studied in the most 
general way by four universal charge form factors. Next, by assuming the running of the 
charge form factors to be governed by SM physics alone, the electroweak parameters S, 
T, U were determined. Finally, the SM itself was confronted with experiments. 

It was our aim to render this analysis as transparent as possible by developing the theo- 
retical formalism in full detail and by presenting the results in figures and parametrizations 
in a form useful for appreciating consequences of future improvements in the experimental 
data. 

The analysis proceeded in two steps. First, the information in the whole body of 
electroweak precision data has been condensed in the 9 electroweak parameters : m w and 
m z , e 2 (0), s 2 (0), g|(0) and ^(0), s 2 (m z ), g z {m z ) and 5 2 (m 2 z ). At the present time no 
direct information exists for e 2 (m|). In order to keep the analysis flexible e 2 (m|) and 
also the QCD coupling constant a s have been treated as external parameters in the fit 
procedure. Second, this universal set of quantities with the complete covariance matrix 
has been interpreted within the electroweak theory at three qualitatively distinct levels. 

The main result is that the data can be consistently interpreted at all levels, in par- 
ticular there is nowhere evidence against the SM. This conclusion is not affected, when 
the new precision measurements of the left-right asymmetry from SLD |3]J is included. 
The fits to the universal charge form factors or that to the universal S, T, U parameters 
work well and do not hint at a violation of the SU(2)l X U(1)y universality, nor at an 
anomalously large non-standard vertex/box corrections. Generally speaking, the inclu- 
sion of the SM vertex/box corrections improves the fit to the data, while the improved 
Born approximation gives a poor fit to experiments. The ratio Rb = 0"°/cr° measured by 
the LEP experiments turned out to be in poor agreement with the large Zb L b L vertex 
correction predicted by the SM. The fit to the S, T, U parameters gives us information on 
spontaneous symmetry breaking. The T parameter is essentially determined by the charge 
form factor g§(^|), and positive value is favored. The 5* parameter is then fixed mainly 
via s 2 (m 2 z ), and hence its best-fit value is affected by the asymmetry data. A negative 
S value is favored by the new left-right asymmetry from SLD, and the naive technicolor 
models are disfavored f|. Due to strong correlation between the fitted S and T values, the 
region of the (S, T) plane with relatively large S and T (— 0.3^S'^ — 0.1 and 0.5<^T^1) is 
consistent with the SM prediction for 150 GeV^m^200 GeV and 50 GeV^m H ^200 GeV 
(see Figs. ^6|(b) and |30](b)). The U parameter is measured only via g^(0), and it is 
consistent with zero. 

The analysis showed that the experimental precision required to detect a deviation 
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from the SM is still insufficient. For instance, the running of the charge form factors 
can be tested presently only for s 2 (q 2 ) and gf (g 2 ) and is limited by the precision of the 
low energy neutral current experiments. Nevertheless, the data are precise enough to 
show that their consistent description within the SM is only guarantied, if the top quark 
mass exceeds about 145 GeV. This low mass bound of m t is nearly independent of a s , but 
changes by about =f5 GeV due to the uncertainty ±0.1 in the hadronic vacuum polarization 
contribution to S a = l/a(m|) — 128.72. Note that the SM top-bottom contribution to 
the form factors have been calculated by using the 0(aa s ) two- loop formula [Q-Q. 
Perturbative ti threshold effects |60-|62| will affect these corrections, and the predicted 
m t value may shift upwards by as much as a few GeV [62]. 

The near future promises a clarification of the value of the left-right asymmetry pub- 
lished by the SLD group and the ratio Rb from LEP experiments. The precision scan 
around the Z resonance performed 1993 by the four LEP experiments will further im- 
prove substantially the Z resonance parameters. It would be advantageous to publish the 
data without the subtraction of the Z-7 interference contribution. Eagerly awaited is the 
definitive observation of the top quark. If its mass turns out to be compatible with the 
electroweak analysis of the 1-loop effects there is hope to constrain for the first time the 
elusive Higgs sector. 

By introducing the QCD coupling strength Q s ( m z)MS an< ^ ^ ne sn ^ ^ a = — 128.72 
as external parameters in the fit, we have made clear the significance of their precise 
measurements. Unless these parameters are accurately measured, the search for effects 
beyond the SM through the electroweak radiative effects gets increasingly limited. 
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App. A SM radiative corrections at one-loop order 

In this appendix the propagator, vertex and box corrections of the standard model 
(SM) are presented, all at one-loop level and partly at two- loop level for the 0(aa s ) terms 
of the gauge boson propagators. All the Green's functions are calculated in the 't Hooft 
Feynman gauge in the dimensional regularization and renormalized in the MS scheme. 
Definitions of the scalar one-loop integrals, A, B, C, D functions, are given in Appendix 
D. Vector boson propagators are given in A.l, the vector boson fermion vertex functions 
and the fermion wave function corrections follow in A. 2, while the box corrections are 
listed in A. 3. All the one-loop calculations are done independently and we reproduce the 
known results of ref. |3"2] , |TT3| for two-point functions and those of ref. [4*T-fP|j for the 
three- and four-point functions. 

A.l Propagator corrections 

There are four vector boson two-point functions contributing to processes with external 
light quarks and leptons at one-loop order. They can be parametrized by [Q] 

W?{q 2 ) =e 2 Il?V), (A.la) 
(g 2 ) = eg z {H?V) - s 2 Tf T Q ( q 2 )} , (A.lb) 

nf (g 2 ) = gl {Tf T \ q 2 ) - 2S 2 n?V) + * 4 n? V)} , (A.ic) 

nrV)=S 2 lf (g 2 ), (A.ld) 



with the coupling factors 



9z = ~ = ^, (A.2) 



c sc 

and the use of the compact notation 

s 2 = 1 - c 2 = sin 2 6 W , (A.3) 

throughout the appendix. These two-point functions and the coupling factors are renor- 
malized in the MS (the modified minimal subtraction) scheme, and hence they depend 
on the 't Hooft unit of mass fi which appears explicitly in the B functions as defined in 
appendix D. The coupling factors of ( |A.2| ) and ( |A.3| ) also depend implicitly on the unit 
of mass /i. The subscripts T in eqs. ( |A.1| ) denote the transverse part of the polarization 
tensor 

IV(g 2 ) = {-g, v + ^)lW) + ^n L (g 2 ). (A.4) 
The longitudinal parts HiiQ 2 ) do not contribute to processes with light external fermions. 
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With the help of the four B functions, Bq^B^^B^ and B 5 (see appendix D), all SM 
contributions to the above two-point functions are expressed compactly. Il r (g 2 )'s is 
decomposed into the bosonic and the fermionic contributions, 

Tlf(q 2 )=Tlf(q 2 ) B + U^(q 2 ) F , (A.5) 

and the expressions are given separately. 

A. 1.1 Bosonic contributions 



The bosonic contributions with pinch terms are given by |34 

1 

Air 2 



Ii Q T Q {q 2 ) B = Y^ Q {q 2 ) B ~-^q 2 B (q 2 ; W, W) 



1 



(A.6a) 
(A.6b) 
(A.6c) 



Tf T \q 2 ) B = Uf(q 2 ) B ~-^(q 2 - m 2 w )B (q 2 ; W, W) 

4 7T Z 

U 1 T \q 2 ) B = U 1 T \q 2 ) B -^- 2 (q 2 -ml v ) [c 2 !^ 2 ; W, Z) + s 2 B (q 2 ; W, T )j , (A.6d) 

4 71 J 

where the short-hand notation 

B n (q 2 ; A, B) = B n (q 2 ; m A , m B ), (A.7) 

is introduced for the B functions. Each Ht(<7 2 ) function without overline is calculated in 
the 't Hooft-Feynman gauge, whereas the n T (g 2 ) functions with pinch terms are gauge 
invariant ||34[| . The explicit expressions are 



n T (q ) B 



n 



•3Q/ 2^ 



16 7T 2 

-,2 



q 



n 



33, 2 ^ 



)B 



16 7T 2 
1 

16 7T 2 
1 

16 7T 2 
1 

16 7T 2 
1 

16 7T 2 

fi2 



5S + 12B 3 
^B + 10B 3 



(q 2 ; W,W) + ^ 



(q 2 ] W,W) + 



m 2 z B + ^B 5 



1 

4" 



(q 2 ;Z,H) 



23 
T 



g 2 - 2m 2 w ^B + 9q 2 B 3 
W,H) 



m 2 w B$ + - A B^ 



1 

4" 



(g 2 ;^,W) 



9 



24 7T 2 



JcY - (1 - 4s 2 )m 2 y - m 2 z )B - ( ^ - 2s 2 ]£., 



■ 7T Z 



(4g 2 - 2m 2 w )B - £? 5 



(g 2 ;^, 7 )- 



24 7T 2 



(A.8a) 
(A.8b) 

(A.8c) 

(q 2 ;W,Z) 

(A.8d) 



At one-loop order of the minimal SM, the first terms in eqs. ( |A.8c[ ) and ( |A.8d| ) are the 
only ones in the transverse component of the vector boson propagators being dependent 
on the Higgs boson mass (m H ). 
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A. 1.2 Fermionic contributions 



The fermionic contributions to the gauge boson propagators are known to 0(aa s 
two-loop level : 

,2 



n?V)F 



i 

T C g J2 Q 2 f {8B 3 tf-JJ) + ^B> v (q 2 -JJ)}, 

f=Ui,di ° 



16 7T 2 
,2 



16 7T 2 



(A.9a) 



i 

T C q ]T Q f lJ^B 3 ( q 2 ;f,f) + l^B' v ( q 2 ;f,f) 

f=Ui,di 



16 7T 2 

f_ 

16 7T 2 



(A.9b) 



16 7T 2 



£ (/ 3 /) 2 [VS 3 -2mJi?ol(g 2 ;/,/) 



16 7T 



1 <?* E (/3/) 2 i[4g 2 53-2mj5 ](g 2 ;/,/)+ 1 ° 



f=lH,0 



3 7T 



5y + 5, 



(A.9c) 



^V)^ = TTT^ Ef 2 ^ - B^q 2 ; v^) 

16 7T Z ~^ L 

+J^CflE|^ 'I [2g 2 B 3 - B<\ (q 2 ; u h d 3 ) 



16 7T 2 



'■J 



5y + B, 



(g 2 ; dj-] 



(A.9d) 



The summation over i, j extends over the three generations of lepton and quark flavors, 
0i,^2,^s) = (uejj/^Ur), (£i,£ 2 ,h) = (e,/i,r), (ui,u 2 ,u 3 ) = (u,c,t) and (d 1 ,d 2 ,d 3 ) = 
(d, s,b). C q = 3 is the color factor, Qf the electric charge of the fermion / in units of the 
proton charge, I 3 f the weak isospin 



'3/ 



+| for / = Vi or u iL , 
-\ for / = i iL or d iL , 
for / = £ iR , u iR or d iR , 



(A.10) 



while V^. are the Kobayashi-Maskawa quark mixing matrix elements. The 0(aa s ) cor- 
rections in perturbative QCD |52, 55|,|5q| are given by the functions By and B A '■ 



B v (q 2 ;m,m)=q 2 B' v (q 2 ;m,m), 

B A (q 2 ] m, m) = q 2 B' A (q 2 ; m, m) + B^(0; m, m) 
B v {q 2 - m, 0) = q 2 B' v (q 2 ; m, 0) + B v (0; m, 0) , 
B A {q 2 - m, 0) = q 2 B' A (q 2 ; m, 0) + B A (0; m, 0) , 



(A.lla) 
(A.llb) 
(A.llc) 
(A.lld) 
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where 



and 



^(g 2 ;m,m) =ln^ + — -4C 3 + , 



4m 2 ' 

/ / 9 \ , U 2 55 , 4m 2 r . . q 2 
B' A {q 2 ;m,m) = hi At + — _ 4 Ca + — [^i ' 



m 2 12 
^(g 2 ;m, 0) = ^(g 2 ; m, 0) 
■ 2 55 



m 2 12 



4m 



4m 2 



- M 



i?4(0; m, m) 



m 



12 ln 2 4 + 22 ln4 + - 



B v (0;m,0) = B A (0;m,0) 



m 



n , 2 fi 2 11 . /i 2 u 35 
3 In 2 ^- + — In ^— + C 2 + — 



m 



(A.12a) 
(A.12b) 

(A.12c) 

(A.13a) 
(A. 13b) 



Here ( 2 = 7] " 2 /6, C3 = 1-2020569, and the complex functions Vi, A 1 and Fx are given in 
refs. i55|)p6|. The following limits are useful : 



^(M 2 ;0,0) = ^(M 2 ;0,0) 

1 /W 2 55 

= In -^r H 4( 3 + ^7^ . 

M 2 12 S3 

11 2 1 Fi M 2 
B' v (M 2 ;m,m)=ln^ + - + 0(—) } 

^(M 2 ;m,m)=ln^ + | + 0(^), 



36 

B' V (M 2 ; m, 0) = ^(M 2 ; m, 0) 

-ln^l + — --C +Of— ) 
m 2 36 9 ; 



(A.14a) 
(A. 14b) 
(A.14c) 

(A.14d) 



A. 2 Vertex correction 

The vertex form factors r{(g 2 ), T^iq 2 ), r{(g 2 ) and r{(g 2 ) appearing in the helicity 
amplitudes (|2.2| ) contribute to the 7// and Zff vertices as follows : 

r -r// (g 2 ) = _^ Qf[1 + T f (q 2 )] + j 3/ r^( g 2 )}, (A.15a) 

T Z// (g 2 ) = -$ Z {(I 8 , - Q/S 2 ) [ 1 + T{(q 2 ) } + I 3/ [ C 2 V[ (q 2 ) + r{(g 2 ) ] + r{(g 2 )} . (A.15b) 

It should be noted that the functions T{(q 2 ) and T 2 (q 2 ) are common to the 7// and 
Zff vertices, and that 1^3 (g 2 ) and T{{q 2 ) are additional contributions to the Zff vertex. 
These vertex functions depend on the chirality of / and their explicit forms at one-loop 
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level of the SM are 



r?(g 2 ) = r 3 ^( g 2 ) = rf(^) = o 



r fV) = (^r) r{ z (g 2 ) + E 



47T 



riV(g 2 ) 



rf(g 2 ) 



-2E 



Air 



V 



rV(g 2 ) 



r 3 /L (g 2 ) 


= E 




/' 


r{ L (g 2 ) 


= o, 



47T 



(A.16a) 
(A. 16b) 

(A.16c) 

(A.16d) 

(A.16e) 
(A.16f) 



with the gauge boson coupling convention 



9l j = 9r =eQ f , 



-.iff _ „jff 
9l" ' = 9z(hf-Q,'s 2 ) 



A" 



Wff 

9l 



~9z Qf 



9 

V2 



V ff ' 



(A.17) 



where 



T{ z {q 2 ) = T^q 2 ; f, Z, f) - E'(m 2 ; f,Z) , 

?(w{q 2 ) = (Tx + r lm ) (g 2 ; /', W, f) - £'(mj; /', W) 



(A.18a) 
(A.18b) 



r^(g 2 ) = (r\ + r lm ) ( q 2 -, f, w 7 /') - r 2 ( g 2 ; w, f, w) + 2Re B ( q 2 ; w, w) , (A.isc) 

Ti w {q 2 ) = r lm (g 2 ; /', W, /') + T 2m (q 2 - W, f, W) . (A.18d) 



Here 



S'(g 2 ; m,M) = -(2 + ^)Bi(q 2 ; m, M) — 1 



(A.19) 



is the external light fermion self energy correction, and the last term in eqs. (|A.18cj) of 
T 2W is the pinch term @,|34| which is subtracted from the vertex functions as calculated 
in the 't Hooft-Feynman gauge. The remaining vertex functions in eqs. ( |A.18j ) are 



ri(q' 2 ; m, M, m) 
T lm (q 2 ;m,M,m) 



m 



2q 2 (C n + C 23 )+4C 24 -—C 



(g 2 ; m, M, m) 



ffP 1 ( r -i I 

W2 { [q 2 (C 12 + C 23 ) + 2C 24 - 2M 2 C ] (g 2 ; m, M,m) — - 



(A.20a) 
, (A.20b) 



r 2 (g 2 ; M, m, M) =2 [g 2 (C n + C 23 ) + (e + ^) C 24 + (g 2 - m 2 )C 



m 



T 2m (q 2 ;M,m,M) = 2— 



2M 2 C - C 24 



(q 2 ;M,m,M) 



(q 2 ;M,m,M) - 2, 
(A.20c) 

(A.20d) 
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with the shorthand notation for the C functions of appendix D: 

Ci(q 2 ; mi, m 2 , m 3 ) = Cj(0, 0, q 2 ; mi, m 2 , m 3 ) 



(A.21) 



In the limit of the diagonal KM matrix elements V,^ = Sij, which is assumed in 
all our numerical results, the internal fermion mass m — mp is non-negligible only for 
/ = = 0- Otherwise we can set m = at high energies (m 2 /q 2 ~ 0) and find 



/i 2 A M 2 , o M 2 x -q 2 -it 



r l( g 2 ;0,M,0) = ln^-4-2— + (3 + 2— )ln. 



M 2 q 

M 2 -^2 



M 2 



(A.22a) 



r / L - V m 2 

T 2 (q 2 ; M, 0, M) = 3 In £ + 2 - 2^ + (l + 2^)/3L + 2^( 2 + ^)l 2 ,(A.22b) 



9 



T lm (g 2 ; 0, M, 0) = T 2m (q 2 ; M, 0, M) = , 



and 

S'(0;0 > M)=ln^-i, 

R eJ B (g 2 ; M, M) = In + 2 — f3L . 

Here Sp(z) = — J 2 ^ is the complex Spence (dilogarithm) function, and 

13 = yjl - 4(M 2 - ie)/q 2 , 
/3+1 



L = In 



(3-1 



(A.22c) 

(A.23) 
(A.24) 

(A.25a) 
(A.25b) 



At low energies, light fermion masses may not be neglected as compared to the mo- 
mentum transfer q 2 . In the limit of |g 2 |/m|<^l and m 2 /m|-<l, but at fixed m 2 /q 2 , the 



vertex functions reduce to 





-q 2 

m| 


M q 2 ;f) + 0(^)], 


rUf) = 


-q 2 

m 2 w 


Mq 2 ;f') + o(^)' 


r f 2W (q 2 ) = 


-q 2 

m 2 w 


Jw(q 2 ;f') + o(^)' 



where 



J v (q 2 -J)=4F 3 (q 2 ;fJ)-^\nm 2 v -±, 
J w (q 2 ; f) = 4F 3 (q 2 ; /, /) - \ \nm 2 w - \ - \ 



(A.26a) 
(A.26b) 
(A.26c) 

(A.27a) 
(A.27b) 



The function F 3 is defined in appendix D. The last 1/3 factor is the pinch term. 
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A. 3 Box correction 



Box corrections for the process e\e\ — > f a f a are expressed by B^[, where A, a = — 1 is 
used for left-handed fermions and A, a — +1 for right-handed fermions. 



+ 



i 



167T 2 
5\L 



16tt 2 
1 



16tt 2 



n WeuWff 

9l 9l 



9l ee 9 Z -{ f 



[Ii(u,s;m z ,m f ) - I 2 (t, s;m z ,m f )] (A. 28a) 

2 ^ f +I 1 (u,s;m w ,m f ,) for hf = +\ (/ = v t ,Ui) 
-I 2 (t,s;m w ,mf>) for I 3f = -\ (f = £,di) 



x 



[I 2 (u, s; m z , m f ) - h(t, s; m z ,m f )\ , 



with 



s = ( Pe+Ps ) 2 = (p f + p f -) 2 , 

t = (Pe-Pf) 2 = (Pe-Pf) 2 , 
U = (Pe-Pf) 2 = (Pe-Pf) 2 , 



(A.28b) 



(A.29a) 
(A.29b) 
(A.29c) 



and pi being the 4- momentum of particle i . The internal fermion mass rrifi is non- 
negligible only for f — b L , for which the top quark contributes in the limit of the diagonal 
KM matrix elements. The functions h(u, s ; m v ,mf) and I 2 (u, s ; m v ,mf) are expressed 
in terms of the D functions of appendix D: 

h(u, s; m v ,m f ) = -2u (D n + D 12 - D 13 + 2D 24 ) - 4tD 25 - 4sD 26 - 16D 27 , (A.30a) 
I 2 (u, s; m v ,m f ) = -2u (D n + D 24 - D 25 ) - AD 27 , (A.30b) 

where 

Di = A(0, 0, 0, 0, u, s; 0, m v ,m f , m v ) , % = 0, 11-13, 21-27 . (A.31) 
After reduction of the higher D functions 



Ii(u, s; M, m) 
I 2 (u, s; M,m) 



-2Cf ) -2Cf 4) + 2(«- m 2 ) 4 



(1234) 



(A.32a) 



s + u 



(13) 



B 



(24) 



+ 



U 



(s + 2 u + 2 M 2 - m 2 ) [cf 23) + C { 



.(134) 



S + u) 2 

^— - ( s ( s + 2 u - 2 M 2 + m 2 ) + 2m 1 



{s + u 
1 



^ 124) + cj- 



(234) 



u — m 



..... ,,)s 2 -(2M 4 -2M 2 m 2 +m 4 + 
(s + -u) 2 I v 7 v 

+2 (M 2 + u ) (M 2 -m 2 + u) u\d£ 2U) , 



m 2 u 



2 m 2 ) s 
(A.32b) 
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is obtained. For the case / 7^ b L the limit m — > can be carried out 
J 1 ( M , S ;M,0) = -4C( 124) + 2n^ 1234) , 



(A.33a) 



I 2 {u,s;M,0) 



s + u 



d(13) r (24) 



(s + u) 2 
2 

(s + w) 2 



(123) 



(s + 2 u + 2 M 2 ) c£ 



s(s + 2!i-2M 2 ) + 2u 2 )d 124) 



— ^--| S 2 u - 2 s (ilf 4 - u 2 ) + 2 u (M 2 + n) 2 }^ 1234 ) . (A.33b) 



Eqs. ( |A.33| ) agree exactly with ref. j|l 



Am v ' 



In the low energy limit, only .D27 survives: 
D 27 (0,0,0,0,0,0;0,m y ,0,m y ) - ' 
^27(0, 0, 0, 0, 0, 0; 0, m w , m t , m w ) 
and hence 



4(m 2 - m%) 



m. 



m, 

9- In — ~ 1 

m t - ™w m w 



+ 



5a £/ -a(0,0) 
for / 7^ bL, and 

5 L e ,i(o,o) = 

for / = b L . 



16tt 2 

16vr 2 
1 



2 3 



TO; 



n WeuWff> 



for J; 



x 



m 



3/ 



JV 



for h f = -\ (f = £,di) 



W 



16rr 2 



2 3 



2 • 



r .Zee r .Zbb 

9l 9l 



n Wev n Wbt 
iJL iJL 



16rr 2 m| 



167T 2 m 2 — m 



IV 



1 m t 1 
In —| 1 



m 2 - m 2 ^ m 2 ^ 



(A.34a) 
(A.34b) 



(A.35a) 
(A.35b) 



(A.36) 
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App. B Renormalization group and hadronic contributions 

The effective charges of the SU(2) x U(l) theory are expressed in terms of the MS 
couplings by 



g 2 (q 2 ) g 2 (fi) 
where the 577(2) effective charge 



1 1 +RenS(g 2 ), 



(B.l) 
(B.2) 



s 2 (q 2 



(B-3) 



is introduced for convenience. The expressions ( B.l|) and (|B.2| ) are explicit solutions of 
the renormalization group (RG) equation in the MS scheme : 



D [effective charges] = 0. 



(B.4) 



with the RG operator 



D 



fi 



d/j, 2 
_d_ 



+&(-. 



d 



16tt^ <9(e 2 /16vr 2 ) M9 ^16n 2 ^ d(g 2 /167v 2 ) ' 







The MS /3-functions in the minimal SM read at one-loop order 

4. 



-7 



-HfCfQ 



16tt 2 



6 



(B.5) 



(B.6) 
(B.7) 



where Cf = 1(3) for / = £(g). The two-loop 0(aa s ) contributions are accounted for 
by replacing C q — > C q {l + 24 ) in eqs. ( |B.6|) and ( P-7|) . Note that the effective charges 
e 2 (q 2 ) and g 2 {q 2 ) behaves similarly to the MS couplings at asymptotically high energies, 
\q 2 \ m 2 ^, since the functions IT T (g 2 ) and (q 2 ) do not have large logarithms at 
fi 2 ~ |g 2 | rnly 0. This is enabled by adding the pinch terms @,|34| in the self energy 
Il(g 2 ), and our e 2 (q 2 ) and s 2 (q 2 ) are equivalent to the corresponding ^-charges up to 
the imaginary parts and the two-loop corrections. 

Although the MS couplings e and g could be adopted directly in our analysis, we prefer 
the effective charges of (|B.1|) and (|B.2|) as quantities to be used when confronting theory 
with experiment. We give two reasons, one being associated with the non-decoupling of 
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heavy particles in the MS scheme, the other being related with the treatment of non- 
perturbative hadronic contributions to the electroweak parameters. 

Traditionally, the appearance of large logarithms of heavy particle masses (non- 
decoupling) in the MS scheme is avoided by adopting the effective field theory | 114| , |115 l, 



where the heavy particle fields are integrated out in the action. The couplings of the 
effective theories are then related to each others by matching conditions ensuring that 
all effective theories give identical results at zero momentum transfer, since the effects of 
heavy particles in the effective light field theory must be proportional to <? 2 / m heavy 

In general, the two MS couplings e 2 (/x) e // and g 2 (/x) e // of the effective light particle 
theory can be obtained by the matching conditions 

where only the light particles at the scale n contribute to the two-point functions at the 
right-hand side. In the minimal SM, one may, for instance, employ an effective theory of 
particles of mass up to the scale [i: 

m=i^r f + ln ^ - - ( 7 ln i + 1^ - m ^ (R10) 

16tt 2 2 n i v /43 2, 

Such a scheme is often adopted in quantum chromodynamics (QCD), but leads to a dis- 
continuity at = m w of the effective MS coupling constants. The appearance of the 
discontinuity in the unphysical MS couplings is not really a problemP], but the appearance 
of many quark and lepton mass scales renders the use of these effective couplings imprac- 
tical at the scale /i < m z . Furthermore, direct use of the effective MS couplings at lower 
energies leads to expressions with light-quark masses suffering from large non-perturbative 
QCD corrections. 

These two problems of the MS scheme can be overcome simultaneously by adopting 
the effective charges ( |B.1|) and (|B.2|) as expansion parameters at 

< \q 2 \ < 0(m 2 z ) (B.12) 

when confronting with experiments. The connection with a high energy theory, e.g. at 
q 2 = m 2 z , can then be made free from light quark mass ambiguities by the use of the 
7 In fact the discontinuity can be evaded by using yet another unphysical effective 



coupling, the so called dimensional reduction DR scheme ||116 | 
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manifestly RG invariant expressions (|B.1|) and ( |B.2| ). In the region the effect lve 

charges at two different q 2 are related by dispersion relations. 

The light hadron (first 5-quark, or "5g") contributions to the differences 



1 


1 


a(q 2 ) 


a 


s 2 {q 2 ) 


s 2 (0) 


a(q 2 ) 


a 



An 



An 



e 2 (q 2 ) e 2 (0) 
An Air 



g 2 (q 2 ) g 2 (0) 



An 
■ An 



Reng(g 2 ) - ng(o) 



Relf T Q Jq 2 ) - IC1(0) 



3Q 



(B.13) 
(B.14) 



have been parametrized in the region < \q 2 \ < m 2 z as follows. For the photon vacuum 
polarization function, we use 

-f(q 2 



An 



n?>) 



5q 



-/( 9 2 ) + ^E[R<(? 2 )-n?,?H 2 

q=c,b 



for - m| < q 2 < 0, (B.15a) 
for < q 2 < m|. (B.15b) 



Here the results of the dispersion integral analyses P7|,|2"8"| are parametrized by 



f(q 2 



' 1.096 ln(l + |g 2 |) 
0.3261 ln(l + 3.927|g 2 |) 
0.2486 + 0.4009 ln(l + |g 2 |) 



for 0.0 < 
for 0.3 < 
for 3.0 < 



3.878 + 0.4084|ln ^ + 0.005696(^ - 1)} for 50 < 
I 1 s s J 



g 2 |(GeV) < 0.3 

GeV) < 3 
g 2 ] ( GeV) < 50 

g2j( GeV) < m z 

(B.16) 



with so = (91.176 GeV) 2 . The parametrization (p.!6| ) is copied from ref. |27] for GeV < 
\J\q 2 \ < 50 GeV and smoothly connected to the most recent estimates of ref. at 
g 2 = m | ; /(m|) = (0.0283 ± 0.0007)/a. In the time-like region (0 < q 2 < m|), the 
second term in ( |B.15b| ) is added in order to account approximately for the threshold 
contributions of the charm and bottom quarks. 

Hadronic contribution to the photon-Z mixing two-point function can then be esti- 
mated as 



An 



Rel4> 2 ) - 11^(0) 



TT 3< 2 



5q 



where 



A^,(g 2 



2vr 



3g 2 



RelT^(g 2 ) - I%%(0)] + A^(q 2 ) + A c (q 2 ) + A b (q 2 ) 



2 ) 

B.17) 



2T(u 



e 1 e 



e 1 e 



4a 2 | m^iml + |g 2 |) rn^rrv^ + \q 2 \) 



(B.18) 
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is an estimate [29| for the extra contribution from the u, d, s quarks, and 

1 



67T 



C q QA ,83(0; m q , m q ) - B 3 (q ; m q , m q ) 



+ 



67T 



B' v (0; m q , m q ) - B' v (q 2 ; m q , m q 



T3.19) 



for q = c and b, are calculated perturbatively. Note that in the m u = rrid = m s limit, the 
identity A u + + A s = holds. Thus, the term gives an estimate [2£] of the flavor 
SU(3) violation effect. Contributions of leptons, the top quark and any other new parti- 
cles, as well as the light 5-quark contributions at \J\q 2 \ > rn z are treated perturbatively. 

The light quark masses to be used in the region | q 2 1 > m| are determined by requiring 
continuity of the two effective charges at q 2 = m|. The left-hand sides of eqs. ( |B.15| ) and 



( |B.17| ) are evaluated perturbatively, and equated with the estimate at q 2 



rrtr 



An 
4ir 



ReU 



■3Q 



n?,?(o) 

T 3Q 



5q 



-f(m 



2 ) 



Ren^(m|) - n£(0) 



5q = -7>f( m z) + A^(m|) + A c (m|) + A 6 (m|) 



where the mean value of the estimate P8 
- /(m 



4) 



-0.0283 ±0.0007 



-3.88 + cW ; 5 



had 



0±0.1 



(B.20) 
(B.21) 

(B.22) 



is taken at m z = 91.187 GeV. Note that the additional term at the right-hand side of 



( |B.15b| ) is less than 0.001 and the discontinuity at q 2 = m 2 z is negligibly small. With the use 
of the expressions ( |A.9a| ) and ( |A.9b[ ), the two matching conditions can be approximated 
by: 



J2Q 

5g 



In 



6 

m d m s 
m 2 



1 m 
-In — 



mi 



^m z ix 



18ttA, 



0.152. 



7r m\ 



(B.23) 
(B.24) 



Taking the charm and bottom quark masses 



m c = 1.4 GeV, 
m b = 4. 7 GeV, 



(B.25a) 
(B.25b) 



and including 0(a s ) corrections one finds for 8h&d = 



a s (m z ) 





0.11 


0.12 


0.13 


m u = md( GeV) 


0.055 


0.089 


0.093 


0.097 


m s (GeV) 


0.064 


0.104 


0.108 


0.113 



(B.26) 
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Our program finds appropriate light quark masses for arbitrary a s (m z ), c\ad ; and 
input values by solving the continuity conditions ( B.20|) and QB.21Q . It should be pointed 



out here that these light quark masses are fixed merely to ensure the continuity of the 
effective charges at q 2 = m| and that they do not have a direct physical significance. 
At \q 2 \ > m|, where those quark masses are used, the mass effects are suppressed by 
m^/m 2 z and never become significant. Whenever the light quark mass values play a 
physically significant role, their values must be chosen independent of those of eq. ( |B.26|) 
by appropriate physics arguments. 

In Fig. |2| the SM predictions for the effective charge 4ir/e 2 (q 2 ) and the effective weak 
mixing angle s 2 (q 2 ) are shown in the region 1 MeV < \J\q 2 \ < 1 TeV for m t = 100, 150, 
200 GeV and m H = 100, 1000 GeV with 5had = 0. The solid lines show the space-like 
(q 2 < 0) effective charge, whereas the dashed lines the time-like (q 2 > 0) effective charge. 
The top-quark effect at q 2 = m| can be parametrized by 

i = 128.71 + ^ + 0.024(1 + 5 ^)(^X) 2 , (B .27) 

a ( m z)sM V tt j \ rnt J 

s 2 (ml) s 2 (0)l 5 had i a s x/100GeV\ 2 

z = -3.09 + — + 0.009(1 + 5—) , B.28 

vr 2 v ?t J \ m t J 



a(m 2 z ) a 



for m t > 100 GeV representing typical contributions of a heavy particle to the running of 
the effective charge form factor a(q 2 ) and s 2 (q 2 ) between q 2 = and m%. 

When constraining new physics contributions the value of a{m 2 z ) is required, but only 
a = «(0), that is, the fine structure constant, is precisely measured. When new physics is 
contributing significantly to the running of the effective charge form factors between q 2 = 
and q 2 = m|, its value can deviate from the SM prediction ( |B.27p . In order to account 



for both such new physics contributions and future improvements in the measurement of 
5had, the parameter 

128.72. (B.29) 



is introduced. For instance, in the SM one finds from (|B.27|) 



[<y S M = 4ad + 0.024(1 + 5^) f 10 ° GeV y - 0.01 . (B.30) 
The last two terms are close to zero for m t = 150-200GeV, such that within the SM : 

sm w S had . (B.31) 
In general, new physics contributions can be accounted for by 

S a = [5 a ]sM + 47r[Ren^(m|)-n?-(0)l M ^ . . (B.32) 
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An example of the extra term is found in ref. [!1|, where consequences of the gauge- 



invariant dimension six operators JT(|] have been studied in detail. 



The MS couping constants e 2 (/i) and g 2 (fi) are determined from the identities eqs. ( |B.1| ) 
and (|B.2|) evaluated at large \q 2 \, say at q 2 = m|. The magnitude of e 2 {m z ) depends on 



m t and the assumed a s (m z ) value, and that of s 2 (m z ) depends also on the s 2 (m 2 z ) value 
as observed at LEP/SLC. For a s = 0.12 one obtains 

-0.12 for m t = 100 GeV 
128.00 + <) had + < |nnfi _onnn„v » ( B - S6 > 



&{m z ) S M~ e 2 (m z ) SM naa ) +0.08 for m t = 200 GeV 

+0.15 for m t = 250GeV 



and 



r 0.0007 for m t = 100 GeV 

.2/ n -2r 2-. a , 0.0009 for m< = 150GeV , nnA . 

s (m z ) SM = s K) + -<W + Q 001Q for ^ = 20Q ^ , (B.34) 

[ 0.0011 for m t = 250 GeV 

The relatively large m t dependences above, as opposed to those of eqs. ( |B.27|) and ( |B.28|) , 



result from the non-decoupling of the heavy top quark due to the logarithmic m t depen- 
dence of the MS renormalized two-point functions in eqs. (p.l|) and (|B.2| ), as explained 
earlier. In the presence of many new particles at the TeV scale, such as in the super- 
symmetric standard model, all new particle contributions are suppressed by their inverse 
mass-squared as demonstrated for a heavy top quark in ( |B.27| ) and ( |B.28| ) for the effec- 
tive charges, while the magnitude of e 2 (m z ) and s 2 (m z ) are affected strongly. One should 
then either adopt the effective light particle theory for the MS couplings |5], |9], [IT], ^3[ or 
use the above effective charges below TeV scale. 

Finally, it is worth mentioning that the expressions for the running of the remaining 
two charge form factors: 

UT^T = Jt^ + ^Tf T 3 z (q 2 ) - 2s 2 ReTf T Q z (q 2 ) + s'ReT^% 2 ) , (B.35) 



' + ReU^ w {q 2 ) , (B.36) 



9w(l 2 ) 9 2 (n) 

are not the exact solution of the one-loop RG equations of the MS couplings, but that the 
0(1) terms at the right-hand sides remain small at all q 2 , provided the renormalization 
condition 

e 2 (/i) = e 2 {Q 2 ) , (B.37a) 
S 2 (/i) = s 2 (Q 2 ), (B.37b) 
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is chosen with 



Q 2 



m 2 z if \q 2 \ <> 0(m 



q 2 if \q 2 \ > m 2 z . 
Therefore, g z (q 2 ) and glyi^q 2 ) can be safely calculated from 

4vr s 2 (Q 2 )c 2 (Q 2 ) , 1 , ^ 



(B.38) 



47T 



^(g 2 ) «(Q2) 4 



(B.39) 
(B.40) 



where the two quantities 



S z (q 2 ) = 16vrRe{4 Q 7 (Q 2 ) - lf T 3 z (q 2 ) - 2s 2 (Q 2 )[lif n {Q 2 ) 



+s\Q 2 )m®(Q 2 ) -Tl Q T % 2 ) 



S w (q 2 



16vrRe 



n^W 2 ) - n^ V y(g 2 ) s -2 =5 2( Q 2) 



(B.41a) 
(B.41b) 



< oo. 



remain small (free of large logarithm) at all q 2 , < | q 2 

In principle, the parametrization ( |B.16| ) can be used to account for the hadronic con- 
tributions to the n T Z (g 2 ) and U T W (q 2 ) terms at \q 2 \ < m z with the help of the CVC and 
PCAC hypotheses. However, we find that the contribution of light hadrons are negligible 
at low momentum transfers \q 2 \ <C m z , and hence the perturbative expressions ( |A.9| ) with 
the light quark masses as obtained by the matching conditions ( |B.20| ) and ( |B.21| ) are used 
when evaluating these functions. 

It is important to note that the expressions ( |B.39| ) and ( |B.40| ) are valid in the sense 
of a perturbative expressions, and therefore the scale Q 2 has been chosen such that the 
Sz{q 2 ) and Sw{q 2 ) terms remain small. The typical scale of the charge form factors g z (q 2 ) 



and Qiviq 2 ) are Q 2 = m 2 z rather than Q 2 
and U parameters then follow 



q 2 for \q 2 \ <C m 2 z . Our definitions of the S 



S = S z (0) = 167rRe 
S + U = S w (0) = 16vrRe 



IT^ Q 7 (m 

riy i7 (m 



=33 



- n T , z (o) 



I) 



;b.42) 

^B.43) 
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App. C SM contributions to S, T, U and 5b(m 2 z ) 

This Appendix deals with the SM contributions to the universal electroweak param- 
eters S, T, U and the Zb^L vertex form factor ^(m^) used as free fit parameters. The 
complete analytic formulae are given at one-loop level and the two-loop corrections are 
also included as far as they are known. We adopt the perturbative order aa s [^6j|5^ , |5^ -|5^ 
corrections at a s = ct s { m z)Ms m evaluating the S, T, U and 5b(m 2 z ) parameters, since it 
allows the readers to reproduce our results unambiguously and straightforwardly. The ef- 
fects due to non-perturbative threshold corrections [§3H§3 snou ld be evaluated carefully, 
and one can obtain more precise predictions of the SM from our formulae by adjusting 
the effective top-quark mass to produce the same S, T, U and 5f,(m|) values. 

CI Ssm 

The S parameter in the SM can be expressed as a sum of three pieces: 

Ssm = Si + S q + Sb , (C.l) 

where the indices denote contributions from the leptonic, hadronic and the bosonic (that 
is, W, Z, H) sectors of the SM, respectively. Each term is separately finite. St and Sb are 
given at one- loop order, whereas the hadronic contribution S q with the two- loop 0(aa s ) 
correction [|52],[55|,|56|]. 

The leptonic contribution is a sum of three terms 



i=i 

where each generation (z/j,£j) contributes 



GUx) = -Ulnx + (l + 5x)A(x) - 10a;}. (C.3) 
6 ' 

The real function A is 



1 + VI - 4x . n 1 



2y/l - Ax In -= for < x < 

A(x) = \ 1 4 (C.4) 

2V4x — 1 tan 1 = for x > — . 

V4x - 1 4 

For the case of charged lepton masses much smaller than the Z mass one finds 

G e s (x) = (2 + hnx)x + 0{x 2 ) . (C.5) 
The [y T , r) doublet contribution is hence St ~ —0.0002. 
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The hadronic contribution calculated up to 0(aa s ) two-loop level is fl52|, p5|. |56 



S q = Sf + S® . (C.6) 
The one-loop contribution is again a sum of the three terms 

^9 r q ( 1 J^h 



sf^-EGO,^ . (C7) 

with C g = 3, where each quark generation contributes 

y) = 4{ ln - - (1 + llx)A(x) + (1 - 7y)A(y) + 22a; + 14y} . (C.8) 
When both quarks are light as compared to Z, one finds 

G Q s (x, y) = ( \ + \ lnx)x + ( \ + \ In y)y + 0(x 2 , y 2 ) . (C.9) 
When only the down- type quark is light (y <C 1), one finds 

y) = — { 22x - lnx - (1 + llx)A(x) + 3y(4 + 31ny) + 0(y 2 )} , (CIO) 
18 ^ ^ 

= l{_ lnx _I + _L + 3y ( 4 + 31ny) + 0(1 y 2 )} . (C.11) 
18 L 6 20 x or J 

For large m t , the quark contribution 5^ becomes negative and its magnitude grows loga- 
rithmically. 

The two-loop 0(aa s ) correction [52, ^5 , 56] can be expressed in terms of the B' v and 
B' A functions ( |A.12| ) of appendix A : 

= C q ^ £ { 1 B' v (m 2 z ; u h Ui ) + 1 ^(m|; d,, d,) 

--^ B A( m h u ii u i) - -^ B A( m h d ii d i)} > (c.i2) 

where the quark label stands for its mass as in appendix A. It is easily seen that the right- 
hand side of the above equation is independent of the unit-of-mass /i for each generation, 
and that they are in fact a function of the ratios m^./m 2 ? and m 2 d Jm 2 z . The contributions 
from the first two quark generations are again negligible. The two-loop term S^> is hence 
dominated by the (t, b) doublet contribution, which can be approximated by 

=°.p{-s b "' + i< i+ 5 + <)}' (CU3) 
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with x t = m 2 /m 2 z . The expression QC.13 ) agrees with ref. [52|. The following table shows 
the full hadronic contribution S q for several values of m t in lowest order (a s = 0) and 

0.12: 



with the 0(aa s ) corrections for a 



m t ( GeV) 


s q 


a s = 


a s = 0.12 


100 


-0.008 


0.010 


120 


-0.033 


-0.017 


140 


-0.052 


-0.038 


160 


-0.069 


-0.055 


180 


-0.083 


-0.070 


200 


-0.095 


-0.082 



(C.14) 



using for the bottom mass mj = 4.7 GeV. The two-loop correction is important for 
relatively small m t values. 

The bosonic contribution is expressed as 

,2 ,^2 



7i I v m 7 ' v mi 



where 



Fo(c 2 ) 



1 i 2 7 14 2 n 2 i ] A{ 2 , 

mc c + - c Ac 1 

12 8 3 V3 



12 



# s (x) 



-1.451 for c 2 
3 1 2 

- X X + 

8 12 



80.24 \ 2 
91.187 



3 — x x 2 

+ 777 + 



4 



24 4 1-i 



xlnx + 



x x 
3 + 12. 



Here v4(x) is given by (|C.4|) , and 



B(x) 



For large m H one has 



r , _ a 

yx(4 — x) tan 1 y 1 



x(x — 4) In 



x 
2 



\fx + \Jx — 4 



for < x < 4 



for x > 4 . 



i3 



J- j >-" 1.1 lib g _ 

12 n 7^| + 36 _ 48m^" + 



37 17 mi 



m 



The total bosonic contribution S# is tabulated below for several m H values 



m H ( GeV) 


Sb 


50 


-0.234 


100 


-0.166 


200 


-0.107 


400 


-0.061 


1000 


-0.008 



(C.15) 



(C.16a) 



B{x). (C.16b) 



(C.17) 



(C.18) 



(C.19) 
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C2 T SM 

The T parameter in the SM can be expressed as a sum of three individually finite 
pieces: 



SM 



Te + T a + T) 



B 



(C.20) 



where the indices denote the leptonic, the hadronic and the bosonic (that is, 7, W, Z, H) 
contributions, respectively. Te and Tb are evaluated at one-loop order, whereas T q contains 
irreducible two- loop contributions in aa s order [p^,pl|-|o^[] and in the m\ order B^ , |57 , 
Reducible higher-oder contributions are taken account of by the identity ( ^.36a| ). 
The leptonic contribution is a sum of three terms 



2y/2ir 2 a ^ V ' mi 



(C.21) 



where 



G T (x,y) 



x + y 



xy 



2(x-y) 



x 



(C.22) 



The leptonic contribution of the first three generations is hence negligible; even the iv T , r) 
doublet contributes to Tg only about 0.00005. 

The hadronic contribution is calculated including the 0(aa s ) |52|,p^-|56| and the irre- 
ducible 0(mf) ]|7|,|57|,0 two-loop corrections: 



T Q = T (0) + T« + rj 2) . 



(C.23) 



The one-loop contribution is a sum of the nine terms 



T (o) 



GF ' nl C, ± |K„ 



2V2 



G 1 



tc* a 



m 



2 ' 



(C.24) 



with C q = 3. The function Gt(x, y) is found in (|C22|) . In the limit of the diagonal KM 
matrix elements Vy = the contributions from the light quarks of the first and second 
generations can be neglected. 

The two-loop contributions are only important for the t-b doublet : 



a s 3 + tc 2 



G F m\ 



2^2 



tc* a 



T^ = C q a { , 
1 \Amr 



2 / n m 2 \ 2 



2^ 



tc* a 



(C.25) 
(C.26) 
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where terms of order m\jm\ are neglected. The function [m H /m t ) gives [47| : 
p^\m H /m t ) = -0.74,-4.72,-6.95,-11.70,-10.74, for m H /m t = 0,|,1,5,10. The nu- 
merical value of the 'expansion parameter' in the above expressions is GFm 2 z /2\/2ii 2 a = 
0.4761. 

The following table shows the contributions from each term in eq. ( |C.23|) for several 
values of m t , the lowest order contribution and the 0(aa s 



contribution with 



a. 



0.12, and 0(mf) contribution with 



m 



H 



100, 1000 GeV: 



mt ( GeV) 


T (0) 
9 


<? 


T( 2 ) 
i 


a s = 0.12 


m H = 100 


m H = 1000 


100 


0.419 


-0.047 


-0.003 


-0.005 


120 


0.607 


-0.068 


-0.006 


-0.011 


140 


0.830 


-0.092 


-0.010 


-0.020 


160 


1.087 


-0.120 


-0.016 


-0.035 


180 


1.379 


-0.152 


-0.024 


-0.055 


200 


1.705 


-0.188 


-0.034 


-0.084 



using for the bottom mass m& = 4.7 GeV in T®. 
The bosonic contribution is 

Gf tv? z 



b 
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7r z a 
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where 
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and s 2 = 0.2312 
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ha(x/c 



2\ 1 
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c 2 = 



m 



w 



TO; 



For a heavy Higgs boson (m 2 H ^> m 2 z ) one finds 



Hn 



( m \) 


3 ■ 


\m z J 


~ _ 4 . 



(l_ c 2 )ln !^ + c 2 lnc 2 + ( (1 _ c 4 ) ln !^ + c 4 lnc 2|!4 +0 



The total bosonic contribution Tb is tabulated below for several m H values : 



(C.27) 



(C.28) 



(C.29) 



(C.30) 



rrfr 



m 



(C.31) 



m H ( GeV) 


T B 


50 


-0.227 


100 


-0.257 


200 


-0.314 


400 


-0.396 


1000 


-0.529 



(C.32) 



93 



C.3 Usm 

The U parameter in the SM can be expressed as a sum of three pieces: 

U su = U e + U q + U B , (C.33) 

where the indices denote the leptonic, the hadronic and the bosonic (that is, 7, W, Z, H) 
contributions. Each term is separately finite. Ug and Ub are given at one-loop order, 



whereas the hadronic contribution U q is given with the two- loop 0(aa s ) correction [ |5i2| 
The leptonic contribution is a sum of three terms 

^ = -E^(0,^|V (C.34) 



7T S \ m z 



where the contribution of each generation (fi,£i) is 

Gu(x, y) = -^±1 - ^A(x) - l -^-A{y) + /^(x/c 2 , y/c 2 ) . (C.35) 



Here A(x) is given in ( |C.4j ), and 



a; + y (x - yf {x - y) 4 - 3 (x 2 + y 2 ) y 

trjix.y) = ; r In — 

JuK 12 6 I2(x-y) x 

(x-y) 2 + x + y-2 

P(x, y)L(x, y) , (C.36) 



6 



with 



and 



L{x,y) 



(3(x, y) = ^\l-2(x + y y + (x- y y\ , (C.37) 



i, i-x-y + p(x,y) , 1 1— 1 — 1 /-, /-^i 

-In- — r for \y/x- y/y > 1 or y/x + y/y < 1, 

2 1 -x-y- /3(x,y) 

, 1 — x + y y . — 

^ ^^ r + tan" for \^-^[<K^+^y. 



In the limit of vanishing lepton mass one has 

GV(0,0) = ^lnc 2 , (C.39) 

for c 2 = m^z/rnl, and hence the contribution of the first three lepton generations is 

1 / 80 24 \ 2 

U t « - In c 2 = -0.0814 for c 2 = . (C.40) 

7T V91.187/ v ; 
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The hadronic contribution is calculated up to two- loop 0(aa s ) level fl52|, p5|, |56| : 

U q = C/f + 17« . (C.41) 
The one-loop contribution is again a sum of the three terms in the limit of diagonal KM 



matrix elements Vij = 5if. 



^(o) = C, j, ^ 
11 i=i 



mt . m 



d, 



1 ' 2 

m z m z 



(C.42) 



with C q = 3, where Gjj{x,y) has been given above in (|C.35|) . For the first two quark 
generations the approximation (|C.39| ) holds. For the contribution of the (t , b) doublet 
the following approximation is useful: 



Gu(x,y) = - lnx + — - -y \ny + — - -c + -y)- + 0[—,y 



6 



12 



1 1 



40 8 



x 



(C.43) 



The two-loop 0(aa s ) correction [52, 55,56] can be expressed by 



(C.44) 



i=l 



where the two-loop function Gjj is given in terms of the By and B A functions ( |A.11| ) of 
appendix A: 



G { u\m u ,m d ) = Re^-^[B' V + B' A ](m 2 z ; m u , m u ) + j^[B' v + B' A ](m 2 z ] m d ,m d ) 

- - [B' v + B' A ] {m 2 w ;m u ,m d ) 
o 

.(1) : 



(C.45) 



It is readily seen that the function G v is independent of the unit-of-mass fi. The contri- 
bution of the first two quark generations can be approximated by 



G[ 7 1) (0,0) = ilnc 2 



1 



(C.46) 



just like in eq. ( |U.39| ) for the one-loop contribution. The top-bottom contribution can be 
approximated by 



In 4 

m z 




MO 


+ 


L 

3m| 


Am 2 




2 


3m^ 


mi 



m: 



Am 2 
-F 1 (0) 
4 



Ami 



-Mo) 



m. 



(C.47a) 
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The expressions ( |C.46|) and (|C.47| ) agree with ref. [52]. The following table shows the 
total hadronic contribution U q for several values of m t in lowest order (a s = 0) and with 
0(aa s ) corrections for a s = 0.12: 



m t ( GeV) 



100 
120 
140 
160 
180 
200 



U a 







-0.118 
-0.034 
0.029 
0.079 
0.122 
0.159 



O.s 



0.12 



-0.148 
-0.057 
0.009 
0.063 
0.108 
0.147 



(C.48) 



using for the bottom mass m& = 4.7 GeV in evaluating the lowest order Gy , while the 
0(aa s ) correction Gjp is calculated in the limit of vanishing bottom quark mass. 
The bosonic contributions are given as 



7i I v mi 7 v mi J J 



(C.49) 



where the constant term Fu{rn'^/m 2 z ) is found to be 



1 

+ "7 



1 /l 



s' + -7 r 



1 



9 7 



1-c 2 c 4 6cV 4 V3 1-c 2 c 2 c 4 2 c' 
29\ f /_ 4 1 \ 2 s 2 /_ 2 3 

rr 



lnc 2 



+ (, + f)^, + { ( 6 + i4 



-3 



2 \s 



3 c 2 / c 2 



2 c 2 



3 
3 



2 



'4 



19 

8c 2 + J^ 4 



(C.50) 



, / 80.24 \ 2 , „ 
1.043 for c 2 = 7 —r^z and s 2 = 0.2312 , 
V91.187/ 



and the m H dependence is given by 



Hu{x H ) = -H s {x H ) + H s { ^ I = 



ni 



2 ' H 



(C.51) 



The function Hs(x H ) is defined in eq. (|C.16b|) . In the large mass limit (m H ^> m z ) the 
leading logarithm (\nm H ) of the function H$ (see eq. (|C.18|) ) cancels in Hu of eq. (|C.51|) , 
and hence one has 

MSK-'-^-^Q)- (c - 52) 

Note, however, that the m H dependence is very small as seen from the table below showing 
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the total bosonic contribution U q for several m H values: 



m H ( GeV) 


U B 


50 


0.345 


100 


0.344 


200 


0.341 


400 


0.340 


1000 


0.339 



(C.53) 



C.4 5 b (m 2 z ) S M 

The Zb L b L vertex form factor 5b(m|)sM in the SM is expressed by: 

UO = ^Vl) + ^Vl) + SVl) • 

The one-loop contribution 

4° VI) = It (m|) + c 2 r^(m|) + r^(m|) 

is calculated using the vertex functions of appendix A, which can be approximated by 

'm t + 36 x - 



(C.54) 



(C.55) 



-0.00076 - 0.00217 



(C.56) 



V 100 

for 100 GeV < m t < 250 GeV. The second term at the right-hand side of eq. ( p.54| ) is the 
0(a s m 2 ) two-loop contribution [f46| : 

^ 2 G F m\ 



7T v O 

The last term is the 0{mf) two-loop contribution 



8V2 



7T Z 



GfVTI 

8^2 TP 



2 \ 2 



-(2), 



m H /m t ) 



(C.57) 



(C.58) 



where the function T^ 2 \m H /m t ) is given in ref. [57]. For m H /m t = 0, |, 1,5, 10, it gives 
T^(m H /m t ) = 5.71, 2.46, 1.47, 3.69, 7.92. 

The following table shows the contributions from each term in eq. (C.54), (^(m^), 
^^'(ml) with a s = 0.12, and S^\m 2 z ) with m H = 100, 1000 GeV, for several m t values: 



mt ( GeV) 


SVl) 




jp 


(m|) 


a s = 0.12 


= 100 


m H = 1000 


100 


-0.00481 


0.00018 


-0.00000 


-0.00002 


120 


-0.00603 


0.00026 


-0.00001 


-0.00003 


140 


-0.00746 


0.00036 


-0.00002 


-0.00005 


160 


-0.00908 


0.00047 


-0.00003 


-0.00007 


180 


-0.01089 


0.00059 


-0.00005 


-0.00010 


200 


-0.01285 


0.00073 


-0.00009 


-0.00013 



(C.59) 
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App. D One loop scalar functions 



In this appendix explicit analytic expressions for the B functions are given, as well as 
the reduction of higher C and D functions to Co and Do functions. 

D.l A, B, C and D functions 

Following Passarino and Veltman [3^] the A, B, C and D functions are defined by : 

A{ mi ) -- 
[B , £?H (ij) = 

[c , c», cn (ijk) = 

[Do, £>", D^] (ijk£) = 



dPk 1 

Wk [1, fc**, W] 



N l N j 



d D k [1, jf, fc^fc"] 



iTT 2 A^iVfc 

Wk [1, fc", A^Jfe"] 



m 2 NiNjNkNe 



(D.l 
(D.2 
(D.3 
(D.4 



where £) = 4 — 2 e, 



d D fc = r(l-e) (ir/j, 2 ) e d D k 



is the MS regularization |7T), 117| , and the propagator factors are 



N x = k 2 - ml + ie , 

N 2 = {k+pxf -m 2 2 + ie , 

N 3 = (k + px+ P2) 2 - ml + ie , 

Ni = (k + px + P2 + Pz) 2 - ml + ie . 



(D.5 

(D.6a 
(D.6b 
(D.6c 
(D.6d 



(D.7a 
(D.7b 

(D.8a 



The vector/tensor functions are reduced to scalar functions as 

£"(12) =p?fli(12), 
B^(12) = Kp^ 21 (12) + ^£22(12) , 

for the two-point functions, 

C**(123) = rfCn(123) +^Ci 2 (123) , 
C^(123) = p^C 21 (123) + p&£C 22 (123) +pfe } C 23 (123) +^C 24 (123) , (D.8b 

for the three-point functions, and 

D"(1234) = rf J D 11 (1234) + p 2 i J D 12 (1234) + p 3 1 J D 13 (1234) , (D.9a 
L>^(1234) = p^£> 21 (1234) +^D 22 (1234) +p^L> 23 (1234) 

+pS^2 } ^24(1234) +pl^ } D 25 (1234) + jft£ } £> 28 (1234) 

+^D 27 (1234) , (D.9b 
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for the four-point functions. Higher rank tensor functions do not appear in our applica- 
tions in the 't Hooft-Feynman gauge. 

The scalar functions B iy Ci, Di, are defined by 

B < (12) = J B i (p?;mi,m 2 ) (D.10) 

for i = 0,1,21 and 22, 

(^(123) = Qiplpl, (p! + p 2 ) 2 ;m 1 ,m 2 ,m 3 ) (D.ll) 
for i = 0, 11, 12 and 21-24, and 

A (1234) = Di(pl,pl,pl, (P1+P2+P3) 2 , (P1+P2) 2 , (P2+P3) 2 ;rn u m 2 ,m 3 ,m 4 ) (D.12) 
for 2 = 0, 11-13 and 21-27. The basic scalar functions Bo, Cq and Dq were obtained by 



't Hooft and Veltman | 118 |. The Fortran code FF | 119 | is used for the general form of 
Co and D functions. Reductions of higher B, C, D functions are given in the following 
subsections. 

D.2 B functions 

It is convenient to introduce the following four scalar B functions in addition to Bq 
and B\ above : 

B 2 (q 2 ;m 1 ,m 2 ) = B 21 (q 2 ;m 1 ,m 2 ), (D.13a) 

B z (q 2 ;m 1 ,m 2 ) = -B 1 (q 2 ;m 1 ,m 2 ) - B 2 (q 2 \m x ,m 2 ), (D.13b) 

B,(q 2 ; mi ,m 2 ) = -mfB^q 2 ; m 2) mi) - m^g 2 ; n»a), (D.13c) 

B 5 (q 2 ; mi ,m 2 ) = A{m x ) + A(m 2 ) - 4 B 22 (q 2 ; m u m 2 ). (D.13d) 



All two-point functions of the standard model and its supersymmetric extension [|120|| are 
expressed compactly in terms of the above six B n functions (n = 0,1,. .,5) being only 
logarithmically divergent. The ultra-violet singular factor is parametrized by : 

A = ^ + ln/x 2 . (D.14) 

In the MS (modified minimal subtraction) renormalization scheme the singular piece A 
in these functions is simply replaced by a logarithm of the unit of mass [i: 

A^ln/i 2 . (D.15) 
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The six B n functions are then expressed by 



B (q ; mi, m 2 ) = A - F (q ; m l , m 2 ), 
B 1 (q 2 ;m 1 ,m 2 ) = --A + Fi(g 2 ; mi, m 2 ), 

^2(g 2 ; mi, m 2 ) = ^A - F 2 (g 2 ; m 1 , m 2 ), 

5 3 (<? 2 ; mi, m 2 ) = ^A - F 3 (q 2 ; m 1 , m 2 ), 

-B 4 (<? ; mi, m 2 ) = A - F A {q ; m 1 ,m 2 ) 

q 2 

B 5 (q 2 ; mi, m 2 ) = yA - F 5 (q 2 ; m 1 , m 2 ), 



where the finite parts F n have the following Feynman parametrizations : 

F (q 2 ; mi, m 2 ) = / dxlnH, 
Jo 

F 1 (q 2 ;mi,m 2 ) = / dxxlnH, 
Jo 

F 2 (q 2 ;mi,m 2 ) = dx x 2 \nH, 
Jo 

F 3 (q 2 ;mi,m 2 ) = / rfx x(l — x) IniJ, 

JO 

FAq 2 ] mi, m 2 ) = / [(1 — x)m 2 + xm 2 ] ln_ff, 
Jo 

F 5 (g 2 ;mi,m 2 ) = / rfx [(1 - 2x)(m\ - m 2 ) + (1 - 2x) 2 q 2 ] InH, 
Jo 



with 



H = 



(D.16a) 
(D.16b) 

(D.16c) 

(D.16d) 

(D.16e) 

(D.16f) 

(D.17a) 
(D.17b) 
(D.17c) 
(D.17d) 
(D.17e) 
(D.17f) 

(D.18) 



(1 — x)m\ + xm 2 , — x(l — x)q 2 — it . 

Among the six F n functions four (n = 0, 3, 4, 5) are symmetric under the exchange of the 
two masses. It is useful to introduce the antisymmetric F function 



F A (q 2 ;m 1 ,m 2 ) = Fi(g 2 ; m 2 , m^ - F 1 (q 2 ;m 1 ,m 2 ). 



(D.19) 



In terms of the two symmetric functions F and F 3 and the antisymmetric function Fa 
all the remaining F n functions can be expressed compactly : 



F 1 (q 2 ;m 1 ,m 2 ) = - [F - F A ] (g 2 ;mi,m 2 ), 



F 2 (q 2 ;m 1 ,m 2 ) 
F 4 (g 2 ;mi,m 2 ) 



(q ;mi,m 2 ), 



-(F — F A ) — F 3 

m\ + m\ m\-m\ 

fo H 7, ?A 



(g 2 ;mi,m 2 ), 



(D.20a) 
(D.20b) 

(D.20c) 



2 u 2 

F 5 (g 2 ;mi,m 2 ) = q 2 (F - 4F 3 ) + (m 2 - m 2 2 )F A (g 2 ;mi,m 2 ). (D.20d) 
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Therefore it is convenient to give closed analytic expressions for the three functions, F , F 3 
and Fa'- 



F (q 2 ; mi, m 2 ) = ln(mim 2 ) — 5 In — - — 2 + f3L, 

mi 

2 1 3a — 2<5 2 m 2 

F 3 (q ; mi, m 2 ) = - hi(mim 2 ) 5 In - 

o o mi 1 



5 a-5 2 



l + a-25' 
6 



F A (q 2 ; mi , m2 ) 
where 



[a-5 2 )\n — + 5CI-/3L), 
mi 



a 
5 



1 2 

m\ — m 2 



{l-2a + 5 2 )-2 
i(2o-5 2 -l)^ 



for q 2 < (nil — rn 2 ) 2 or <f > (mi + m 2 ) 2 
for (mi — rn-i) 2 < q 2 < (nil + m 2 ) 2 , 



(D.21a) 

-PL, 

(D.21b) 
(D.21c) 



(D.22) 
(D.23) 
(D.24) 



and the function L is defined as 

,2 



L(q ;mi,m 2 ) 



2 l-B-a 
2 111 l-/3-cr 



— «7r for g 2 > (mi + m 2 ) 2 , 

for q 2 < (nil — i~n 2 ) 2 , 



(D.25) 



7 (tan 1 + tan 1 t#) for (mi — m 2 ) 2 < g 2 < (mi + m 2 ) 2 . 
Also the derivatives Fq, F% and are needed for certain applications. One finds 



F'(q 2 ; m x , m 2 ) = -A 1 + Sln^ - (5 2 - cr)§} 
(7 L mi p J 



F3(g 2 ;mi,m 2 ) 
F' A (q 2 ;mi,m 2 ) 



i ri 



(D.26a) 



g 2 16 



(7 



(a - 25 2 ) In — - 25 + 5(1 - 3cx + 25 2 " L 



(D.26b) 
(D.26c) 



m x ' j3 . 

The phase factors in f3 and L in eqs. ( |D.24|J[X2^| ) are required to obtain correctly the 
ratio L//3. In terms of the above three functions all the other F^ functions are expressed 
compactly: 



F[(q 2 ;mi,m 2 ) = ~[Fq - F A ](q 2 ;m 1 ,m 2 ), 



F 2 (q ;m 1 ,m 2 ) 



[^-^-^](g 2 ;mi,m 2 ), 



(D.27a) 
(D.27b) 

2 -„ . 2 (D.27c) 
F^(q 2 ; mi, m 2 ) = [F - 4F 3 + g 2 ^ - 4F 3 ') + (m 2 - m 2 )F^](g 2 ; m x , m 2 ), (D.27d) 



Fi(g 2 ;mi,m 2 ) = [^^ + ^^(9^1,^ 
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The derivative of the B Q functions is found to be : 

B' n (q 2 ; m 1 ,m 2 ) = -F^(q 2 ; m 1 ,m 2 ) for n = 0, 2, 3, 4 
B' 1 (q 2 ;m 1 ,m 2 ) = F[(q 2 ;m 1 ,m 2 ), 
B' 5 {q 2 ; m u m 2 ) — -A - F^{q 2 ; mj, m\). 

D.3 C functions 

The higher C functions 

d = Ci(pl,p 2 2 , (pi +p 2 ) 2 ;m 1 ,m 2 ,m 3 ) 
for i = 11, 12, and 21-24 are given in terms of the Co, B and B\ functions as 



where 



and 



C\ 2 



x- 1 



sh - 4 i3) + AO) 



Bf ] + ACq 



r - 1 _i_ 1 r( 23 ) . m i r & r $ 2 r 

<_/24 — T + T-^O ' — 7T 7 7 ( - / 12 



C 2 \ 

C*23 / 

C23 
C22 



/ d(13) 1 75(23) \_ f (-1 o r 1 

-1 £>\ + -DQ + Jl^ll — ^24 



(12) 



-1 



51 



(13) 



Si 13) + AC n 

fif^ + AC12 
- 5^ ^ + j 2 C\ 2 — 2 C24 



/ 2p? 2p lP2 \ 
* ~ \ 2p lP2 2p 2 ) ' 



(D.28a) 
(D.28b) 
(D.28c) 



(D.29) 



(D.30a) 
(D.30b) 
(D.30c) 

(D.30d) 
(D.31) 



Pi 



r 2 2 

fi=m 2 - m 1 
f 2 = m\-m\- (pi + p 2 ) 2 + p? 



(D.32a) 
(D.32b) 



Here the shorthand notations 



B\ 12) = B i (p 2 l ;m 1 ,m 2 ) , 

3) = Bi((pi + p 2 ) 2 ; mi, ma) , 



5f^i^2>2,m3), 



(D.33a) 
(D.33b) 
(D.33c) 



are used for B and B\ functions. 
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D.4 D functions 

The higher D functions 

A = AOi,P2,P3, (pi +P2 + P3) 2 , (pi + P2) 2 , (pi + p 3 ) 2 ;mi,m 2 ,m 3 ,m 4 ) 



(D.34) 



for i = 0, 11-13 and 21-27 are expressed in terms of the Do, Co and Bq functions as 
follows: 



where 



and 

















= m\D 


(D 21 \ 












(D 24 \ 




D 22 








(D 25 ) 




D 26 




\D 23 J 





Co 



(124) 



-Co 



(134) 



+ /2A 



y(234) 



fxD 11 + f 2 D 12 + f 3 D 



1 

/c| 34 )+cr 2) +/iAi-2D 27 \ 



13 



,-r(124) 
°11 

V c^ 23) 

(cW-CgV + hDu 



^(124) 

. ^(123) 
\ °12 



2D 



27 



'12 
^(124) 
°12 



Cn 3 ^ + /2 -D12 

r (124) , n 
^12 +73-^12 

U 12 + J 2 -^13 

c 1 ( 2 24) + /3 A3-2 J D 27 y 



/ 2p\ 2pip 2 2pip 3 
2pip 2 2p\ 2p 2 p 3 



\ 2p!p 3 2p 2 p 3 



2p|. 



(D.35a) 
(D.35b) 
(D.35c) 

(D.35d) 

(D.35e) 



(D.36) 



ft 


to to 


— 


h 


CO to 


- m\ 


h 


= ml 


-m\ 



(pi +P2) 2 +p 2 



(D.37a) 
(D.37b) 
(D.37c) 



The higher C functions in eqs. ( |D.35| ) are written in terms of the Co and Bq functions in 
analogy to the previous subsection : 



c1i 23) \ 
c[fV 



' 2p 2 2p!p 2 \ (B { 
2 PlP2 2p\ J \b, 



(13) 


(12) 



Bt ] + fiCo 



(123) ~ 



^ 13) + hCh 



(123) 



2p\ 2p 1 (p 2 +p 3 y 

2pi(p 2 +p 3 ) 2(p 2 +p 3 ) 2 



B, 



(14) 


(12) 



+ (f 2 + / 8 ) C 



(124) 



(D.38a) 



,(D.38b) 



103 



^(134)' 



'11 



£.(134) 



'12 



£,(234) - 
£(234) 



'12 



with 



' 2( Pl +p 2 ) 2 2{p 1 +p 2 )ps 
2(pi +p 2 )P3 



2j9 2 P3 



2P2P3A 

2p| ) 



2p\ 



.B, 



(14) 



(13) 



Bt ] + (fi + / 2 ) CJ 



(134) ~ 



B^+fsC^ 



(134) 



,(34) 



5 f ) _ + (/ 3 + 2pm)Q 



,(24) 



(234) 1 
3 

(234) 



5, 



(12) 



B 



Boipj; m 1 ,m 2 ) , 
B ((Pi + p 2 ) 2 ;m 1 ,m 3 ) , 
#o((Pi +P2+ P3) 2 ] mi, m 4 ) , 
B (pl;m 2 ,m 3 ) , 

#o((P2 +P3) 2 ;^2,^4) , 

(34) = %^3,m 4 ). 



(23) 




Br = 



, (D.38c) 



(D.38d) 



(D.39a) 
(D.39b) 
(D.39c) 
(D.39d) 
(D.39e) 
(D.39f) 
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Table captions 

Table 1 Universal electroweak parameters of the spontaneously broken generic SU(2)l x U(1)y 
theory. Column 2 list the 10 universal parameters, the masses of the weak bosons 
and the 4 charge form factors at two q 2 scales, and m 2 z . Column 3 contains three 
precisely measured quantities (the fine structure constant a, the Fermi coupling con- 
stant Gp and the Z boson mass mz) together with their relation to the universal 
charge form factors. The factor 5 G is explained in the text: see eq. ( j2.27| ). The last 
column lists those parameters which are used in fits. The 'star' marks parameters 
for which no direct experimental information is available. 

Table 2 Three types of fits are considered. For each sector the free parameters are listed. 

External parameters in the fits in addition to the precisely known fine structure 
constant a, the Fermi coupling constant Gp and the Z boson mass m z are listed 
separately. The quantity 5 a is defined as 5 a = l/a(m z ) — 128.72 Q2.31]) . The 
parameters S, T, U are defined in eqs. ( |2.33|) and 5b in eq. Q2.22Q . 



Table 3 Vertex form factors T{(q 2 ), T f 2 {q 2 ) and r{(g 2 ) in the SM at q 2 = m 2 z . The definitions 
of the form factors are given in eq. (|2.17|) and their explicit forms in appendix A. 2. 



Table 4 Box form factors B(e a , f a ) = B^f a (s,t) for the process e a e^ — > f a f a in the SM at 
s = —It = m 2 z . The definitions of the form factors are given in eq. ( ]2.24j ) and their 
explicit forms in appendix A. 3. 

Table 5 Partial and total Z widths in MeV units in the minimal SM for m z = 91.197 GeV, 
a s (m z ) = 0.12, 5 a = 0, 5 G = 0.0055, m b = 4.7 GeV and m c = 1.4 GeV. See 
eqs. ( |3.15| )-( |3.18| ) in section 3.1 for details. 



Table 6 Summary of all the electroweak data used in the fit, and the fit results. The Z 
boson parameters are studied in section 4.1, the low energy neutral current experi- 
ments are studied in section 4.2, and the charged current experiments are studied in 
section 4.3. In addition, we use the fine structure constant a datum which fixes the 
charge form factor e 2 (0). y 2 h as been calculated by taking the model-independent 
parametrizations of the original data as the inputs of our analysis: eqs. (|4.1| )- (|4.2j ) for 
the Z parameters (section 4.1), eq. (|4.15|) for the v^-q scattering experiments (sec- 



tion 4.2.1), eq. (|4.19|) for the u^-e scattering experiments (section 4.2.2), eq. (|4.22|) 
for the atomic parity violation experiments (section 4.2.3), eq. ( (4.33 ) for the e-D 



polarization asymmetry measurements (section 4.2.3), and eq. ( f4.38|) for the W mass 
measurements (section 4.3). 
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Table 7 The SM predictions for the electroweak parameters. The column 'no-EW is ob- 
tained by dropping all radiative corrections except in the running of d(q 2 ) and s 2 (q 2 ) 
due to light quarks and leptons. The column 'IB A' is obtained by dropping all ver- 
tex and box corrections except 8 G . In both 'no-EW and 'IBA' cases, corrections 
due to the tree-level propagator effects and the external QED/QCD corrections are 
kept. When the predictions depend on a s (mz), we show three representative cases 
for a s (mz) =0.11, 0.12 and 0.13 from top to bottom. The x 2 values are obtained by 
taking account of the correlations among the errors that are presented in the text 
(see section 4). The total number of the data is 22 by counting also (a, G F , m z ), 
while the above three parameters are used as inputs of the SM analysis. The degree 
of freedom of the fit is hence 22 — 3 = 19. 
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Table 1 



Electroweak 
propagators 


Universal 
parameters 


Precisely known 
parameters 


Fit 

parameters 


7 w\0wv 7 


e 2 (0) e 2 (m|) 


a = e 2 (0)/4vr 


* 


7 m0wv Z 


s 2 (0) s 2 (m|) 




s 2 (0) s 2 (m|) 


Z m0wv z 


™z #1(0) 9 2 z ( m z) 


m z 




jy m0wv 


m w 9w(fy 9w( m z) 


AV2G F = %P (1 + * G ) 





Table 2 



Experimental 
inputs 


6-parameter fit 


4-parameter fit 


2-parameter fit 


a, Gp, m z 


(input) 


(input) 


(input) 


low energy 
neutral currents 


AO) 9l(0) 




S T 


5 a 


m t m H 


6 a 


Z parameters 


s 2 (™|) 9z(m 2 z ) 5 b (m 2 z ) 


a s 


S T S b (m%) 


5 a a s 


m t m H 


5 a a s 


m w 






S T U 


5 a 


m t m H 


5 a 


number of 
fit parameters 


6 




4 




2 




number of 
external parameters 




1 




2 




2 



Table 3 



/ 


r{ 


1 2 


1 3 




vl 




0.00252 + 0.00431 i 


-0.00680 - 0.00565 % 










0.00185 + 0.00325 i 


-0.00680 - 0.00565 i 










0.00020 + 0.00032 % 




















0.00203 +0.00354i 


-0.00680 - 0.00565 i 






Ur 




0.00009 +0.00014i 
















d L 




0.00225 + 0.00389 i 


-0.00680 - 0.00565 i 






da 




0.00002 +0.00004i 














b L 


(m t = 


100) 


0.00176 + 0.00107i 


-0.00402 + 0.00000 i 


-0.00347 + 0.00000 i 


b L 


(m t = 


150) 


0.00141 + 0.00107i 


-0.00261 + 0.00000 % 


-0.00763 + 0.00000 % 


b L 


(m t = 


200) 


0.00126 +0.00107i 


-0.00179 + 0.00000 % 


-0.01270 + 0.00000 i 



Table 4 



/ 


sB(e L J a ) 


sB(e R J a ) 




vl 





00109 


+ 


00000 i 


-0 


00006 


+ 


00000 i 




£l 


-0 


00005 


+ 


00000 % 


-0 


00002 


+ 


00000 i 




£r 


-0 


00002 


+ 


00000 % 





00001 


+ 


00000 i 




u l 





00104 


+ 


00000 % 


-0 


00003 


+ 


00000 i 




Ur 


-0 


00001 


+ 


00000 % 





00001 


+ 


00000 i 




d L 


-0 


00001 


+ 


00000 % 


-0 


00005 


+ 


00000 % 




dn 





00000 


+ 


00000 i 





00000 


+ 


00000 i 


b L 


(m t = 100) 


-0 


00002 


+ 


00000 i 


-0 


00005 


+ 


00000 i 


b L 


{mt = 150) 


-0 


00001 


+ 


00000 i 


-0 


00005 


+ 


00000 i 


b L 


(mt = 200) 





00001 


+ 


00000 i 


-0 


00005 


+ 


00000 i 



Table 5 
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m t ( GeV) 


150 


150 


175 


175 


200 


200 


m H ( GeV) 


100 


1000 


100 


1000 


100 


1000 


9z( m z) 


0.55516 


0.55405 


0.55641 


0.55523 


0.55784 


0.55656 


s 2 (m 2 z ) 


0.23119 


0.23245 


0.23040 


0.23170 


0.22952 


0.23086 


5 b (m 2 z ) 


-0.00789 


-0.00792 


-0.00994 


-0.00999 


-0.01226 


-0.01230 


r„ 


166.95 


166.61 


167.32 


166.97 


167.75 


167.37 


r e = r, 


83.81 


83.59 


84.04 


83.80 


84.30 


84.04 


r T 


83.62 


83.40 


83.85 


83.61 


84.11 


83.85 


r„ 


299.20 


297.94 


300.41 


299.09 


301.76 


300.35 


r c 


299.14 


297.88 


300.35 


299.03 


301.70 


300.30 


r d = r s 


382.65 


381.28 


383.77 


382.34 


385.09 


383.56 


r b 


376.90 


375.51 


376.25 


374.79 


375.55 


374.01 


r h 


1740.54 


1733.87 


1744.55 


1737.59 


1749.18 


1741.78 


r z 


2492.63 


2484.27 


2498.44 


2489.70 


2505.15 


2495.82 



Table 6 
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Z parameters measurements 


(section 4.1) 


data 


m z , F z , cr°, -R^, A F ' B , P T , A LR , A FB , A FB , i? 6 


fit parameters 


m z (input), s 2 (m|), g z (m 2 z ), 5 h {m 2 z ) 


d.o.f. = 10-4 


external parameter 


a s 


g\(m\) = 0.5542 - 0.00030 ± 0.0017 / 
s 2 (m|) = 0.2313 + 0.00008 ± 0.0007 p corr = 
6 b {m%) =-0.0061 - 0.00430^-^ ±0.0034 \ 

a // , r n r ( a s ~ 0.1017\ 2 1 / 
xL./(d.o.f.) = [2.48 + ( 00127 )j/6 


' 1 0.14 -0.36 ^ 
1 0.20 

1 J 




Low energy neutral current experiments (section 4.2) 


data 


(fli i 9ri^\i ^r)i {pve , s^ e ), <5ty, (2Ci„— , 2C 2u — C 2 d) 


fit parameters 


s 2 (0), fl|(0) 


d.o.f. = 9-2 


£|(0) = 0.5462 ±0.0036 | 
s 2 (0) = 0.2353 ±0.0044 j 

xLn/(d.cf.) = 2.22/7 


Pcorr — 0.53 


Charged current experiments 


(section 4.3) 


data 




fit parameters 


G F (input), ^(0) 


d.o.f. = 2-2 


^(0) = 0.4225 0.0031 &G °'°° 55 ± 0.0017 
xLn/(d.o.f.) = 0/0 



Table 7 
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data 


no-EW 


IBA 


Exact 


SM 


rot (GeV) 






175 


175 


175 


175 


150 


150 


m H (GeV) 






100 


100 


60 


1000 


60 


1000 


S 






-0.2323 


-0.2323 


-0.2832 


-0.0749 


-0.2638 


-0.0555 


T 






0.8869 


0.8869 


0.9174 


0.5875 


0.6136 


0.3002 


U 






0.3577 


0.3577 


0.3591 


0.3529 


0.2993 


0.2931 


8g 






0.0055 


0.0055 


0.0055 


0.0055 


0.0055 


0.0055 


l/a(ro|) 




128.85 


128.71 


128.71 


128.71 


128.71 


128.72 


128.72 


s 2 (ro|) 




0.2312 


0.2304 


0.2304 


0.2301 


0.2317 


0.2309 


0.2325 






0.5486 


0.5564 


0.5564 


0.5564 


0.5552 


0.5552 


0.5540 


<5 6 (m|) 








-0.0099 


-0.0100 


-0.0100 


-0.0079 


-0.0079 


s 2 (0) 
S%(0) 




0.2388 


0.2389 


0.2389 


0.2386 


0.2401 


0.2394 


0.2408 




0.5486 


0.5492 


0.5492 


0.5493 


0.5480 


0.5481 


0.5468 


9w(0) 




0.4218 


0.4242 


0.4242 


0.4245 


0.4224 


0.4229 


0.4208 


T z (GcV) 


2.489 ± 0.007 


2.481 


2.519 


2.493 


2.494 


2.484 


2.488 


2.479 






2.487 


2.524 


2.498 


2.499 


2.490 


2.493 


2.484 






2.493 


2.530 


2.504 


2.505 


2.495 


2.499 


2.490 


a° h (nb) 


41.56 ± 0.14 


41.53 


41.53 


41.52 


41.52 


41.52 


41.50 


41.51 






41.47 


41.47 


41.46 


41.46 


41.47 


41.45 


41.46 






41.42 


41.42 


41.41 


41.41 


41.42 


41.39 


41.40 


Re 


20.763 ± 0.049 


20.734 


2\).7A7 


20.689 


20.693 


20.665 


20.701 


20.673 






20.801 


20.814 


20.756 


20.760 


20.732 


20.769 


20.741 






20.869 


20.880 


20.823 


20.827 


20.799 


20.836 


20.808 




0.0158 ± 0.0018 


0.0167 


0.0182 


0.0167 


0.0171 


0.0144 


0.0157 


0.0132 


Pr 


-0.139 ± 0.014 


-0.149 


-0.156 


-0.148 


-0.150 


-0.138 


-0.144 


-0.132 


A L r 


0.1637 ± 0.0075 


0.1494 


0.1557 


0.1480 


0.1500 


0.1378 


0.1438 


0.1318 


Rb 


0.2203 ± 0.0027 


0.2183 


0.2182 


0.2157 


0.2156 


0.2157 


0.2165 


0.2165 






0.2183 


0.2182 


0.2157 


0.2157 


0.2157 


0.2165 


0.2166 






0.2183 


0.2182 


0.2157 


0.2157 


0.2157 


0.2165 


0.2166 


^FB 


0.099 ± 0.006 


0.105 


0.109 


0.104 


0.105 


0.096 


0.101 


0.092 




0.105 


0.109 


0.104 


0.105 


0.096 


0.101 


0.092 






0.105 


0.109 


0.104 


0.105 


0.097 


0.101 


0.092 


^FB 


0.075 ± 0.015 


0.075 


0.078 


0.074 


0.075 


0.069 


0.072 


0.065 




0.075 


0.078 


0.074 


0.075 


0.069 


0.072 


0.065 






0.075 


0.078 


0.074 


0.075 


0.069 


0.072 


0.066 




(a, = 0.11) 


7.65 


26.38 


11.16 


11.00 


19.88 


10.78 


29.21 




(a a = 0.12) 


7.40 


35.10 


10.71 


10.94 


16.35 


10.15 


25.10 




(a a = 0.13) 


12.87 


49.38 


15.76 


16.39 


18.31 


15.09 


26.55 


•> 


0.2980 ± 0.0044 


0.2887 


0.2893 


0.2995 


0.2998 


0.2973 


0.2979 


0.2955 


9 

9r 


0.0307 ± 0.0047 


0.0302 


0.0303 


0.0295 


0.0295 


0.0297 


0.0295 


0.0298 


Si 


-0.0589 ± 0.0237 


-0.0588 


-0.0589 


-0.0634 


-0.0634 


-0.0634 


-0.0633 


-0.0632 




0.0206 ± 0.0160 


0.0181 


0.0182 


0.0177 


0.0177 


0.0178 


0.0177 


0.0178 


9 

X 




6.91 


6.09 


0.24 


0.29 


0.25 


0.19 


0.78 


S e// 


0.233 ± 0.008 


0.239 


0.239 


0.230 


0.230 


0.231 


0.231 


0.232 


Pe// 


1.007 ± 0.028 


1.000 


1.001 


1.013 


1.013 


1.011 


1.011 


1.009 


x 2 




0.61 


0.60 


0.18 


0.21 


0.06 


0.11 


0.02 


Qw 


-71.04 ± 1.81 


-74.89 


-74.98 


-73.21 


-73.17 


-73.31 


-73.17 


-73.30 


x 2 




4.52 


4.74 


1.43 


1.39 


1.57 


1.38 


1.57 


2Ci„ — Cid 


0.938 ± 0.264 


0.709 


0.709 


0.723 


0.724 


0.717 


0.720 


0.713 


2C2u — Cid 


-0.659 ± 1.228 


0.081 


0.080 


0.104 


0.105 


0.096 


0.101 


0.092 


x 2 




1.96 


1.94 


1.27 


1.23 


1.51 


1.40 


1.69 


m w 


80.24 ± 0.16 


79.95 


80.39 


80.39 


80.42 


80.22 


80.27 


80.08 


x 2 




3.23 


0.91 


0.91 


1.28 


0.02 


0.03 


1.06 


Xtot 


(a s = 0.11) 


24.87 


40.66 


15.20 


15.40 


23.29 


13.88 


34.33 




(a s = 0.12) 


24.62 


49.38 


14.74 


15.34 


19.76 


13.26 


30.22 




fa s = 0.13) 


30.10 


63.65 


19.79 


20.78 


21.72 


18.20 


31.66 
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Figure captions 



Fig.l The SM predictions for the (S, T, U, 5 b ) parameters defined in eqs. (|2.34 ) and 



( [2. 22 ) are shown as functions of m t for selected m H values. Their closed analytic 



expressions are given in appendix C. a s (m z ) is set to 0.12 in the two- loop 0(aa s ) 
corrections for S, T, U [|4j] and 5 b (m 2 z ) 

Fig. 2 The four charge form factors in the minimal SM as functions of the momentum 
transfer scale. The SM predictions are given for m t = 100, 150, 200 GeV and 
m H = 100, 1000 GeV. The parametrization [^] of the hadronic vacuum polariza- 
tion contribution is used in the space-like region (— m 2 z < q 2 < 0). In the time-like 
region (0 < q 2 < m 2 z ) only the heavy quark (c, b) threshold corrections are taken 
into account. The light quark contributions at \q 2 \ > m 2 z are calculated in per- 
turbative QCD by requiring continuity at q 2 = m z . See appendix B for details. 

Fig. 3 The Z total width T z as function of the universal charge form factors g z {m 2 z ) and 
s 2 (m 2 z ) for S b (m 2 z ) = 0(a), 5 b {m 2 z ) = -0.01(b) and 6 b (m%) = -0.02(c). Three 
cases of a s (m z ) (0.11,0.12,0.13) are shown for each 5 b {m 2 z ). 

Fig. 4 The hadronic cross section on the Z-pole cr°(LEP) as function of s 2 (m 2 z ) and 5 b (m 2 z ) 
for a s (m z ) = 0.11,0.12,0.13. The solid (dashed) lines are obtained for g 2 z {m 2 z ) = 
0.55 (0.57). Here a°(LEP) is defined by 127iT e Th/m 2 z T 2 z in ref. PSfl ; see discussions 
in the text. 

Fig. 5 The ratio Re = cr°(LEP)/(T°(LEP) of the on-Z-pole cross sections as function of 
s 2 (m 2 z ) and Sb(m 2 z ) for a s (m z ) = 0.11,0.12,0.13. The solid (dashed) lines are 
obtained for g z (m 2 z ) = 0.55 (0.57). 

Fig. 6 The ratio R b = u^/a^ of the 6-quark production cross section to the hadronic cross 
section on the Z-pole as function of s 2 (m|) and Eb{m 2 z ) for a s (m z ) = 0.11, 0.12, 0.13. 
The solid (dashed) lines are obtained for g 2 z {m 2 z ) = 0.55 (0.57). 

Fig. 7 The asymmetries on the Z-pole as functions of the effective charge s 2 (m|): the lep- 
tonic forward-backward asymmetry Ap^ (a), the left-right beam polarization asym- 
metry A^ R (b), the 6-quark forward-backward asymmetry Ap^ (c), the forward- 
backward asymmetry of the c-quark A FR (d). The solid lines ('Full') are obtained 
from the full helicity amplitudes (|2.2|) including the 7 and Z exchange as well as the 
box contributions (which are negligibly small). The long dashed lines ('Full —7') 
are obtained by subtracting from the full amplitudes (|2.2|) the real and imaginary 



119 



parts of the 7-exchange contribution QiQj [e 2 (m|) — i e 2 A 77 (m|) ]/s. The thick 
dashed lines ('Z-only') are obtained by retaining only the Z-pole term, the term 
multiplying the Z-propagator factor in eq. (|2.2|) . The dotted lines ('Z-only(IBA)') 
are obtained by using the improved Born approximation to the Z-exchange ampli- 
tudes. The thick dashed lines ('Z-only') are used in the present analysis. QCD 
corrections to A^ B and Ap^ are calculated for a s = 0.12 



Fig. 8 m t dependence of the SM predictions for the electroweak Z boson parameters. Pre- 
dictions for three values of m H are shown by dashed lines (60 GeV), solid lines 
(300 GeV) and by dash-dotted lines (1000 GeV), all calculated for a s {m z ) = 0.12. 
Also shown by straight lines are the mean (dotted lines) and the 1-er allowed ranges 
of the experimental data [pE, 3~I| (see sections 4.1 and 6.2). 

Fig. 9 a s (m z ) dependence of the SM predictions for the electroweak Z boson parameters. 
Predictions for three values of m t are shown by dashed lines (100 GeV), solid lines 
(150 GeV) and by dash-dotted lines (200 GeV), all for m H = 100 GeV. Also shown 
by straight lines are the mean (dotted lines) and the 1-er allowed ranges of the 
experimental data |26| (see section 4.1). 



Fig. 10 Relation between the model- independent parameters of the u^—q scattering experi- 
ments (gl, g%) and the two universal form factors (s 2 (0), gf (0)). The 1-er contour 
of the present data |78| is also shown: see eq. ( |4.15|) in section 4.2.1. The p — 1 line 
corresponds to <?f(0) = 4\/2G^m| = 0.5486. 

Fig. 11 Relation between the v{p)—e scattering cross sections per neutrino energy (er(^e) / E u , 
a(ue)/E u ) and the two universal parameters (s 2 (0), gf (0)). The 1-er contour shows 
the experimental constraint: see eq. ( |4.19| ) p4f and eq. ( f4.20| ) in section 4.2.2. The 



p = 1 line corresponds to (?|(0) = 4 y^GVml = 0.5486. 

Fig. 12 Relation between the model-independent parameters (2Ci u — Cid, 2C 2u — C2d) 

of the e-D polarization asymmetry experiments and the two universal form factors 
(s 2 (0), ^|(0)). The 1-cr contour of the present data |!7j is also shown: see eq. (|4.33|) 
in section 4.2.4. The p = 1 line corresponds to g z (0) = ^\[2Gfw? z = 0.5486. 

Fig. 13 The weak charge (Qw) of the cesium atom 55 3 C S in the atomic parity violation exper- 
iments as function of the two universal parameters (s 2 (0), g z (0)). The 1-cr contour 
of the present data [|80| is shown by dashed lines: see eq. fl4.22|) in section 4.2.3. The 



p = 1 line corresponds to g z (0) = 4 \f2Gp rn 2 z = 0.5486. 
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Fig. 14 The SM predictions for m w as functions of m t and m H for S a — 0, S G — 0.0055 and 
a s = 0.12. The 1-er allowed range of the present data [M is shown by thick dashed 



lines: see eq. ( |4.38|) in section 4.3. 



Fig. 15 3-parameter fit to the Z boson parameters: the Z&l&l vertex form factor 5b(m 2 z ) is 
introduced as the third parameter of the fit in addition to the two universal charge 
form factors s 2 (m 2 z ) and g\{j^?z)'- see eq. ( f4.3|) . The 1-er contours are shown for three 
representative a s {m z ) values, 0.11 (dashed lines), 0.12 (solid lines), 0.13 (dot-dashed 
lines). Also shown are the SM predictions in the range 100 GeV< m t <200 GeV 
and 50 GeV< m H <1000 GeV, which are calculated assuming (Ai) hadrons = —3.88 
(5 a = 0) |27| for the hadronic vacuum polarization contribution to l/a(m z ). 

Fig. 16 Fit to the low energy neutral current data in terms of the two universal charge form 
factors s 2 (0) and g|(0). 1-a contours are shown separately for the v^-q data ]78]| , 
the v^-e data |84|| , the atomic parity violation (APV) data |SD|, and the SLAC 



e-D polarization asymmetry data [jS7| : see eqs. ( [4.17a|) , ( [4.21|) , ft4.23|) and ([4. 34|) 



respectively. The 1-a contour of the combined fit, eq. ( f4.35|) , is shown by the thick 
contour. The straight dashed line shows the 'tree' level prediction of the minimal 
SM: p = g z (0)/(AV2G F m 2 z ) = 1, or g 2 z (0) = A^2G F m 2 z = 0.5486. 



Fig. 17 Fit to the SLAC e-D polarization asymmetry data |87j m terms of the model 



independent parameters 2C± U — Cid and 2C 2u — C 2 d @] of the effective weak Hamil- 



tonian ( p.76| ). Uncertainties due to the sea-quark contributions (a), the longitudinal 



to transverse virtual photon cross section ratio R = cxl/^t (b), and the higher twist 
effects (c) have been examined, and the fit (d) is obtained after taking account of 
all the uncertainties. Shown by the solid lines are the tree-level predictions of the 
SM, and the blobs show the predictions at selected sin 2 9 W values. 

Fig. 18 The running of the charge form factor g z (q 2 ), 4tc /g z (m z ) — 47r/gf (0), as functions 
of m H calculated in the SM for 100 GeV < m t < 200 GeV. The 1-a allowed range 
from the neutral current experiments on the Z-pole and at low energies, eq. QS.lp , 
is also shown for comparison. 

Fig. 19 The running of the charge form factor g z {q 2 ), s 2 (q 2 ) and e 2 (q 2 ) as expected from 
the one-loop contribution of the wino (fermionic partner of the W in the supersym- 
metric SM) to the three neutral current propagators, (a) 4ir /g z (m 2 z ) — 47r/g|(0) ; 
(b) s 2 (m 2 z )/a(m 2 z ) - s 2 (0)/a; (c) 5 a = l/a(m 2 z ) - 128.72. The SM contributions 
are shown for m t = 100, 200 GeV and m H = 60, 1000 GeV. The singularity at 
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^■wino i^Z 

/2 in (a) reflects |g5| the deviation of the Z line-shape from the stan- 
dard Breit-Wigner form that has been assumed both in the experimental fit and in 
our theoretical formula. The 1-er allowed ranges from the neutral current experi- 
ments on the Z-pole and at low energies, eq. Q5.1J) , are also shown for comparison. 
There is no direct measurement of S a . 

Fig. 20 2-parameter fit to the combined low energy neutral current data and the Z param- 
eters. The latter fit ('LEP+SLC') is copied from Fig. |1^ for a s (m z ) = 0.12. The 



low energy combined fit of Fig. [16] has been rescaled to the m z scale by assuming 
the SM running of the two charge form factors, s 2 (q 2 ) and g\[c^\ which depend on 
m t and m H . Uncertainties due to mt and m H in the SM predictions for the running 
of the form factors are illustrated by drawing the results for m t = 100, 200 GeV 
and m H = 60, 1000 GeV in the same figure. The 1-er contour of the combined 
fit, eq. ( |5.5| ), is given by the thick contour, for which the above uncertainties give 
negligible effects. 

Fig. 21 Global fit to the (S, T, U) parameters for three a s values and S a = l/cv(m|) — 
128.72 = and 5 G = 0.0055. Four 1-er contours are obtained for each a s , by using 
m t = 150, 200 GeV and m H = 100, 1000 GeV in evaluating the running of the 
charge form factors: see eq. Q5.8|) for a parametrization of the fit for m t = 175 GeV 
and m H = 100 GeV. The fourth parameter of the fit, the Zb L b L vertex form factor 
5b(m z ), is allowed to take an arbitrary value, free from SM constraints. The SM 
predictions with 5 a = and 5 G = 0.0055 are also given for 100 GeV < m t < 200 GeV 
and 50 GeV <m H < 1000 GeV. 

Fig. 22 Electroweak constraints on (m t , m H ) in the minimal SM, for three selected a s values 
at S a = 0. Dashed lines show the best m t values for a given m H , and the solid 
contours are for \ 2 — Xmin + 1 anc ^ X 2 = Xmm + 4.61. The minimum point of \ 2 is 



marked by "x". The region m H < 63 GeV is excluded by LEP experiments ||101 



Fig. 23 Total \ 2 °f the SM fit to all the electroweak data as functions of m t for m H = 60, 
300, 1000 GeV and a s {m z ) = 0.11, 0.12, 0.13. The uncertainty S a in the hadronic 
vacuum polarization contribution to the effective charge l/cv(m|) is shown for three 
cases, 5 a = —0.1 (a), (b), +0.1 (c). The dotted lines are obtained by using the 
approximate formula (|5.9|). The degree of freedom is 19. 



Fig. 24 Total \ 2 of the SM fit to all the electroweak data as functions of m H for m t = 
120,140,160,180,200 GeV and a s (m z ) = 0.11, 0.12, 0.13. The hadronic vacuum 
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polarization contribution to the effective charge l/a(m 2 z ) is fixed by setting 5 a = 0. 



The dotted lines show our approximation (5.9). The degree of freedom is 19. 



Fig. 25 The universal weak mixing form factor s 2 (m|) as determined from various asymme- 
try measurements on the Z-pole: the lepton (e, /x, r) forward-backward asymme- 
tries |[26|| , the r polarization asymmetry | 26fl , the left-right asymmetry and the 



quark (b,c) forward-backward asymmetries p6f : see eqs. ( |4.4| ), ( p.2|) and the foot- 
noted. Also shown is the deviation (that is, % = ((^ 2 ( m D) — 0.2302)/a(s 2 (m|))) 
for each fit individually, where (s 2 (m|)) and a(s 2 (m 2 z )) denote mean and stan- 
dard deviation of each fit, respectively. At the bottom the above x _vames are 
histogrammed. 



Fig. 26 Impact of the left-right asymmetry data []3TJ] by the SLD collaboration. The band 
(mean (dashed line) and the l-a (solid lines)) represents the constraint from the new 
left-right asymmetry data alone. The constraints from the fits with and without (see 
Figs. [15] and ^T]) the new asymmetry data are shown by thick and thin lines, re- 
spectively, for a s (m z ) = 0.11,0.12,0.13. The Zb L b L vertex form factor Sb(m 2 z ) is 
allowed to take an arbitrary value, free from SM constraints. 

(a) : The fits are shown in the s 2 {m 2 z )-g 2 z (m 2 z ) plane: see eq. ( |6.5|) . The SM predic- 
tions are obtained by assuming 5 a = in the range 100 GeV < m t < 260 GeV and 
1 GeV <m H < 1000 GeV. 

(b) : The fits are shown in the S-T plane, where 5 a = and 5 G = 0.0055 are assumed, 
and m t = 174 GeV and m H = 100 GeV are used to calculate the SM running of the 
charge form factors between q 2 = and q 2 = m|: see eq. ( |6.5|) for parametrization 
of the result. The results are insensitive to the actual (m t , m H ) values in the region 



m t > 100 GeV and m H > 50 GeV (see Figs. [18] and [19]). The SM predictions are 



given in the range 100 GeV < m t < 260 GeV and 50 GeV <m H < 1000 GeV. 

Fig. 27 Electroweak constraints on (m t , m H ) in the minimal SM, including the new left-right 
asymmetry data [^TJ, for three selected a s values a s {m z ) = 0.11, 0.12, 0.13, and for 
(a) 5 a = —0.1, (b) 5 a = 0, and (c) 5 a = 0.1. Dashed lines show the best m t values 
for a given m H , and the solid contours are for \ 2 — Xmin + 1 an d X 2 — Xmin + 4.61. 
The minimum point of \ 2 is marked by "x". The region m H < 63 GeV is excluded 



by LEP experiments [|101 



Fig. 28 Total \ °f the SM fit to all the electroweak data including the new left-right asym- 



metry data as functions of m t for m H =60, 300, 1000 GeV and a s (m z ) =0.11, 



0.12, 0.13. The uncertainty 5 a in the hadronic vacuum polarization contribution to 
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the effective charge l/a(m|) is shown for three cases, 5 a = —0.1 (a), (b), +0.1 
(c). The dotted lines are obtained by the approximate formula Q6.6Q . The degree of 
freedom is 19. 

Fig. 29 Total \ 2 of the SM fit to all the electroweak data including the new left-right asym- 
metry data |31] as functions of m H for m t = 100 — 200 GeV, for three selected a s 
values (a) a s {m z ) = 0.11, (b) a s {m z ) = 0.11 and (c) a s {m z ) = 0.11, at 5 a = 0. 
The dotted lines show our approximation ( |6.6| ) obtained by fitting the \ 2 values in 
the region 63 GeV < m H < 1000 GeV. The degree of freedom is 19. 

Fig. 30 2-parameter fits to the Z boson parameters, where in (a) s 2 (m|) and gz( m z) are 
free parameters, and in (b) S and T are free parameters. In both cases the Zb L b L 
vertex correction is assumed to be dominated by the SM contribution, and the m t 
value in the vertex correction is treated as external parameter in the fit. The 1- 
o contours are shown for three representative a s (m z ) values, 0.11 (dashed lines), 
0.12 (solid lines), 0.13 (dot-dashed lines). Also shown are the SM predictions in 
the range 100 GeV < m t < 200 GeV for 1 GeV < m H < 1000 GeV (a), and for 
50 GeV < m H < 1000 GeV (b). The SM predictions in (a) and the 1-a contours in 
(b) are obtained by assuming 5 a = l/a(m|) — 128.72 = 0. 

Fig. 31 Electroweak constraints on (m t , m H ) in the minimal SM, including the new left-right 
asymmetry data |3l[ and the constraint mt = 174 ± 16 [|106|1 , for three selected a s 



values at 5 a = 0. Dashed lines show the best m t values for a given m H , and the solid 
contours are for \ 2 — Xmin + 1 an d X 2 = Xmm + 4.61. The minimum point of \ 2 is 
marked by "x". The region m H < 63 GeV is excluded by LEP experiments ||101]1 . 

Fig. 32 Constraints on the Higgs mass in the SM from all the electroweak data includ- 



ing the new left-right asymmetry data [31] . Here the top mass m t is considered 



as external parameter with negligible uncertainty. Upper (solid lines) and lower 
(dashed lines) bound of the Higgs mass at 95% CL are shown as functions of m t for 
a s (m z ) = 0.11,0.12,0.13. The hadronic vacuum polarization contribution to the 
effective charge l/a(m|) is set by 5 a = l/«(m|) — 128.72 = 0. 

Fig.33 The contributions to \ 2 from each sector of the analysis in the SM: (a) from the Z 
parameters including the new left-right asymmetry data pH] , (b) from the low energy 
neutral current experiments and (c) the m w measurements. They are calculated 
as functions of m H for mt = 100-200 GeV, at a s (m z ) = 0.12 and 5 Q = 0. The 
degree of freedom is 9 for the Z parameters (a), 9 for the low energy neutral current 
experiments (b), and 1 for m w (c). 



